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APPLICATION  OF  RESEARCH  TO  THE  NEEDS 
OF  THE  U.S.  NAVY 


Mechanical  face  seals  are  used  in  numerous  applications  in 
Naval  machinery.  These  applications  range  from  propeller  shaft 
seals  to  boiler  feed  pump  stcls.  In  such  equipment  the  mechani¬ 
cal  seal  plays  a  vital  role.  When  such  seals  fail,  repair  is 
costly  both  in  terms  of  lost  time  and  direct  costs,  so  any 
improvement  in  seal  life  and  reliability  would  be  of  significant 
benefit. 

As  more  advanced  equipment  is  designed,  it  is  sometimes  dif¬ 
ficult  to  achieve  desired  performance  in  more  severe  service 
environments  with  the  present  state  of  the  art  of  seal  design. 
Thus,  an  improvement  in  seal  technology  would  serve  this  impor¬ 
tant  application. 

One  objective  of  the  research  herein  is  to  further  the 
understanding  of  rtechanical  face  seal  lubrication  phenomena. 
Another  objective  is  to  develop  the  capability  of  designing  con¬ 
tacting  face  seals  having  a  longer  life,  greater  reliability,  and 
for  extreme  environments.  The  immediate  objective  herein  is  to 
apply  the  knowledge  gained  to  the  design  of  a  small  scale  sub¬ 
marine  type  seal.  Thus,  the  objectives  of  this  research  are  com¬ 
patible  with  mechanical  face  seal  needs  for  Naval  machinery. 
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CHAPTER  1 
INTRODUCTION 


Mechanical  Face  Seals 

Applications  of  face  seals  range  from  water  pumps  to  com¬ 
pressors,  to  power  pumps,  to  propeller  shafts.  In  many  applica¬ 
tions  the  reliability  of  the  mechanical  seal  is  of  the  greatest 
importance  to  the  reliability  of  the  equipment  itself. 

Mechanical  face  seal  technology  has  been  steadily  improving 
over  the  past  several  decades.  However  there  still  remain 
demands  for  seal  performance  which  have  not  been  met.  One  such 
application  of  note  is  the  submarine  propeller  shaft. 

Such  demands  on  a  particular  technology  can  often  be  satis¬ 
fied  by  first  improving  the  technology.  In  the  case  of  mechani¬ 
cal  face  seals,  the  main  barrier  to  advancement  has  been  that  the 
mechanics  of  seal  operation  are  not  well  enough  understood  to  be 
able  to  reasonably  anticipate  seal  performance  as  a  function 
design  parameters. 

During  the  past  five  years  of  this  research  program,  much 
has  been  learned  about  controlling  the  hydrostatic  and  hydrody¬ 
namic  mechanisms  which  enhance  face  seal  operation.  In  this 
work,  this  knowledge  is  applied  to  the  design  of  an  improved 
small  scale  submarine  shaft  seal.  Theory,  design,  and  test 
results  are  presented.  The  results  show  promise  that  significant 
improvements  in  submarine  shaft  seal  design  are  possible,  and 
the  application  of  such  designs  could  greatly  increase  the  length 
of  trouble-free  shaft  seal  service. 

Seal  Lubrication 

As  background,  the  mechanical  face  seal  consists  basically 
of  two  annular  rings  which  rotate  relative  to  each  other  and 
which  are  pressed  together  by  spring  and  fluid  pressures  (see 
Figure  1-1).  In  conventional  seals,  the  surfaces  that  rub 
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together  are  generally  manufactured  as  flat  as  possible  initially 
so  as  to  minimize  leakage.  The  effective  gap  between  the  faces 
is  ideally  quite  small  (order  of  1  ym)  so  that  leakage  flow 
across  the  faces  will  be  quite  small.  The  difficulty  in 
designing  a  mechanical  seal  is  in  maintaining  the  gap  at  a  very 
low  value  while  at  the  same  time  providing  a  definite  lubricant 
film  between  the  faces. 

The  load  that  must  be  supported  at  the  faces  of  a  mechanical 
seal  is  due  primarily  to  loading  caused  by  the  sealed  pressure. 
The  load  support  at  the  faces  is  derived  from  fluid  pressure  and 
mechanical  pressure.  If  the  fluid  pressure  at  the  faces  is  large 
enough  to  support  all  of  the  load,  then  there  will  be  no  contact 
and  no  adhesive  wear.*  If  none  of  the  load  is  supported  by 
fluid  pressure,  the  load  must  be  carried  by  mechanical  contact, 
and  the  wear  rate  will  be  large. 

In  practice,  conventional  seals  often  operate  at  one  of  two 
extremes.  At  one  extreme,  a  large  gap  will  be  created  by  hydro¬ 
static  or  hydrodynamic  pressure  or  distortion,  all  of  the  load 
will  be  supported  by  fluid  pressure,  and  the  seal  will  leak  a  lot 
and  wear  very  little.  At  the  opposite  extreme,  the  gap  will 
close  completely.  Leakage  will  be  low  but  only  a  fraction  of  the 
load  will  be  carried  by  fluid  pressure,  and  wear  and  heat  genera¬ 
tion  will  increase. 

Based  on  the  above,  it  can  be  concluded  that  an  effective 
seal  should  operate  between  these  two  extremes--having  both  ade¬ 
quate  fluid  pressure  load  support  and  low  leakage.  The  seal 
should  operate  so  that  it  just  touches  to  minimize  leakage  but 
such  that  the  load  is  carried  by  fluid  pressure.  To  do  this 


♦There  may  still  be  abrasive  or  corrosive  wear  even  if  the  sur¬ 
faces  do  not  touch. 


requires  that  any  fluid  pressure  generation  mechanism  used  to 
provide  load  support  to  the  seal  must  be  very  carefully  con¬ 
trolled.  At  present  in  commercial  seals,  this  is  left  primarily 
to  chance  and  sometimes  seals  operate  at  one  of  the  undesirable 
extremes  mentioned. 

In  this  research  program,  attention  has  been  devoted  toward 
studying  the  effects  of  waviness  as  a  source  of  controlled  hydro- 
dynamic  and  hydrostatic  load  support.  Waviness  was  selected 
because  it  is  controllable.  In  this  present  work,  waviness  has 
been  applied  as  the  basis  for  the  design  of  an  improved  small 
scale  submarine  shaft  seal. 


Background 

ONR-sponsored  research  on  mechanical  seals  had  been  con¬ 
ducted  for  five  years  prior  to  the  beginning  of  the  submarine 
seal  design  phase  described  herein.  The  work  being  reported 
evolved  from  various  discoveries  over  this  five  year  period,  and 
the  past  work  must  be  reviewed  to  better  understand  the  nature  of 
the  current  work. 

As  a  starting  point  for  this  Navy  research  program,  the 
effects  of  waviness  on  seal  performance  were  modeled  in  some 
detail.  In  the  first  annual  report  for  this  project.  Reference 
[1],  this  general  problem  was  solved  using  a  one-dimensional 
theory.  In  the  second  annual  report,  [2],  the  much  more  complex 
two-dimensional  solution  to  the  above  problem  was  solved.  The 
effects  of  waviness,  roughness,  asperity  contact,  wear,  cavita¬ 
tion,  and  elastic  deflection  were  included  in  this  model.  Using 
this  model,  predictions  were  made  for  the  relative  wear  rate, 
friction,  and  leakage  as  a  function  of  roughness,  waviness, 
speed,  size,  pressure,  viscosity,  and  material. 

A  number  of  conclusions  were  reached  based  on  these  first 
two  annual  reports: 


1)  The  effects  of  roughness  on  hydrodynamic  lubrication  are 
not  completely  understood.  Certain  fundamental  ques¬ 
tions  remain  concerning  the  roughness  model  used. 

2)  As  to  the  potential  of  utilizing  hydrodynamic  effects 
caused  by  parallel  face  waviness  to  advantage  by  design, 
the  results  show  that  wear  rate  and  friction  can  be 
greatly  reduced  while  maintaining  leakage  at  acceptable 
levels. 

3)  While  a  comparison  of  predicted  results  to  experimental 
results  given  in  the  literature  is  generally  good,  data 
contained  in  the  literature  is  incomplete,  so  more  com¬ 
plete  experimental  data  are  needed  for  comparison. 

4)  In  low  viscosity  or  heavily  loaded  applications  where 
some  touching  is  expected  to  occur,  waviness  will  wear 
away  with  time  and  any  benefit  derived  will  be  lost 
unless  something  is  done  to  counteract  this  effect. 

5)  Based  upon  data  for  some  commercial  seals  and  using  the 
model,  it  was  determined  that  there  was  insufficient 
accidentally  caused  waviness  to  produce  significant 
hydrodynamic  effects  in  water.  One  cannot  generalize  to 
say  that  such  effects  do  not  occur  in  commerical  seals. 
However,  using  the  model  the  question  can  be  answered  on 
a  case  by  case  basis. 

Item  1)  was  treated  extensively  in  the  third  annual  report 
[3].  Even  after  this  analysis  certain  fundamental  questions 
remain  concerning  how  to  deal  with  roughness  in  lubrication  prob¬ 
lems.  However,  this  thorough  analysis  led  to  conclusions  allow¬ 
ing  certain  simplifying  assumptions  discussed  in  the  fourth 
annual  report  [4], 

Item  2)  was  also  treated  extensively  in  the  third  annual 
report  [3].  A  methodology  for  the  design  of  a  wavy  face  seal  was 
developed  and  applied.  Theoretical  results  showed  large  reduc¬ 
tion  in  friction  and  wear  rate  compared  to  conventional  designs 
whereas  leakage  could  be  controlled. 


Concerning  Item  3),  the  second  and  third  annual  reports 
[2,31  describe  a  test  apparatus  designed  to  test  the  wavy  seal 
theory.  This  apparatus  has  been  in  operation  for  more  than  three 
years  and  many  tests  have  been  conducted.  These  test  results  are 
reported  in  the  fourth  and  fifth  annual  reports  [4,5]. 

Early  in  the  test  program  it  was  observed  that  the  type  of 
waviness  which  can  be  practically  applied  is  not  of  the  radially 
parallel  type.  Waviness  generally  consists  of  alternating  tilt 
plus  radially  parallel  waviness.  Based  on  these  considerations, 
a  new  model  for  predicting  performance  was  developed  a^d  appears 
in  the  fourth  annual  report  [4], 

Concerning  Item  4)  above,  a  solution  to  this  problem  was 
first  proposed  in  the  first  annual  report  [1],  It  was  proposed 
to  move  the  waviness  slowly  around  the  seal  so  that  whatever  wear 
occurred  would  be  uniformly  distributed.  Then  the  shape  of  the 
wave  would  be  preserved  and  tests  using  a  constant  wave  could  be 
made.  The  concept  is  illustrated  in  Figure  1-2.  This  concept 
was  incorporated  into  the  test  apparatus  and  is  described  in 
detail  ir,  References  [2]  and  [3].  In  the  fifth  annual  report 
[5],  a  new  concept  to  move  the  waviness  with  no  internal  moving 
parts  is  described  in  detail. 

In  the  fifth  annual  report  [5],  additional  experimental 
results  using  waviness  are  presented,  the  wavy  seal  model  is 
further  improved,  and  theory  and  experiment  are  compared.  The 
results  using  a  new  concept,  that  of  a  self-generating  seal  pro¬ 
file,  are  reported.  Results  from  tests  of  the  effects  of  radial 
taper  and  high  temperature  environment  are  also  reported  and  com¬ 
pared  to  theory. 

In  summary,  both  the  theoretical  and  experimental  basis  for 
applying  waviness  to  a  mechanical  seal  to  reduce  friction  and 
wear  were  well  established  during  these  first  five  years.  In  the 
present  work,  the  wav i ness  concept  is  applied  to  the  design  of  a 
long  life  submarine  shaft  seal  as  described  later. 
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Figure  1-2.  Moving  Waviness  Concept. 
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The  results  developed  to  date  and  discussed  in  the  five 
annual  reports  are  also  presented  in  several  papers  and  theses 
[6-24]. 


Wavy  Face  Seal 

The  concept  of  waviness  is  that  the  film  thickness  varies  in 
some  fashion  circumferenti al ly  around  the  seal.  Generally  speak¬ 
ing,  film  thickness  may  vary  radially  as  well  as  tangentially. 

h  =  h(r,  9)  .  (1-1) 

In  the  present  work  interest  is  focused  upon  film  thickness 
shapes  of  the  following  functional  form 

h  =  hg  +  f (r)  cos  ne  .  (1-2) 

At  any  particular  radius  r  the  film  shape  is  periodic  with  n 
waves  around  the  seal  and  is  therefore  wavy.  However,  film  shape 
can  also  vary  in  some  general  manner  with  r. 

If  f(r)  =  const  then  the  faces  are  always  radially  paral¬ 
lel.  This  component  of  film  thickness  variation  is  commonly 
termed  waviness.  If  f(r)  *  const,  then  the  faces  are  not  in  gen¬ 
eral  radially  parallel.  f(r)  is  referred  to  as  tilt.  Thus,  the 
film  thickness  shapes  of  interest  are  combinations  of  waviness 
and  tilt.  Since  at  any  radius  the  film  thickness  is  wavy,  the 
combination  of  waviness  and  tilt  defined  above  will  also  be 
called  a  wavy  film  shape. 

The  reason  for  choosing  film  shapes  as  described  by  Equation 
(1-2)  as  a  subject  for  study  is  that  these  shapes  can  conven¬ 
iently  be  generated  by  planned  mechanical  distortions  in  a  seal 
ring,  and  the  shapes  also  include  common  modes  of  unplanned  dis¬ 
tortion  found  in  operating  seals.  For  example,  generally  seals 
undergo  a  uniform  tilt  due  to  pressure  and  thermal  deformation. 
When  rings  are  loaded  by  any  nonaxisymmetr ical  load  they  become 
wavy  as  described  by  Equation  (1-2). 


Now,  for  the  sake  of  illustration,  assume  a  seal  has  a  wavy 
film  thickness  shape  given  by  Equation  (1-2)  where  f(r)  is  a  lin¬ 
ear  function  of  r.  This  gives  a  film  shape  as  shown  in  Figure 
1-3  where  for  each  of  the  n  periods  the  seal  touches  all  across 
at  one  point  and  is  radially  convergent  n/n  radians  away.  At  any 
radius,  the  seal  is  wavy  circumferential ly.  The  waviness 
enhances  hydrodynamic  effects  and  the  radial  convergence  enhances 
hydrostatic  effects. 

It  is  this  shape  which  previous  research  [1-5]  has  shown  has 
the  greatest  potential  to  reduce  friction  and  wear  while  holding 
leakage  at  very  low  levels.  This  is  the  film  shape  used  for  the 
design  of  the  small  scale  submarine  seal. 

Submarine  Seal 

The  general  objectives  of  the  present  three-year  contract 
phase  are: 

1)  to  conduct  further  experiments  using  moving  waviness  to 
further  the  understanding  of  this  concept. 

2)  to  refine  mathematical  models  already  developed  so  as  to 
be  able  to  better  predict  performance. 

3)  to  demonstrate  in  a  practical  way  the  use  of  the  concept 
of  moving  waviness  to  reduce  friction  and  wear  in 
mechanical  seals. 

4)  to  undertake  the  development  needed  to  be  able  to  demon¬ 
strate  the  applicability  of  the  concept  to  submarine 
shaft  seals  and  to  design  a  full  scale  long  life  seal. 

5)  to  explore  other  uses  of  the  waviness  concept  such  as 
for  gas  seals  and  to  create  better  methods  by  which  the 
concept  can  be  applied. 

From  these  objectives,  the  following  specific  tasks  have 
been  defined. 

1)  Design,  fabricate,  and  test  a  nine  wave  optimum  seal. 

As  shown  in  the  previous  annual  report  [5],  nine  waves 
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Figure  1-3.  Mechanical  Seal  with  Waviness  and  Tilt 


are  needed  to  minimize  changes  in  tilt  with  changes  in 
pressure  and  speed. 

2)  Operate  the  seal  for  an  extended  period  of  time  under 
simulated  submarine  operating  conditions  including 
seawater  and  changing  speed,  pressure,  and  alignment. 

3)  Evaluate  compliancy  of  existing  and  proposed  submarine 
seals. 

4)  Conduct  friction  and  wear  tests  on  carbon  and  hard  face 
materials. 

5)  Seek  alternate  methods  of  applying  moving  waviness. 

6)  Perform  analysis  of  a  moving  wave  gas  seal. 

Progress  has  been  made  on  all  of  thes.e  tasks,  and  the  report 
to  follow  details  the  research  efforts  and  achievements  for  the 
period  December  1,  1980  through  June  30,  1982. 
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CHAPTER  2 

EXPERIMENTAL  RESULTS 


A  nine-wave  seal  has  been  designed,  fabricated,  and  tested 
during  this  reporting  period.  The  test  apparatus  has  been  modi¬ 
fied  to  automatically  change  operating  conditions,  and  a  75  per¬ 
cent  balance  ratio  test  has  been  completed.  In  this  chapter 
these  results  will  be  reported. 

Nine-Wave  Seal 


Design--The  first  item  to  be  reported  is  the  design  of  the 
nine-wave  seal  (the  design  procedure  is  discussed  in  Chapter  3; 
the  design  itself  is  discussed  here).  Figure  2-1  shows  the 
assembly  drawing  of  the  nine-wave  seal.  Starting  at  the  left, 
one  of  three  sinusoidally  varying  pressures  generated  by  the 
waviness  drive  unit  [5]  is  directed  into  two  (180°  apart)  of  six 
pressure  channels  (T) in  the  waviness  cylinder  (2).  The  pressur¬ 
ized  fluid  then  passes  through  the  pressure  coupler  (3)  and  is 
ported  to  one  of  three  circumferential  channels  on  the  left  end 
of  the  waviness  adapter  (4).  Eighteen  smaller  passages,  con¬ 
nected  to  each  of  the  tnree  circumferential  channels,  terminate 
at  pressure  pockets  (?)  located  on  the  inside  and  outside  diam¬ 
eters  of  the  waviness  adapter  (4)  (right-hand  end).  The  54  pres¬ 
sure  pockets  (three  sets  of  eighteen)  apply  a  waviness  pressure 
to  the  54  "fingers"  of  the  nine-wave  seal  (?) . 

Figure  2-2  shows  a  modified  isometric  cross  section  of  the 
nine-wave  seal.  The  pads  shown  are  labeled  with  a  (?)  ,  (?)  or  (?) 
depending  on  the  termination  point  of  the  three  sinusoidally 
varying  pressures,  i.e.,  one  pressure,  say,  Plf  exerts  a  wave 
producing  load  on  all  pads  labeled  with  a  (T) .  The  load  on  all 
outer  pads  is  directed  outward  on  the  inside  radius  and  the  load 
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Figure  2-2.  Nine-Wave  Seal  -  Modified  Isometric  Cross 


on  all  inner  pads  is  directed  inward  on  the  outside  radius  of  the 
corresponding  pad.  With  this  configuration  three  sets  of  nine 
waves  each  can  be  imposed  on  the  face  of  the  carbon  insert.  As 
discussed  previously  [5],  by  sinusoidally  varying  the  three 
separate  pressures  (and  waves)  with  time,  the  result  is  one  set 
of  nine  waves  which  moves  circumferentially  around  the  seal  with 
time. 

The  irregular  shape  of  the  nine-wave  seal  (cross  section)  is 
due  in  part  to  the  adjustments  made  on  the  geometry  so  as  to 
obtain  a  zero  moment  design,  i.e.,  there  is  no  rotation  of  the 
seal  about  a  circumferential  axis  due  to  pressure  variation. 

Other  reasons  include  the  need  for  a  low  stiffness,  so  that  the 
wave  could  more  easily  be  transmitted  to  the  carbon  insert  (7), 
and  the  need  for  the  centroid  to  be  placed  at  such  a  location  so 
that  the  wave  contact  point,  between  the  carbon  insert  (7)  and 
the  silicon  carbide  seat  (D ,  is  radially  located  for  optimum 
seal  performance. 

The  secondary  seal  (ll)  is  located  at  the  left  end  of  the 
nine-wave  seal  on  the  inside  diameter.  Springs  to  the  left  of 
the  seal,  housed  in  the  seal  ring  spring  retainer  (0),  provide 
preload  through  the  spring  seat  (^2)  .  The  balance  ratio  for  the 
design  shown  is  unity. 

The  carbon  insert  (7)  is  a  Pure  Carbon  P658RC  material  which 
is  epoxied  into  the  nine-wave  seal  using  3M  1838  B/A  adhesive. 

It  should  be  pointed  out  at  this  time  that  the  carbon  insert,  as 
shown  in  Figures  2-1  and  2-2,  has  a  small  lip  on  the  inside  rad¬ 
ius.  This  lip  has  a  slight  relief  ground  into  it,  approximately 
76  un  deep  in  the  axial  direction  and  500  pm  deep  in  the  radial 
direction,  for  purposes  of  wear  measurements. 

The  drive  ring  (lo)  is  designed  so  that  the  primary  ring  can 
float  and  align  itself  but  cannot  rotate  about  the  shaft  axis. 

The  nine-wave  seal  therefore  takes  its  alignment  from  the  face  of 
the  rotating  secondary  ring.  The  secondary  ring  is  a  Carborundum 
KT®  silicon  carbide  material  which  is  also  of  zero  moment  design. 


Mod  if ic  at ions--Mod ifications  to  the  overall  design  were  made 
in  light  of  the  results  obtained  from  preliminary  testing.  Fig¬ 
ure  2-3  shows  the  original  design  cross  section  of  the  nine-wave 
seal.  The  original  cross  section  shape,  as  mentioned  earlier, 
was  designed  to  meet  specific  requirements,  and  as  a  result,  does 
not  adhere  closely  to  the  assumptions  used  for  the  ring  equations 
used  to  predict  waviness.  Bench  tests  to  measure  waviness  on  the 
carbon  face  showed  that  less  than  half  of  the  needed  waviness  was 
present.  Torsional  deformation  is  predominant  and  for  this  par¬ 
ticular  shape  and  mode  of  deflection,  warping  becomes  of  consid¬ 
erable  importance  and  as  a  result  increases  the  stiffness.  The 
addition  of  a  warping  factor  to  the  design  analysis  did  indeed 
show  that  the  original  cross  section  was  much  too  stiff.  To 
reduce  the  stiffness  and  maintain  a  zero  moment  design  the  cuts 
as  shown  in  Figure  2-4  were  made  to  the  seal. 

Another  problem  encountered  was  that  of  the  epoxy  bond 
between  the  carbon  insert  and  the  seal.  The  original  design 
called  for  a  0.0254  cm  (0.01  in.)  bond  thickness  around  the  car¬ 
bon  insert.  Static  tests  done  with  100  percent  sealed  pressure 
and  with  a  convergent  taper  (400  ym/m)  lapped  into  the  carbon 
face  showed  that  the  epoxy  bond  "creeped"  in  24  hours.  Final 
traces  of  the  radial  profile  showed  no  taper  at  all.  As  dis¬ 
cussed  in  detail  later,  several  bonding  techniques  were  experi¬ 
mented  with.  These  included  the  use  of  Loctite  Speedbonder  319, 
with  and  without  a  primer,  Loctite  Superbonder  420  with  the  addi¬ 
tion  of  two  stainless  steel  rings  for  the  inner  and  outer  diam¬ 
eter  of  the  carbon  insert  to  give  a  near  zero  clearance  fit,  and 
the  use  of  the  3M  1338  B/A  epoxy  also  with  the  same  stainless 
steel  rings.  Static  and  dynamic  tests  on  the  seal  under  water 
pressure  and  independent  tests  on  the  strength  of  each  adhesive 
showed  that  the  3M  1338  8/A  epoxy  with  the  stainless  steel  spacer 
rings  was  the  best  method  of  bonding. 
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By  lowering  the  stiffness  another  problem  developed  and 
hence  a  design  change  was  needed.  Figure  2-5  shows  the  original 
0-ring  configuration.  By  lowering  the  stiffness,  the  "fingers" 
were  subject  to  greater  deflections  which  resulted  in  the  0-rings 
extruding  out  from  under  the  pad  sections.  Even  with  recommended 
clearance,  0-ring  extrusion  was  discovered  to  be  a  problem  simply 
because  the  cyclic  motion  of  the  clearance  proved  to  aid  extru¬ 
sion  by  somewhat  of  a  ratchet  mechanism.  To  remedy  this  problem, 
the  waviness  rings  were  ground  down  in  thickness  to  incorporate  a 
spacer,  and  a  piston  pad  with  a  smaller  0-ring  cross  section  were 
added  (Figure  2-6).  This  now  allowed  the  "fingers"  to  flex  and 
the  0-rings  to  be  contained  without  leakage  and  without  extrud¬ 
ing. 


Waviness  Measurements--As  an  aid  to  understanding  the  prob¬ 
lem  with  creating  sufficient  waviness,  a  means  of  measuring  wavi¬ 
ness  outside  the  seal  environment  was  devised.  A  fixture  was 
made  which  can  pressurize  one  set  of  18  pockets  at  one  time  using 
gas  pressure.  With  one  set  of  pockets  under  pressure  the  seal 
ring  was  placed  on  a  precision  rotary  table.  The  table  rotates 
and  the  stylus  of  the  surface  analyzer  was  placed  so  as  to  read 
axial  displacement  near  the  0-0  of  the  carbon.  The  signal  was 
digitized  and  a  Fourier  analysis  of  waviness  components  was 
made.  Using  this  techique,  the  nine-wave  component  was  easily 
picked  out.  Table  2-1  shows  the  results  of  the  many  measurements 
taken  to  analyze  the  waviness  problem. 

Test  Results--A  series  of  six  seal  tests  were  run  with  the 
new  nine-wave  seal.  The  tests  occurred  at  different  points  in 
the  modification  process  mentioned  before.  As  a  result,  the 
tests  must  be  carefully  compared  to  one  another.  Only  the 
performance  of  Test  No.  117  can  be  evaluated  for  comparison  to 
theory  as  a  final  design  seal. 
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Figure  2-6. 


Final  0-Ring  Configuration  with  Piston  Pads. 
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Table  2-1 

Nine-Wave  Amplitude  Study 
Waviness  Amplitude  Near  Seal  0.0. 


Date 

h 

9 

P 

test 

h  * 

9 

@1500  psi 

Conditions 

1 

6/28 

33 

750 

66 

Close  gap  epoxy  bond,  54-0.030  pads 

6/28 

59 

1390 

64 

Close  yap  epoxy  bond,  54-0.030  pads 

6/15 

80 

1150 

104 

Close  gap  epoxy  bond,  18  0-rings  only 

‘■m 

6/16 

65 

1050 

93 

Close  gap  epoxy  bond,  18  0-rings  only 

6/4 

123 

1175 

157 

Groove  bottom,  18  0-rings 

6/4 

129 

1200 

161 

Groove  bottom,  18  0-rings 

6/2 

29 

1500 

29 

Super  Glue  -  after  test  run, 

90  Durometer 

6/2 

23 

1500 

23 

Super  Glue  -  after  test 

70  Durometer  0-rings 

I 

4/7 

72 

1500 

72 

Carbon-thick  bond  epoxy  after  final 
x-section  modification 

4/6 

73 

1500 

73 

Carbon-thick  bond  epoxy  after  final 
x-section  modification 

• 

3/27 

54 

1500 

54 

Carbon-thick  bond  epoxy  before  final 
x-section  change 

3/15 

16 

1500 

16 

After  first  test  -  thick  epoxy 

• 

*Two 

t  imes 

the  effective  waviness  pressure  in  operation. 

« 
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Table  2-2  shows  the  results  of  the  tests  performed.  Test 
No.  112  (Figure  2-7),  first  of  the  series,  was  run  where  the 
carbon  insert  had  a  bond  line  of  0.0254  cm  3M  1836  B/A  adhesive. 
The  initial  torque  levels  were  quite  high,  approximately  14  N*m. 
There  was  one  shutdown  at  16  hours  into  the  test  due  to  excessive 
torque  levels.  The  performance  of  this  test  was  not  as  expected. 

It  was  at  this  point  that  the  stiffness  of  the  seal  ring  was 
re-analyzed  and  as  a  result,  a  warping  factor  was  incorporated 
into  the  design  computer  program.  Computer  results  showed  that 
to  obtain  the  desired  waviness,  the  stiffness  would  have  to  be 
reduced  by  at  least  50  percent.  The  seal  was  then  modified  and 
Test  No.  113  (Figure  2-8)  was  run.  This  test  had  the  identical 
start-up  behavior  of  Test  No.  112.  The  test  was  stopped  after 
only  4.7  hours. 

The  seal  was  removed  from  the  test  apparatus,  traces  made  of 
the  face,  and  then  a  300  ym/m  convergent  taper  was  lapped  into 
the  face.  The  seal  was  then  replaced  in  the  test  apparatus  and 
pressurized  to  100  percent  pressure  for  a  static  test.  The  seal 
was  left  overnight  under  these  conditions.  After  disassembly, 
traces  showed  no  convergent  taper  at  all,  but  rather  a  flat  face 
condition.  This  result  showed  that  the  epoxy  bond  creeps  and 
allows  the  seal  face  to  rotate  under  pressure.  The  seal  was  then 
re-installed  with  this  flat  face  condition  and  run  without  wavi¬ 
ness  pressure.  The  results,  Test  No.  114  (Figure  2-9),  show  that 
the  performance  was  quite  similar  to  the  two  previous  tests. 
Clearly,  any  waviness  imparted  to  the  seal  was  not  being  trans¬ 
mitted  to  the  carbon  insert.  Other  adhesives  were  investigated 
and  tested  to  try  and  remedy  this  problem. 

The  carbon  insert  was  machined  out  and  a  new  one  installed 
using  Loctite  Speedbonder  No.  319  as  the  adhesive.  The  bond 
thickness  on  the  back  of  the  carbon  was  reduced  to  zero  to 


stiffen  the  connection.  Test  No.  115  (Figure  2-10)  was  then 
run.  It  showed  a  torque  level  which  was  about  half  of  that  of 
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Table  2-2 
Waviness  Tests 

Ph20  =  100%  pressure,  1800  rpm 


p  * 

Wave 
(MPa) 

Test 

No. 

Test 

Ourat  ion 
(h) 

Average  Face 
Temperature 

rc) 

Average 

Torque 

(N*m) 

Average 
Leakage 
(cm3 /min 

5.2 

112 

23.5 

39.6 

11.0 

0.00 

5.2 

113 

4.7 

39.8 

10.2 

0.00 

0 

114 

13.2 

41.1 

11.7 

0.00 

5.2 

115 

67.7 

39.2 

6.4 

0.00 

5.2 

116 

47.8 

38.8 

4.9 

0.00 

5.2 

117 

99.8 

38.4 

2.6 

0.19 

♦Equivalent  waviness  pressure 


Test  #112  -  Nine-Wave  Seal,  100%  pressure,  1800  rpm. 


rpm. 


rpm. 


previous  tests.  However,  there  was  no  leakage  and  some  was  to  be 
expected.  The  test  was  ended  after  47  hours  and  the  carbon 
insert  was  again  machined  out.  There  were  indications  that 
incomplete  bonding  had  occurred  between  the  seal  and  insert  at 
certain  locations,  mainly  on  the  outside  diameter. 

The  next  step  was  to  reduce  the  inner  and  outer  diametrial 
clearances  around  the  carbon  insert.  This  was  accomplished  by 
making  two  stainless  steel  spacer  rings  for  the  inner  and  outer 
diameters.  The  rings  were  0.0254  cm  thick  and  as  a  result  made 
the  carbon  insert  a  slight  press  fit  into  the  seal.  Because  of 
the  press  fit  another  type  of  adhesive  was  needed.  Loctite 
Superbonder  No.  420  was  used  because  of  its  low  viscosity  proper¬ 
ties  and  post  assembly  application  ability.  The  seal  was  bench 
tested  for  waviness  and  then  re-installed  in  the  test  apparatus. 
Test  No.  116  (Figure  2-11)  was  run  for  more  than  47  hours  and 
showed  a  performance  similar  to  Test  No.  115.  Some  leakage  had 
occurred  between  4-1/2  hours  and  14-1/2  hours  into  the  test  but 
at  no  other  time. 

The  test  was  stopped  after  the  torque  readings  became 
erratic,  varying  as  much  as  5  N«m.  The  carbon  ring  was  machined 
out  and  again  showed  indications  of  incomplete  bonding,  this  pos¬ 
sibly  being  the  result  of  water  weakening  the  adhesive  strength. 
It  was  decided  to  test  the  strengths  of  the  various  adhesives 
used  by  independent  experiment.  Table  2-3  shows  the  results. 

The  results  show  clearly  that  the  Loctite  Superbonder  No.  420  is 
weakened  by  water  contamination  and  that  3M  1838  B/A  was  the 
strongest  even  after  soaking  in  water. 

Based  on  the  previous  findings  a  new  carbon  insert  was 
cemented  into  the  seal  along  with  the  stainless  steel  spacer 
rings  using  the  3M  1838  B/A  adhesive.  Because  of  the  high 
viscosity  of  the  epoxy,  installation  required  careful  application 
of  the  expoxy  and  pressing  in  the  carbon.  The  0-rings  in  the 
pressure  adapter  were  also  modified.  Because  of  the  reduced 
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Table  2-3 


Carbon  to  Stainless  Steel  Bond  Strength  Tests 


Test 

Specimen 


Adhesive 

Load  at 
Failure 
(lb) 

Comments 

Loctite  420 
(Super  bonder) 

3.45 

Under  static  line  water  pressure 
for  2  days--clean  break  at  glue 
line  (one  side) 

Loctite  420 
(Super  bonder) 

6.85 

Under  water  0  psi  for  2  days-- 
clean  break  at  glue  line  (one 
side) 

Loctite  420 
(Super  bonder) 

10.25 

Moisture  f ree--c lean  break  at 
glue  line  (both  sides) 

Loctite  420 

— 

Not  tested 

Loctite  319 
(Speed  bonder) 

-- 

Carbon  too  porous  for  adhesion-- 
no  bond 

Loctite  319 
(Speed  bonder) 

— 

Same  as  above 

3M  1838  8/ A 
Epoxy 

29+ 

No  break 

3M  1838  B/A 
Epoxy 

29.05 

Carbon  broke  (not  at  glue  line) 

1/16  in.  — 
Stainless  Steel 


stiffness,  O-ring  extrusion  was  becoming  a  problem,  so  modifica¬ 
tions's  were  made  to  use  smaller  0-rings  along  with  piston  pads 
and  Delrin  retaining  rings.  The  seal  assembly  was  bench  tested 
for  waviness  and  showed  good  results. 

The  seal  was  re-installed  and  Test  No.  117  (Figure  2-12) 
started.  The  test  showed  a  start-up  performance  similar  to  Test 
Nos.  116  and  115.  The  torque  was  falling  to  around  5  N*m  when 
water  contamination  from  the  building's  air  compressor  clogged 
the  hydraulic  pump  to  the  waviness  generator.  The  test  had  to  be 
stopped  after  four  hours  of  operation  and  the  pump  cleaned.  A 
temporary  separate  air  compressor  to  operate  all  air-assisted 
devices  on  the  test  rig  was  added.  The  test  was  re-started  and 
ran  for  almost  100  hours.  Leakage  began  about  eight  hours  after 
start  and  continued  throughout  the  test  at  about  0.2  cm3/min. 

Test  117  shows  torque  and  leakage  values  approaching  those 
expected  for  the  nine-wave  seal  design  and  would  seem  to  indicate 
that  the  design  is  successful  on  a  preliminary  basis.  However, 
radial  taper  measurements  after  the  test  show  that  a  much  greater 
than  expected  radial  taper  is  present  in  the  seal.  This  question 
is  being  investigated  at  this  time. 

Test  Machine  Modification 

Automatic  operation  of  the  test  apparatus  under  variable 
operating  conditions  required  certain  modifications  of  the  test 
apparatus  and  control  hardware  and  software.  These  changes  are 
described  here. 

Speed  Controller--The  test  apparatus  is  equipped  with  a 
3700  W  (5  Hp)  belt  type  variable  speed  drive  ac  motor.  Speed 
change,  as  equipped  from  the  factory,  was  accomplished  by  means 
of  a  hand  crank  which  varied  belt  pulley  ratios.  By  this  method, 
speeds  from  380  to  4000  rpm  were  obtained.  However,  for  the  sim¬ 
ulated  submarine  operation  test  program,  it  was  desired  that 


rpm. 


the  computer  system  would  control  the  speed,  thus  providing  an 
automatic  operation  so  the  operator  could  save  time  and  do  other 
things  while  running  a  test. 

Modification  began  with  the  addition  of  a  Superior  Electric 
Co.  Slo-Syn  stepping  motor.  Model  No.  SS  250-1027,  and  a  2:1 
step-down  gear  ratio  in  place  of  the  hand  crank  mechanism.  This, 
coupled  with  a  Slo-Syn  ST  1800B  stepping  motor  translator  and  a 
H.P.  69335A  stepping  motor  control  card  for  use  in  the  H.P.  6940B 
multiprogrammer  unit,  was  the  basic  hardware  of  the  speed  control 
apparatus.  Since  the  stepping  motor  operates  by  the  input  of  a 
pre-programmed  pulse  train,  each  pulse  of  which  will  rotate  the 
stepping  motor  shaft  by  1.8°,  a  relationship  between  stepping 
motor  shaft  rotation  and  the  main  drive  motor  rpm  was  found. 

Motor  speed  control  was  achieved  by  input ing  the  desired  rpm  into 
the  computer  which  then  converts  this  to  the  appropriate  number 
of  pulses.  A  data  word  is  sent  to  the  stepping  motor  control 
card  in  the  mult i programmer  unit  which  contains  the  desired  num¬ 
ber  of  pulses.  The  stepping  motor  control  card  in  turn  outputs 
the  pulse  train  to  the  translator  wnich  then  sends  the  stepping 
motor  the  required  number  of  pulses  which  in  turn  changes  the 
ratios  of  the  variable  speed  drive. 

Two  other  factors  had  to  be  considered  for  complete  motor 
control.  The  first  is  the  capability  of  forward  and  reverse  con¬ 
trol  and  the  second  is  minimum  speed.  Forward  and  reverse  con¬ 
trol  was  accomplished  by  use  of  two  reed  relays  on  the  relay  card 
in  the  multi  prog rammer  unit,  two  DPDT  relays  and  two  motor 
starter  control  units.  Figure  2-13  shows  the  schematic  diagram 
of  the  directional  control  circuit.  Directional  control  is 
achieved  by  sending  a  control  word  from  the  computer  to  the  relay 
card  to  close  either  Relay  Nos.  9  or  10  depending  on  the  desired 
motor  direction.  A  desirable  feature  of  this  circuit  is  that  it 
is  impossible  to  energize  both  forward  and  reverse  motor  starter 
units  simultaneously  by  closing  both  Relay  Nos.  9  and  10  on  the 


relay  card.  A  rotary  type  switch  was  added  so  as  to  obtain 
manual  directional  control  external  to  the  computer. 

A  lower  minimum  speed  was  obtained  by  the  use  of  a  Reliance 
Electric  variable  frequency  ac  motor  speed  controller.  This  unit 
has  the  capability  of  changing  the  input  frequency  to  the  ac 
motor  from  6  to  60  Hz,  thereby  changing  the  speed  of  the  motor 
also.  This  device  used  in  conjunction  with  the  stepping  motor 
speed  controller  now  gives  a  range  of  speed  control  from  approxi¬ 
mately  40  to  4000  rpm. 

Pressure  Controller--Another  area  of  computer  control  needed 
for  the  test  program  is  that  of  sealed  water  pressure  regula¬ 
tion.  The  test  apparatus  is  equipped  with  an  accumulator  to 
dampen  out  pressure  surges.  The  accumulator  is  shown  in  Figure 
2-14.  When  the  accumulator  is  charged  with  nitrogen  gas  to  75 
percent  of  the  water  pressure,  then  proper  surge  control  is 
obtained  in  that  75  percent  of  the  original  accumulator  volume 
remains  as  gas  when  the  liquid  and  gas  pressures  become  equal. 

But  once  charged,  if  the  water  pressure  is  changed  significantly, 
up  or  down,  then  the  diaphragm  will  also  change  location  and  as  a 
result  poor  performance  of  the  accumulator  can  be  expected.  To 
operate  at  variable  pressure  conditions,  it  therefore  becomes 
necessary  to  be  able  to  detect  the  location  of  the  diaphragm  and 
correct  its  position  when  needed  for  optimum  performance. 

Figure  2-14  shows  how  this  positioning  capability  is 
obtained.  When  the  diaphragm  is  displaced  either  up  or  down  due 
to  changes  in  water  pressure,  the  guide  rod  will  move  the  perman¬ 
ent  magnet  also.  When  the  permanent  magnet  is  moved  so  that  it 
is  adjacent  to  one  of  tiie  magnetic  reed  switches,  the  magnetic 
fielc  induces  a  contact  closure  in  that  reed  switch.  By  correct¬ 
ly  oositioning  the  reed  switches,  an  upper  and  lower  bound  for 
diaphragm  location  can  be  set  along  with  a  sufficient  dead  zone 
of  operation  and  an  optimum  diaphragm  position. 


The  accumulator  diaphragm  positioning  circuit  is  shown  in 
Figure  2-15.  When  the  permanent  magnet  is  adjacent  to  the  No.  3 
reed  switch  then  the  engaging  coil  of  both  Nos.  6  and  7  dual  coil 
relays  will  be  energized,  thereby  setting  the  switch  in  both 
relays  to  the  engaged  position.  If  the  diaphragm  rises  so  that 
it  is  next  to  either  the  No.  2  or  1  reed  switch  then  the  dis¬ 
engaging  coil  of  No.  6  dual  coil  relay  is  energized  which  then 
registers  an  interrupt  on  Channel  No.  6  of  the  interrupt  card  in 
the  multiprogrammer  unit.  The  computer  recognizes  that  this 
interrupt  indicates  a  high  diaphragm  and  in  turn  operates  two 
solenoid  valves  to  let  water  out  of  the  system  and  add  N2  to  the 
top  side  of  the  accumulator  which  will  result  in  lowering  the 
diaphragm.  Once  the  disengaging  coil  is  energized,  it  will 
remain  so  sending  an  interrupt  to  the  multi  programmer  unit.  This 
will  continue  until  the  diaphragm  moves  down  so  that  the  magnet 
again  closes  the  No.  3  reed  switch  and  as  a  result  energizes  the 
engaging  coil  and  discontinues  the  interrupts.  Likewise,  a  low 
diaphragm  will  register  an  interrupt  on  Channel  No.  5  of  the 
interrupt  card.  This  prompts  the  computer  to  cause  the  pump  to 
pump  more  water  into  the  system  and  let  some  N2  out  by  means  of 
another  solenoid  valve.  This  set-up  therefore  allows  for  pre¬ 
programmed  pressure  changes  during  testing  and  ensures  that  the 
optimum  diaphragm  positioning  is  obtained  at  any  pressure. 

Misalignment  Control--Another  feature  required  for  the  simu¬ 
lated  submarine  operation  test  program  is  a  variable  tilt  and 
offset  of  the  seal.  This  was  accomplished  by  the  set-up  shown  in 
Figure  2-16.  A  Slo-Syn  SS  250-1027  stepping  motor  and  gear 
reduction  arrangement  is  connected  to  the  shaft  of  the  worm  which 
drives  the  waviness  cylinder.  This  waviness  cylinder  has  a  tilt 
and  offset  machined  into  it  which  begins  just  right  of  the  end 
plate.  To  the  left  of  this  point  the  waviness  cylinder  is  con¬ 
centric  with  the  shaft.  Using  the  computer  to  pulse  the  stepping 
motor,  three  different  tilts  with  corresponding  offsets,  as 
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Interrupt 


3  Operating  Point  Magnetic  Reed  Switch 


4  Lower  Dead  Zone  Magnetic  Reed  Switch 

5  Lower  Bound  Magnetic  Reed  Switch 

6  Dual  Coil  Down  Relay 

7  Dual  Coil  Up  Relay 


Figure  2-15.  Accumulator  Diaphragm  Positioning  Circuit. 
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specified  by  the  test  program,  can  be  obtained.  Details  of  the 
geometry  related  to  this  procedure  are  described  in  Chapter  3. 

Seventy-Five  Percent  Balance  Ratio  Test 

Test  No.  Ill,  Appendix  A,  was  run  at  a  balance  ratio  of  0.75 
and  zero  initial  taper  so  as  to  obtain  additional  wear  rate 
data.  The  operation  of  this  test  was  somewhat  similar  to  prev¬ 
ious  tests  run  under  the  same  conditions.  One  difference  was 
that  leakage  did  not  occur  until  some  390  hours  into  the  test. 
Similar  tests  showed  leakage  initially.  Also,  whereas  the  aver¬ 
age  torque  level  was  comparable  to  other  tests,  there  were  large 
fluctuations  associated  with  it  in  the  initial  start-up.  These 
became  lower  as  the  test  proceeded.  The  test  ran  much  more 
smoothly  after  48  hours  of  operation  and  up  to  approximately  412 
hours  where  the  torque  levels  again  began  to  fluctuate.  The  test 
was  eventually  shut  down  due  to  high  torque  at  494  hours. 

This  type  of  operation  at  the  end  of  the  test  is  comparable 
with  Test  Nos.  75  and  76  [5],  and  final  radial  traces  of  the  car¬ 
bon  face  also  showed  a  taper  of  about  -1000  pm/m.  Test  No.  75 
had  a  final  taper  of  -710  pm/m  and  Test  No.  76  -982  pm/m.  This 
behavior  is  to  be  compared  to  that  for  Test  66  which  was  also  for 
a  zero  initial  taper.  Test  66  started  up  with  low  torque  and 
some  leakage  and  continued  about  the  same  for  over  100  hours. 
Clearly,  Test  111,  which  was  outwardly  identical,  started  out  in 
a  different  mode.  Test  66  started  out  taking  advantage  of  a 
thermal  taper  which  never  wore  off  during  the  test.  Test  111 
never  got  into  this  mode.  It  somehow  got  started  in  a  mode  where 
it  wore  rapidly  into  a  nonleaking  parallel  face  condition  with  a 
high  but  worn  away  thermal  taper.  The  last  hours  of  operation  of 
Test  111  show  that  this  mode  is  not  really  stable  but  can  lead  to 
high  friction,  high  wear  operation.  The  intriguing  auestion 
which  requires  still  more  testing  to  answer  is  why  this  test 
started  in  the  high  friction  mode  and  Test  66  started  in  a  low 
friction  mode. 


CHAPTER  3 

NINE-WAVE  SEAL  DESIGN 


Criteria 

In  this  chapter  the  design  procedure  for  the  nine-wave  seal 
is  presented  in  detail.  The  criteria  set  forth  for  the  design 
are  as  follows. 

1)  Very  low  wear--10  year  life, 

2)  Moderate  to  low  leakage--leakage  to  be  consistent  with 
time, 

3)  Low  friction  (to  ensure  low  thermal  distortions) 

4)  Operation  in  the  500  psi,  1800  rpm  at  4  inch  mean 
diameter  range, 

5)  Operation  in  seawater, 

6)  Variable  speed,  pressure,  and  distortion  as  functions  of 
time  to  conform  to  simulated  submarine  operation, 

7)  Seal  components  themselves  must  be  reliable  to  be  compa¬ 
tible  with  1)  above, 

8)  Use  of  carbon  insert  in  a  metal  ring  to  be  compatible 
with  large  scale  design. 

To  satisfy  the  above  conditions  using  a  wavy  seal,  certain 
additional  criteria  become  apparent.  Some  of  these  are  based  on 
conclusions  from  the  previous  report  [5]  and  a  paper  on  this  sub¬ 
ject  [14]. 

1)  The  seal  rings  must  be  designed  so  that  the  pressure 
caused  rotation  is  zero.  This  is  essential  in  a  seal  where  oper¬ 
ating  pressure  changes  frequently.  If  pressure  caused  rotation 
occurs,  then  the  seal  face  must  wear  to  a  new  profile  at  each 
operating  pressure. 

2)  Nine  waves  must  be  imposed  on  the  seal.  References  [5] 
and  [14]  show  that  the  natural  pressure  variation  on  the  seal 
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face  produces  a  tilt  angle  which  is  dependent  on  operating  condi¬ 
tions.  Such  changes  will  lead  to  accelerated  wear  and  must  be 
minimized.  By  increasing  the  number  of  waves  the  relative  tilt 
stiffness  of  the  seal  increases  such  that  operating  conditions  no 
longer  produce  a  change  in  tilt  which  is  significant  relative  to 
the  applied  tilt.  At  nine  waves  the  tilt  changes  only  a  few  per¬ 
cent  with  variable  operating  conditions. 

3)  The  centroid  of  the  cross  section  must  be  located  so 
that  the  sealing  radius  is  near  the  inside  radius  r,.  This 
ensures  that  both  minimum  leakage  and  maximum  load  support  will 
occur.  This  imposes  a  severe  geometrical  constraint  on  the  cross 
section. 

4)  So  that  the  seal  may  easily  comply  with  the  mating  ring 
at  the  lower  harmonics  and  so  that  the  needed  waviness  can  be 
reasonably  imposed,  the  cross  sectional  properties  (moments  of 
inertia)  must  be  maintained  at  low  values  (no  stiff  rings).  This 
requirement  restricts  the  size  of  the  cross  section. 

Assumed  Conf igurat ion 

Figure  3-1  shows  the  assumed  basic  configuration  for  the 
design.  Several  other  configurations  were  evaluated,  but  the  one 
shown  appeared  to  be  the  best.  There  are  54  moment  arms  attached 
to  the  seal  cross  section,  27  on  the  inside  and  27  on  the  out¬ 
side.  This  number  represents  3  sets  of  18  arms  (see  details  in 
previous  chapter).  Eighteen  arms  (nine  inside  and  nine  outside 
are  required  to  produce  one  set  of  nine  waves.  The  force  is 
applied  hydraulically  as  shown.  0-rings  create  a  pressurized 
pocket,  and  this  pressure  acts  on  the  pad.  The  carbon  is  bonded 
into  a  metal  ring.  The  centroid  is  shown  to  be  located  near  the 
inside  radius  of  the  ring.  Since  the  forces  on  the  pads  alter¬ 
nate  in  direction,  there  is  no  net  radial  force  due  to  hydraulic 
pressure  acting  on  the  ring.  The  seal  support  must  drop  down 
after  the  secondary  0-ring  so  that  some  significant  amount  of 


Figure  3-1.  Assumed  Nine-Wave  Seal  Configuration. 
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material  can  be  attached  to  the  cross  section  below  the  balance 
radius  in  order  to  lower  the  centroid. 

The  diameter  of  the  existing  test  vessel  contains  the  out¬ 
side  diameter  of  the  seal  assembly,  and  this  restricts  the  size 
of  the  moment  arms.  The  fact  that  the  hydraulic  pressure  is 
applied  on  the  pressurized  side  of  the  seal  creates  a  problem  in 
that  the  hydraulic  pressure  must  always  be  above  the  sealed  pres¬ 
sure.  If  it  is  not,  the  0-rings  can  collapse  inward  as  shown. 
Even  if  this  problem  is  overcome  by  inserts  inside  the  0-rings, 
it  can  be  shown  that  a  force  discontinuity  will  occur  as  pressure 
drops  below  the  sealed  pressure.  According  to  theory  previously 
developed  [5],  the  average  hydraulic  pressure  (which  is  influ¬ 
enced  by  the  sealed  pressure)  does  not  affect  the  amplitude  of 
the  wave,  so  having  the  hydraulic  pads  in  the  sealed  pressure 
environment  causes  no  theoretical  problem. 

It  was  immediately  recognized  that  because  of  all  of  the 
0-ring  seals  operating  at  high  pressure,  no  significant  corrosion 
of  the  surfaces  could  be  tolerated.  It  was  learned  that  Inconel 
625  was  one  of  the  few  alloys  that  would  not  pit  in  seawater,  so 
it  was  chosen  for  the  metal  parts  of  the  seal  assembly.  Inconel 
625  is  known  to  be  difficult  to  machine  compared  to  stainless 
steel  alloys  and  monel,  but  it  appeared  to  be  the  only  acceptable 
alloy  from  the  corrosion  standpoint. 

Seal  Ring  Design  Solution 

Figures  3-2  and  3-3  show  the  details  of  the  seal  ring  it¬ 
self.  Many  of  the  dimensions  shown  had  to  be  selected  based  on 
the  criteria  discussed.  The  details  of  this  procedure  will  now 
be  described. 

1)  Fix  certain  parameters  (see  Figures  3-2  and  3-3) 

K  =  1.005  (explained  later) 

R-  =  1.9  in.  (consistent  with  test  machine) 


Rb  =  1.9  in.  (balance  ratio  =  1.0) 


I 

t 


1 


\ 


I 


I 


Figure  3-2.  Assumed  Seal  Ring  and  Moment  Arm  Geometry. 
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xx3  =  0.1875  in.  (face  width  is  equal  to  existing  test  seal) 
yy3  =  0.125  in.  (nose  height) 
ti  =  0.010  in.  (bond  thickness) 
t2  3  0.003  in.  (clearance) 

x6  =  0.051  in.  (required  for  proper  squeeze)  (3-1) 

y7  =  0.080  in. 

pQ  =  500  psi  (seal  working  pressure  for  design  purposes) 
n  =9  (number  of  waves) 

p  *  1000  psi  (equivalent  gas  pressure  for  calculation.  Using  a 
9  moving  wave  produced  by  three  sinusoidal  with  time  waves 
requires  that  each  of  the  waves  have  a  pressure  amplitude 
of  667  psi 

E  =  31  x  106  psi  (Inconel  625) 
u  =  0.3 


2)  Choose  the  remaining  set  of  dimensions.  These  will  be 
varied  to  obtain  the  desired  design.  They  are: 

xl,  x2,  *66,  yl,  y2,  y3,  y55,  y6  -  defined  as  shown  in 

Figure  3-2 

d  -  pressure  pad  0-ring  diameter  (3-2) 

<j>  -  desired  seal  from  tilt  amplitude 

3)  Assuming  a  value  for  R{_ ,  calculate  section  properties 
X,  7,  rc,  and  Ix.  These  are  calculated  by  breaking  the  cross 
section  up  into  rectangular  subparts  as  shown  in  Figure  3-3  using 
the  parallel  axis  theorem. 

4)  Calculate  JQ  (torsional  stiffness).  Torsional  theory 
requires  that  the  following  partial  differential  equation  be 
solved  for  the  cross  section. 


i 


\|>  =  0  on  boundary  . 
Once  i>  is  found,  then 


(3-3) 

(3-4) 


(3-5) 


This  problem  was  solved  numerically  by  dividing  up  the  cross  sec¬ 
tion  into  a  grid  as  shown  in  Figure  3-4.  Note  the  grid  pattern 
corresponds  to  the  exact  boundaries  of  the  cross  section.  This 
causes  unequal  spacing  of  the  grid  points.  Equation  (3-5)  was 
solved  using  finite  difference  formulas  which  account  for  unequal 
spacing.  The  system  of  equations  was  solved  quickly  by  relaxa¬ 
tion. 

5)  Calculate  A  =  EJX/GJ9.  Using  this  and  the  previous¬ 
ly  derived  formulas  [5], 


(n‘  -  1) 


?  ’ 


(3-6) 


Pq  d  Rpi  e 
r 

c 


(3-7) 


calculate  the  required  m6o.  Given 
Rp i  =  +  x5  +  x6  +  x7/2 


(3-8) 


which  places  the  pressure  force  near  the  middle  of  the  back  edge 
of  the  seal  ring,  calculate  the  moment  arm  e  needed  to  produce 
the  required  tilt. 


< 


50 


6)  From  before  [5],  calculate 


v 


(1  + 

(?^i) 


(3-9) 


and 


rcontact  <j>  +  rc  '  (3-10) 

actual 

7)  To  insure  proper  sealing  but  also  maximizing  load  sup¬ 
port,  the  contact  radius  must  be  just  inside  the  inside  radius  of 
the  seal  or 

'contact  *  Kri  ■  (3-H> 

desired 

K  was  defined  previously  as  a  number  greater  than  one.  Now  com¬ 
pare: 


rcontact  t0  rcontact  *  (3-12) 

actual  desired 


If  they  are  significantly  different,  adjust  using  a  root 
finding  technique.  Go  back  to  Step  3  and  repeat  all  of  the  cal¬ 
culations  until  Equation  (3-12)  is  satisfied. 

8)  Calculate  Iy.  Calculate  the  moments  due  to  pressure 
as  shown  in  Figure  3-2.  Calculate  pressure  caused  rotation 


pressure 


m  r^ 
9  c 

~TT~ 


(3-13) 


where  m9  is  the  moment  about  the  centroid  in  N*m/m  of  centroi- 
dal  c ircumference. 
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9)  The  above  program  is  rerun  adjusting  the  input  param¬ 
eters  (3-2)  until  pressure  caused  rotation  (3-13)  is  made  suffi¬ 
ciently  small. 

The  computer  program  written  to  perform  the  above  calcula¬ 
tions  is  included  in  Appendix  B. 

Final  Design 

After  many  trials  and  alterations  in  configuration,  the 
final  design  selected  is  described  by  the  output  from  the  program 
shown  as  Table  3-1.  Figure  3-5  shows  the  seal  cross  section 
which  corresponds  to  these  dimensions.  Figures  in  Chapter  2  show 
other  details  of  the  seal. 

Predicted  Performance 

Using  programs  developed  previously  [5],  performance  of  the 
new  seal  design  was  predicted.  Figures  3-6  and  3-7  were  obtained 
using  a  program  which  neglects  elastic  deflection  of  the  seal 
rings  brought  about  by  face  pressure  itself.  As  discussed  pre¬ 
viously  at  n  =  9  this  deflection  is  negligible.  Figure  3-6  shows 
a  significant  friction  torque  at  the  lower  speeds  with  a  signifi¬ 
cant  reduction  at  increased  speed.  Even  at  zero  speed,  predicted 
torque  is  about  one  half  of  what  it  would  be  without  any  wavi- 
ness.  Some  wear  is  expected  at  the  lower  speeds  because  the 
fraction  of  the  load  supported  by  fluid  pressure  drops  to  75  per¬ 
cent  at  200  rpm.  However,  because  the  speed  is  low  where  the 
fluid  pressure  load  support  is  low,  the  amount  of  wear  will  be 
minimized.  At  the  higher  speeds,  the  predicted  fluid  pressure 
load  support  becomes  100  percent  due  to  hydrodynamic  effects  as 
is  suggested  by  the  decreasing  torque  with  increasing  speed 
curve. 

Leakage  is  shown  in  Figure  3-7.  Leakage  is  predicted  to  be 
less  than  1.0  cc/min  under  all  test  conditions.  A  large  value  of 
increases  leakage  considerably  but  reduces  torque  very  little. 


Lr 


Table  3-1 
Seal  Design 


♦ 

= 

320  x 

10 

d 

= 

0.254 

in. 

*1 

= 

0.2100 

in. 

yi 

s 

0.3800 

in. 

*2 

= 

0.2100 

in. 

y2 

= 

0.2200 

in. 

*3 

= 

0.2075 

in. 

y3 

= 

0.1000 

in. 

*4 

= 

0.5010 

in. 

y4 

= 

0.1000 

in. 

*5 

= 

0.5140 

in. 

ys 

= 

0.1700 

in. 

X  6 

= 

0.0510 

in. 

ys 

= 

0.3500 

in. 

X? 

= 

0.1435 

in. 

y7 

= 

0.0800 

in. 

Rpi 

= 

2.4840 

in. 

= 

1.3890 

in. 

X 

= 

0.3518 

in. 

y 

s 

0.3943 

in. 

Ix  =  0.00542  in, 


I  =  0.00916  in.4 
JQ  =  0.00176  in.4 
A  =  7.86 

v  =  8.7  x  10"6  in, 
e  =  0.953  in. 
rr  =  1.9557  in. 


^pressure 


=  6.8  x  10"8  in. /in, 
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Figure  3-7. 


^dieted  Performance 


leakage 
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Thus,  the  41  =  320  was  selected  as  a  compromise  between  reduced 
friction  and  increased  leakage. 

Figure  3-8  shows  a  comparison  between  the  worn  in  shape  of 
the  seal  at  two  extreme  operating  conditions.  Note  the  worn  in 
shape  is  almost  the  same.  These  curves  were  calculated  using 
another  computer  program  which  accounts  for  the  deflections  of 
the  seal  rings.  Thus,  as  expected  at  nine  waves,  there  is  little 
change  in  the  worn  shape  of  the  seal  due  to  a  change  in  operating 
conditions.  Thus,  no  acceleration  in  wear  would  be  expected  in 
the  design  due  to  variable  operating  conditions.  Furthermore, 
the  additional  small  deflection  component  makes  virtually  no  dif¬ 
ference  on  the  predicted  performance. 

Finite  Element  Calculation 

To  provide  an  additional  check  on  the  calculation  for  zero 
pressure  rotation  contained  in  the  design  program,  a  completely 
independent  calculation  for  pressure  caused  rotation  was  made. 
Using  SAP  IV  two-dimensional  ax i symmetrical  elements,  the  cross 
section  shown  in  Figure  3-5  was  modeled  using  161  nodes  and  126 
rectangular  elements.  The  results  showed  that  a  rotation  of  19 
win. /in.  convergent  resulted.  This  is  quite  small.  This  result 
shows  that  the  beam  theory  used  in  the  design  program  is  adequate 
for  calculating  zero  pressure  moment  designs  and  that  the  design 
itself  is  correct  from  the  pressure  rotation  standpoint. 

Strength  and  Deflection  Calculations 

Room  temperature  properties  of  annealed  Inconel  625  are  as 
follows: 

Su  =  120  to  150,000  psi 

Sy  =  60  to  95,000  psi  (3-14) 

Se  =  90,000  psi  at  108  cycles 
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Figure  3-8. 


Several  strength  calculations  were  made.  Points  A  and  B  shown  in 
Figure  3-5  are  highly  stressed  points,  a  =  38,000  was  allowed  at 
these  points.  This  stress  determined  the  pad  thickness  and  the 
moment  arm  connection  system,  both  of  which  are  to  be  as  small  as 
possible  from  a  space  point  of  view.  It  was  assumed  that  the 
moment  arm  itself  does  not  deflect  because  of  the  deep  cross  sec¬ 
tion  of  the  arm. 

Other  Calculations 

Many  other  design  types  of  calculations  were  made.  These 
include: 

1)  Ninth  harmonic  radial  deflection  (effects  0-ring  seals) 

2)  Assembly  bolts  (Monel  and  Inconel) 

3)  0-ring  assembly 

4)  Spring  pressure--spring  design 

5)  Pressure  drop  in  oil  passages 

6)  Mating  ring  rotation 

7)  0-ring  friction 

Tilt  and  Offset 

One  of  the  test  requirements  for  the  simulated  submarine 
operation  test  is  to  introduce  known  and  variable  tilt  and  offset 
misalignments  in  the  seal.  The  purpose  is  to  simulate  the  mis¬ 
alignment  which  results  from  hull  deflection  with  variable  pres¬ 
sure.  In  the  test  apparatus  the  tilt  and  misalignment  are 
obtained  by  the  rotation  of  a  misaligned  mounting  cylinder. 

Figure  3-9  shows  the  geometry  of  the  required  tilt  and  off¬ 
set.  Tilt  and  offset  are  in  the  same  plane.  The  tilt  and  offset 
are  specif  iced  at  each  simulated  test  depth.  Figure  3-10  shows 
schematically  how  the  tilt  and  offset  are  achieved  in  the  test 
apparatus.  The  shaft  is  tilted  and  offset  relative  to  the  end- 
plate  as  shown.  The  mating  ring  axis  coincides  with  the  shaft 
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Secondary  Ring 


Figure  3-9.  Required  Offset  and  Tilt. 


Figure  3-10.  Method  of  Controlling  Offset  and  Tilt. 


axis.  The  seal  axis  is  fixed  to  the  waviness  cyclinder  axis  but 
it  is  in  turn  tilted  and  offset  relative  to  the  endplate.  The 
relative  offset  between  the  shaft  and  the  seal  are  of  interest. 
Thus  the  relative  offset  is  equal  to  the  shaft  offset  vectorially 
minus  the  seal  offset.  The  tilt  is  equal  to  the  shift  tilt  vec- 
torially  minus  the  seal  tilt.  The  angle  cq  represents  the  posi¬ 
tions  of  the  waviness  cylinder.  If  oq  changes  and  all  tilts  and 
offsets  are  held  fixed,  their  vector  sums  change  with  this  rota¬ 
tion.  Thus  by  varying  the  position  of  the  waviness  cylinder,  it 
is  possible  to  cause  a  variable  offset  and  tilt.  The  problem  of 
course,  is  to  find  all  of  the  various  tilt  and  offset  angles 
which  give  the  desired  results. 

Before  the  method  of  solution  is  presented,  first  the  mathe¬ 
matics  of  the  problem  must  be  put  forth.  The  parameters  over 
which  one  has  control  either  by  manufacture  or  adjustment  are: 

3i  -  seal  tilt  relative  to  endplate 
62  -  shaft  tilt  relative  to  endplate 
y  -  angle  of  62  relative  to  ej 
ex  -  shaft  offset  relative  to  endplate 
e2  -  seal  offset  relative  to  endplate 

cq  -  angle  of  e2  relative  to  ex 

Note  61  is  taken  as  perpendicul ar  to  e2.  This  is  done  to  simply 
eliminate  a  variable.  The  needed  offset  and  tilt  are  shown  in 

Figure  3-10  as  vector  sums  of  the  various  components.  Note  that 

5  and  £'  should  be  approximately  equal  so  that  the  tilt  and  off¬ 
set  are  in  the  same  plane. 

Now  given  the  above  values 


2e,eo  cos  a 


(3-19) 


S'  =  cos 


-1 


2  2 
o  +  &2 


258- 


-  y 


Thus,  the  required  parameter  values  may  be  calculated  quite  read¬ 
ily- 

After  some  considerable  trial  and  error,  the  following  cal- 
culational  method  was  devised.  (Note  from  Table  3-2  that  d  and  5 
are  given  as  three  sets  of  values  which  will  be  given  subscripts 
a,  b,  and  c.) 


1)  Choose  ei,  e2,  8i  and  02. 

2)  Calculate  a 2  from  (3-16)  and  a!  from  the  law  of  sines 
for  all  d. 

3)  Find  d5/dcq  =  ( 5C  -  5a)/(aic  “  ala)* 

4)  Find  5b/(d5/dai)  from  above. 

5)  Now,  from  (3-18) 


6  B1  +  82  ■  2s1B2  cos^al  "  y) 
13T"  8^2  sin(a^  -  y) 

daj 


(3-20) 


Solve  this  for  (aib  -  y) . 

6)  Find  y  since  aib  is  known. 

7)  Calculate  actual  Sa5b5c  and  c  and  s'.  Compare  to 
desired  results. 


( 

*  ^ 

►  *  ■ 
t  „ 

l' 

S 


Table  3-2 

Offset  and  Tilt  Results 


Position 

Desired 

Offset 

Desired 

Tilt 

Computed 

Offset 

d 

Computed 

Tilt 

<5 

5'  -  6 

a 

1 

a 

0.058  in. 

0.0002 

0.058  in. 

0.00028 

21° 

56.2° 

b 

0.063  in. 

0.0005 

0.063  in. 

0.00050 

0° 

61.6 

c 

0.067  in. 

0.0007 

0.067  in. 

0.00070 

8° 

66.0 

% 

Computed  results  for:  e!  -  0.063  e2  =  0.060 

6 i  =  0.0029  $2  =  0.0027 

■  y  »  52.1’ 
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specified  by  the  test  program,  can  be  obtained.  Details  of  the 
geometry  related  to  this  procedure  are  described  in  Chapter  3. 
Table  3-2  shows  the  final  results  for  the  selected  parameter 
values  shown.  The  offset  values  can  be  matched  exactly.  It  was 
found  that  there  was  no  way  to  exactly  match  the  tilt  values;  so 
there  is  some  error  in  Sa.  It  is  also  not  possible  to  make  6 
exactly  perpendicular  to  d  at  all  positions.  £'  -  £  shows  the 
error  in  the  first  and  last  positions. 

Given  the  reasons  for  introducing  the  tilt,  it  was  judged 
that  the  above  is  sufficiently  close  to  the  desired  to  assess  the 
effects  of  variable  offset  and  tilt  in  operation.  At  the  same 
time,  the  tilt  and  offset  are  readily  changed  by  rotating  the 
wav i ness  cylinder.  This  can  be  performed  by  computer  control  and 
thus  saves  a  great  amount  of  repetitive  manual  adjusting. 

Summary 

The  preceding  sections  provide  brief  descriptions  for  the 
computational  portion  of  the  nine-wave  seal  design.  Details  of 
the  hardware  itself  are  given  in  Chapter  2. 
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CHAPTER  4 

MOVING  WAVE  GAS  SEAL 


Background 

There  are  several  types  of  gas  seals  being  used  in  industry 
today,  but  due  to  the  demands  of  advanced  technology,  a  more 
reliable  gas  seal  is  still  needed.  Experimental  work  on  the  mov¬ 
ing  wave  liquid  seal  has  demonstrated  some  promising  results  in 
minimizing  leakage,  friction,  and  wear  using  this  concept.  The 
potential  advantages  of  a  moving  wave  gas  seal  are  that  1)  the 
seal  will  not  foul  by  dirt  or  oil  because  the  faces  are  contin¬ 
uously  wiped  by  the  moving  wave  action,  2)  the  seal  will  not  be 
detrimentally  affected  by  liquid  contamination  since  it  can  be 
designed  to  operate  in  liquid  and  gas  environments,  3)  uneven 
wear  caused  by  hot  spots  or  contact  can  be  smoothed  out  by  the 
continual  wearing  in  process  of  the  seal,  and  4)  the  seal  is 
stiff  and  good  tracking  at  high  speeds  would  be  expected.  In 
this  chapter,  this  concept  has  been  applied  to  a  gas  seal.  A 
moving  wave  gas  seal  has  been  modeled  and  the  results  are  pre¬ 
sented  herein. 

State-of-Art  of  the  Gas  Seals--It  was  not  until  the  late  50s 
and  early  60s  that  gas  bearings  and  gas  interface  gas  seals 
respectively  were  developed.  In  order  to  satisfy  different  kinds 
of  applications  and  requirements,  various  types  of  gas  seals  have 
been  tried.  There  are  about  four  main  types  of  gas  seals. 

A  buffered  gas  seal  is  a  double  seal  configuration  (Figure 
4-1)  using  oil  as  a  buffer.  The  purpose  of  the  buffer  oil  is  to 
guard  against  the  escape  of  the  sealed  fluid  to  the  external 
atmosphere  and  also  to  provide  lubrication  between  the  seal 
faces.  The  buffered  seal  can  seal  corrosive  gases  up  to  750  psi 
(5.17  MPa)  with  rubbing  speeds  up  to  350  ft/s  (106.7  m/s).  The 
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1 .  Primary  Ring 

2.  Secondary  Ring 

3.  Spring 

PM 

>  Sealed  Gas'  Pressures 

P2l 

P3  Buffered  Oil  Pressure 


Figure  4-1.  Buffered  Gas  Seal. 


oil  pressure  is  about  30  to  50  psi  (0.21  to  0.34  MPa)  higher  than 
the  system  gas  pressure.  The  buffer  oil  also  acts  as  a  coolant 
to  the  system.  The  major  disadvantages  of  this  gas  seal  are  that 
a  pressure  pump,  a  filter,  and  a  heat  exchanger  are  required  to 
deliver  the  buffer  oil  from  an  external  source  [25]  and  that  the 
oil  may  contaminate  the  process  fluid. 

A  shrouded  Rayleigh  step  gas  seal  is  shown  in  Figure  4-2. 
Basically,  it  contains  a  pumping  land  or  dam  and  a  pocket  which 
is  surrounded  by  a  sealing  area  of  small  gap.  There  are  other 
configurations  as  shown  in  Figure  4-2.  When  the  secondary  seal 
ring  rotates,  the  Rayleigh  pockets  will  scoop  up  the  sealed  gas 
and  raise  the  gas  pressure;  at  the  same  time  the  lands  around  the 
pocket  will  restrict  the  gas  outflow.  So  the  performance  of  this 
kind  of  seal  depends  on  the  land  area  to  pocket  area  trade-off. 
This  implies  that  if  the  land  is  wide,  less  area  is  available  for 
pressure  generation,  but  sealing  is  good.  According  to  Cheng, 
Castelli  and  Chow  [25],  the  leakage  of  this  seal  is  largely 
determined  by  hydrodynamic  action;  the  stronger  the  hydrodynamic 
action  is,  the  more  the  leakage  will  be.  Quantitatively  speak¬ 
ing,  the  leakage  rate  is  slightly  higher  than  that  of  the 
spiral-groove  orifice  gas  seal  discussed  later.  The  change  of  a 
pressure  ratio  (P0/P-j)  has  little  effect  upon  the  stiffness 
as  well  as  the  lead  support. 

According  to  Zuk,  Ludwig  and  Johnson  [27],  NASA  has  spon¬ 
sored  tests  of  the  Rayleigh  seal  for  future  supersonic  tranport 
engines.  The  performance  demonstrated  the  feasibility  of  opera¬ 
tion  at  gas  temperatures  up  to  1200°F  (650°C),  pressure  differen¬ 
tials  across  the  seal  up  to  250  psi  (1.7  MPa)  and  relative  sur¬ 
face  speeds  up  to  450  ft/sec  (140  m/sec).  After  338.5  hours  of 
testing,  the  average  wear  on  the  carbon  seal  was  less  than  5 
yin.  (0.13  urn).  In  the  total  500  hours  test,  the  seal  had 
encountered  over  50  startups  and  shutdowns  and  the  leakage  rates 
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varied  from  11  to  32  SCFM  (0.31  to  0.91  SCMM)  with  an  average 
leakage  about  25  SCFM  (0.71  SCMM). 

Among  the  various  kind  of  gas  seals,  the  so-called  spiral- 
groove  gas  seal  is  well  known.  As  shown  in  Figure  4-3,  the  sur¬ 
face  of  the  secondary  seal  has  been  machined  or  etched  in  the 
form  of  spiral -grooves  which  are  equally  spaced.  The  spiral  does 
not  extend  across  the  full  width  of  the  seal  but  terminates  in  a 
circumferential  sealing  dam.  The  purpose  of  the  spiral  is  to 
pump  the  sealed  gas  inward  while  the  dam  restricts  radial  leak¬ 
age.  Thus,  the  gas  is  pressurized  over  the  whole  area  of  the 
seal  and  the  seal  load  is  supported  by  this  pressure.  According 
to  Sedy  [28],  leakage  rates  as  low  as  0.1  SCFM  at  850  psig  (5.86 
MPa)  gas  pressure  on  4.5  inches  (114.3  mm)  diameter  seal  with 
10,380  rpm  speed  have  been  recorded.  Typically  the  leakage  is 
about  1.9  SCFM  (0.05  SCMM)  in  the  same  test  and  seal.  The  wear 
which  occurred  on  the  seal  face  during  this  test  was  about  4 
uin.  (0.1  un)  after  1000  hours. 

At  the  present  state-of-art,  existing  gas  seals  are  satis¬ 
factory  for  many  applications;  leakage  is  relatively  low  and 
shaft  speeds  up  to  20,000  rpm  can  be  used  for  smaller  seals. 
However,  improvements  are  still  needed.  For  instance,  the  gas 
film  riding  seals  above  are  sensitive  to  the  detrimental  effects 
of  dirt  or  liquid  contamination.  Dirt  or  liquid  will  alter  the 
pressure  distribution  and  abrasive  wear  may  occur.  Seal  leakage 
rate  and  wear  may  dramatically  change.  High  spots  occur  on  most 
mechanical  seal  faces.  As  a  result,  hot  spots  may  occur  due  to 
highly  localized  friction,  and  wear  as  well  as  driving  torque  may 
ultimately  become  higher.  This  behavior  can  lead  to  complete 
seal  failure.  Most  gas  film  riding  gas  seals  have  relatively 
high  leakage  rates.  With  regard  to  these  problems,  further 
improvements  in  the  state-of-art  would  be  most  beneficial.  In 
pursuit  of  such  improvements,  the  wavy  gas  seal  deserves  study. 
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Groove  Gas  Seal 


In  terms  of  the  required  analysis  the  gas  seal  is  somewhat 
different  than  the  liquid  seal  studied  previously  herein.  Thus 
it  is  useful  to  review  the  literature  for  studies  of  other  types 
of  gas  seals  to  better  understand  what  is  important  and  the 
computational  methods  available. 

Literature  Review — In  the  late  1960s,  many  successful  and 
efficient  techniques  of  seal  analysis  were  developed.  First, 
James  and  Potter  [29]  solved  for  the  pressure  distribution  for 
both  spiral-groove  thrust  bearings  and  spiral -groove  compres¬ 
sors.  They  derived  a  "jump"  boundary  equation  between  the  land 
and  the  groove,  and  then  made  use  of  variable  transformations  to 
straighten  out  the  groove  for  easier  computation.  To  facilitate 
further  the  analysis,  they  changed  the  dependent  variable  from  P 
to  Q  where  Q  =  1/2  P2.  Finally,  they  used  an  iterative  technique 
to  find  the  pressure  distribution.  Given  the  pressure  distribu¬ 
tion,  the  gas  flow,  load  support,  friction,  static  stiffness, 
power  consumption  and  efficiency  can  be  calculated. 

Cheng,  Castelli  and  Chow  [26]  made  detailed  analyses  of  two 
high  speed  noncontacting  hybird  seals,  the  spiral-groove  and 
Rayleigh  step  seals.  The  first  part  of  their  paper  shows  the 
performance  of  the  spiral -groove  seal  with  orifices  and  then 
without  orifices.  The  differences  in  performance  between  the 
cases  where  the  grooves  are  located  at  the  high  pressure  and  the 
low  pressure  or  both  sides  were  also  studied. 

In  the  Rayleigh  step  real  analysis  [26],  several  assumptions 
were  made:  laminar  and  isothermal  flow,  Newtonian  fluid,  and 
negligible  inertia.  With  the  change  of  variable  P2  =  Q,  the 
Reynolds  equation  was  expressed  in  finite  difference  form.  The 
result  is  a  set  of  equations  with  the  pressure  (Q)  at  five  points 
operated  un  by  a  linear  algebraic  operator.  These  equations  were 
solved  by  the  so-called  matrix  inversion  technique  developed  from 
Reference  [30].  This  technique  enables  one  to  formulate  the 


m  x  n  unknowns  in  m  x  n  linear  equations  with  known  coeffi¬ 
cients.  However,  because  pressure  P  acts  as  a  variable  coeffi¬ 
cient  (which  is  assumed  to  be  a  constant  for  each  iteration),  an 
iterative  procedure  was  adopted  to  improve  the  guess  values  of  P. 

The  second  part  of  the  Cheng,  Castelli  and  Chow  paper  [26] 
gives  typical  performance  for  a  non-parallel  film  profile  spiral- 
groove-orifice  seal  as  well  as  for  the  Rayleigh  step  seal.  Fin¬ 
ally,  a  comparison  was  made  between  the  two  seal  geometries  on 
their  tolerance  to  tilt  or  coning  under  a  constant  seal  load. 

The  authors  concluded  that  the  difference  in  stiffness  and  leak¬ 
age  between  the  two  seals  is  small,  but  the  spiral-groove-orifice 
seal  shows  small  variations  in  leakage  rate  as  well  as  axial 
stiffness  with  respect  to  tilt  or  coning. 

Gardener  [31]  has  summarized  and  presented  a  numerical 
method  for  solving  the  design  problem  of  a  spiral  groove  noncon¬ 
tacting  seal.  Gardener  considered  pressure  to  be  a  function  of 
radius  only,  so  that  the  analysis  was  simplified.  Some  time 
later.  Gardener  [32]  presented  a  paper  covering  additional  devel¬ 
opments  on  noncontacting  spiral-groove  seals.  He  pointed  out 
several  modifications  to  improve  performance. 

Zuk,  Ludwig  and  Johnson  [27]  presented  an  approximate  quasi- 
one-dimensional  model  for  compressible  fluid  flow  across  face 
seals.  Both  the  flat  face  and  Rayleigh  step  gas  seals  were  con¬ 
sidered  and  tested.  The  model  includes  fluid  inertia,  viscous 
friction  and  entrance  losses.  It  is  valid  for  both  laminar  and 
turbulent  flows.  Both  subsonic  and  choked  flow  conditions  can  be 
predicted  by  using  an  approximate  integrated  model  (an  iterative 
technique).  Results  agree  with  classical  subsonic  compressible 
viscous  flow  theory  for  Mach  number  less  than  l//y»  where  y  is 
the  ratio  of  specific  heats.  If  Mach  number  is  equal  to  1,  sonic 
velocity  and  thus  choked  flow  will  occur.  Results  also  show  that 
near  critical  flow  conditions  (Mach  number  =  1),  the  film  stiff¬ 
ness  is  negative,  that  is,  as  film  thickness  and  Mach  number 


increase,  the  applied  load  actually  increases.  For  choked  flow 
the  analysis  predicts  a  high  positive  film  thickness.  The 
authors  also  found  that  as  the  sealing  dam  width  increases,  the 
critical  pressure  ratio  (ambient  pressure  over  sealed  pressure) 
for  choking  decreases. 

Later,  Zuk  [33]  presented  an  application  of  the  finite  ele¬ 
ment  method  to  compressible  problems  using  a  Galerkin  solution 
technique.  The  most  significant  aspect  of  this  paper  is  the 
inclusion  of  a  fluid  inertia  term  as  well  as  choked  flow  condi¬ 
tions  which  make  the  compressible  flow  equations  nonlinear.  He 
also  compares  the  advantages  and  disadvantages  between  the  finite 
element  method  and  the  finite  difference  method. 

In  the  late  1970s,  Sedy  [28]  found  that  the  performance  of  a 
spiral -groove  gas  seal  can  be  improved  through  enhancement  of 
hydrodynamic  effects.  Similar  results  are  also  applied  to 
Rayleigh  step  gas  seals.  The  methods  of  improving  the  perform¬ 
ance  include  1)  using  a  high  balance  ratio,  2)  widening  the  seal¬ 
ing  dam,  3)  increasing  the  overall  face  width,  and  4)  modifying 
the  carbon  sealing  ring  cross  section. 

Later,  Sedy  [34]  presented  a  new  self-aligning  mechanism  for 
improving  the  stability  of  spiral-groove  gas  seal.  That  is,  when 
originally  parallel  faces  of  the  mating  rings  deflect,  the  change 
in  pressure  profile  tends  to  return  the  deflected  ring  toward  its 
parallel  position.  It  does  not  matter  whether  the  deflection  is 
convergent  or  divergent.  If  the  gap  is  stable,  the  leakage  and 
radial  distortion  are  also  stable.  However,  to  benefit  from  this 
effect  the  pressure  between  the  spiral -groove  and  the  dam  must  be 
higher  than  the  sealed  pressure  and  the  centroid  location  must  be 
controlled.  This  mechanism  not  only  minimizes  angular  deflection 
of  the  primary  low-modulus  seal  ring,  but  also  minimizes  fluctua¬ 
tions  in  the  gap  width. 

Concerning  methods  of  solution  for  the  pressure,  from  a 
mathematical  point  of  view,  gas  seal  analysis  is  more  complicated 


than  liquid  seal  analysis.  Although  gas  seal  analysis  is  not  a 
new  branch  of  research,  relevant  information  and  techniques  to 
solve  for  the  pressure  distribution  are  still  being  developed. 

In  1968,  Caste 1 1 i  and  Pirvics  [35]  reviewed  different  kinds  of 
numerical  methods  for  the  solution  of  the  gas  bearing  problems. 
The  authors  described  some  of  the  techniques  using  explicit, 
implicit  and  semi-implicit  methods  for  the  time-dependent 
Reynolds  equation.  The  Alternating  Direction  Implicit  method  was 
also  mentioned.  Due  to  the  nonlinear  nature  of  the  compressible 
Reynolds  equation,  solution  techniques  such  as  Gaussian  Elimina¬ 
tion  become  just  a  part  of  a  larger  iteration  solution. 

According  to  Fairweather  and  Mitchell  [36],  Alternating 
Direction  Implicit  (A.D.I.)  methods  have  proven  to  be  valuable 
finite  difference  methods  for  various  types  of  partial  differen¬ 
tial  equations,  such  as  the  parabolic  and  elliptic  equations  in 
two  or  three  space  variables.  They  discussed  the  computation 
procedure  of  the  A.D.I.  method  from  a  pure  numerical  point  of 
view.  A  paper  by  Evans  and  Gane  [37]  is  a  good  example  of  the 
application  of  A.D.I.  They  were  successful  in  applying  this 
method  to  transient  heat  conduction  problems.  Also,  Roache 
[38]  has  a  concise  description  on  the  A.D.I.  method.  He  general¬ 
ized  the  experience  of  some  researchers  and  concluded  that 
A.D.I.  methods  do  allow  large  time  steps  and  fast  overall  compu¬ 
tation. 

One  of  the  more  useful  references  for  bearing  analysis  was 
written  by  Shapiro  [39].  He  gives  very  detailed  discussions  of 
numerical  techniques,  presents  several  examples,  and  makes  sug¬ 
gestions  for  future  development.  One  of  his  examples  is  the 
A.D.I.  scheme  applied  to  a  gas-lubricated  bearing.  He  used  a 
so-called  "cell"  approach  which  operates  on  the  continuity  equa¬ 
tion  rather  than  on  Reynolds  equation.  This  method  is  said  to  be 
more  convenient  in  solving  complex  goemetries  than  the  direct 
Reynolds  equation  analysis. 
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Finally,  according  to  Gargiulo  and  Gilmour  [40],  the  A.D.I. 
method  has  been  shown  to  be  efficient,  flexible  and  free  from 
cummulative  round  off  errors  for  the  design  of  externally  pres¬ 
surized  porous  gas  bearings.  They  used  a  systematic  method  to 
replace  the  Laplace  difference  operator  in  the  pressure  field  by 
a  finite  difference  approximation  according  to  the  A.D.I. 
scheme.  They  also  gave  a  set  of  linear  algebraic  equations  first 
in  radial  direction  and  then  in  the  circumferential  direction. 

The  two  tridiagonal  systems  were  efficiently  solved  by  an  elimin¬ 
ation  algorithm. 

Theory 

In  this  section,  a  model  of  a  moving  wave  gas  seal  is  pre¬ 
sented.  The  moving  wave  gas  seal  can  be  described  as  a  gas  seal 
that  has  both  moving  waviness  and  tilt  imposed  upon  the  seal  face 
of  the  primary  seal  ring  (see  Figure  1-3).  This  section  includes 
the  derivation  of  the  difference  equations  from  the  Navier-Stokes 
equations  and  techniques  for  obtaining  load  support,  leakage,  and 
friction.  Mesh  effect,  time  step  effect  and  convergence  criteria 
are  also  investigated.  Finally,  a  flow  chart  of  the  computer 
program  and  its  subroutines  and  example  solutions  are  presented. 

Compressible  Flow  Equations--Compressible  flow  equations  for 
a  thin  film  are  obtained  from  a  simplification  of  the  Navier- 
Stokes  equations.  The  simplifying  assumptions  used  are: 

1)  Gas  viscosity  is  constant 

2)  Perfect  gas 

3)  Body  forces  are  negligible 

4)  Flow  is  laminar,  there  are  no  vortices  or  turbulence 
anywhere  in  the  lubricant  film 

5)  Inertial  forces  are  negligible  compared  with  the 
viscous  forces 

6)  No  slippagae  occurs  at  the  sealing  faces 
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7)  The  curvature  of  the  fluid  film  surface  is  negligible 

8)  Compared  with  the  two  velocity  gradients,  3Vr/3z  and 
3V9/3z,  all  other  velocity  gradients  are  negligible. 

Since  these  assumptions  are  widely  accepted  for  gas  seal  analy¬ 
sis,  further  justification  will  not  be  given  here.  However, 
details  are  available  in  Reference  [23]. 

One  further  assumption  made  is  that  the  seal  surface  rough¬ 
ness  as  shown  in  Figure  4-4  does  not  significantly  affect  the 
resulting  flow  equations  and  that  the  mean  film  thickness  may  be 
used  instead  of  the  actual  film  thickness  for  computational  pur¬ 
poses.  This  assumption  is  discussed  in  detail  in  Reference  [4], 

By  applying  the  above  assumptions,  the  Navier-Stokes  equa¬ 
tions  in  polar  coordinates  can  be  simplified  to  the  following. 
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Figure  4-5  shows  the  applicable  coordinate  system, 
the  above  differential  equations,  the  result  is 
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Applying  the  boundary  conditions  shown  in  Figure  4-5 
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Figure  4-5.  Boundary  Conditions  and  Geometry. 
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we  then  obtain  two  component  equations  of  the  fluid  flow. 
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The  mass  flow  rate  is. 
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The  continuity  equations  for  an  element  of  volume  is 


(4-12) 
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With  Reference  to  Figure  (4-6)  the  continuity  equation  can  be 
written  as, 

3 

TtJ0  0  “  dz  +  Vl/2  *  Vl/2  *  m j - 1/2  +  Vl/2  =  0  ' 


(4-14) 


where  the  1/2  points  refer  to  mass  flow  halfway  between  nodes. 
The  mass  change  inside  the  control  volume  can  be  written  in  dif¬ 
ferential  form  as 
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where  Aa  =  surface  area  of  top  or  bottom  of  the  control  volume. 
Using  Equations  (4  ‘  (4-12)  the  mass  flows  are 
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Substituting  Equations  (4-15)  through  (4-19)  into  Equation 
(4-14),  assuming  a  perfect  gas,  isothermal  conditions  and  3h/at 
0  (no  squeeze),  one  gets 
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♦  te>[(prh))M/2  -  (Prh),_1/2]  =  0  .  (4-20) 
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To  linearize  the  problem  a  change  of  variable  is  made, 

Q  =  P2  (4-21) 

so  that 
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Incorporating  this  change  in  Equation  (4-20)  and  using  central 
differences  to  approximate  the  derivatives, 
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=  0  .  (4-23) 

where  Q  is  the  value  of  Q  at  a  previous  time  step  and  the  remain¬ 
ing  P* 1 s  are  assumed  known  from  a  previous  iteration.  Since  Q 
and  the  remaining  P's  are  assumed  known,  Equations  (4-23)  become 
linear  equations  at  any  time  step  with  variable  Q.  Grouping  all 
the  like  terms  as  coefficients,  the  above  equation  becomes. 
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A  special  Gaussian  elimination  subroutine  is  used  to  solve 
the  above  tridiagonal  equations.  It  stores  the  diagonal  coeffi¬ 
cients  in  a  compact  column  matrix  and  takes  advantage  of  the 
sparseness  of  the  matrix  to  save  computer  time. 


Boundary  Conditions--T'ne  boundary  conditions  for  the  gas 
seal  are  similar  to  those  for  a  liquid  seal.  That  is, 


P(rt)  =  P i  , 

Hr0)  =  P0  ,  (4-28) 

where  P^  and  P0  are  the  inside  and  outside  pressures  of  the 
seal,  respectively.  Within  the  seal  region  pressures  are  set 
initially  at  time  zero  to  values  based  on  linear  interpolations 
between  the  inside  and  outside  pressure. 

Mesh  effect — One  measure  of  solution  accuracy  is  to  compare 
the  leakage  rates  at  the  inside  and  outside  radii  of  the  seal. 
Figure  4-7  shows  this  comparison  as  a  function  of  tilt  for  two 
mesh  sizes.  A  different  grid  size  results  in  different  leakage 
rates,  but  the  difference  is  decreased  by  a  finer  mesh.  The 
error  increases  with  increasing  tilt  because  a  greater  change  in 
gap  across  the  seal  is  produced. 

Presumably,  the  leakages  will  converge  to  each  other  at  all 
tilt  angles  if  one  can  refine  the  mesh  as  small  as  possible.  In 
practice  one  cannot  refine  the  mesh  too  far  because  of  cost. 
Hence,  grid  size  comparison  is  terminated  at  the  19  x  20  mesh. 

All  the  results  obtained  are  based  on  this  mesh  scheme. 

Convergence  and  Time  Step--Using  the  ADI  method,  one  solves 
for  the  pressure  field  as  a  function  of  time.  In  this  steady- 
state  problem,  one  is  actually  solving  for  the  transient  response 
from  an  arbitrary  starting  pressure  field  to  that  which  is  the 
steady-state  solution  keeping  the  film  thickness  constant  with 
time.  A  measure  of  satisfactory  convergence  is  the  total  load 
supported  by  the  pressure.  Presumably,  when  this  load  stops 
changing  with  time,  steady  state  has  been  reached. 

In  order  to  check  the  stability  and  convergence  of  the 
program,  a  series  of  data  was  computed  and  is  plotted  in 
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Figure  4-8.  In  this  figure,  the  maximum  difference  in  absolute 
pressure  (aP)  between  the  old  iteration  values  and  the  new  values 
is  used  as  a  convergence  criterion.  The  effect  of  pressure  error 
and  the  time  step  on  the  fluid  load  support  and  the  total  compu¬ 
tation  time  is  shown.  The  figures  show  clearly  that  the  pressure 
error  can  be  used  to  get  a  satisfactory  convergence  on  load  sup¬ 
port. 

It  is  shown  that  a  larger  pressure  error  criterion  requires 
a  larger  time  step  to  obtain  proper  convergence  on  load.  Less 
computational  time  is  taken  using  this  approach.  If  the  pressure 
error  remains  the  same  while  the  time  step  is  decreased  as  shown 
in  Figure  4-8,  a  false  or  premature  convergence  results.  On  the 
other  hand,  if  the  pressure  error  is  decreased,  very  consistent 
load  support  is  found  at  various  time  steps.  However,  if  the 
pressure  error  is  decreased  even  more  (e.g.,  aP  <  5.5  Pa),  the 
execution  for  the  larger  time  step  (aT  =  5  x  10-5  S)  will  be 
terminated  because  negative  pressures  occur.  Also,  the  program 
loops  without  stopping  for  the  smallest  time  step  (aT  =  5  x 
10" 7  S)  because  the  pressure  error  (AP)  is  too  small  to  meet. 

As  shown,  care  must  be  taken  in  order  to  pick  a  suitable 
pressure  error  and  time  step.  Otherwise,  too  much  computer  time 
or  false  convergence  will  result.  In  this  investigation,  the 
pressure  error  criterion  is  taken  to  be  the  sum  of  the  absolute 
changes  for  all  the  nodes.  This  value  is  taken  at  3.5  KPa  for  18 
x  18  nodes.  This  corresponds  approximately  to  11  Pa  in  Figure 
4-8.  The  time  step  is  chosen  to  be  5  x  10"6  seconds. 

Seal  System--the  moving  wave  gas  seal  can  be  defined  by  two 
sets  of  parameters,  application  parameters  and  design  param¬ 
eters.  The  former  category  includes  the  pressure  differential  to 
be  sealed  (P0  -  P -j ) »  the  angular  speed  (u>)  of  the  shaft, 
shaft  diameter  (D),  seal  centroid  radius  (rc)  and  flu^d  viscos¬ 
ity  (n).  All  of  these  parameters  are  fixed  because  of  the 
requirements  of  a  particular  seal  application. 
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The  design  parameters  consist  of  balance  ratio  (B),  material 
properties  such  as  compressive  strength  (Pm),  seal  surface 
roughness  (c) ,  seal  face  width  (r0  -  rj),  number  of  waves 
around  the  seal  face  (n),  tilt  angle  (<j>),  amplitude  of  the  waves, 
the  ring  stiffness  ratio  (A),  spring  pressure  (P$p)  and  modulus 
of  elasticity  of  the  primary  seal  ring  (E).  The  values  of  the 
above  application  parameters,  design  parameters  as  well  as  some 
constants  are  given  in  Table  4-1.  Most  of  the  system  parameters 
are  shown  in  the  previous  figures,  such  as  Figure  1-3,  Figure 
4-4,  and  Figure  4-5.  The  definitions  used  are  the  same  as  those 
in  previous  reports. 

Film  Thickness  Shape--the  nature  of  the  wavy  tilted  shape  of 
the  gas  film  in  the  seal  is  described  in  detail  in  Reference 
[5].  Film  thickness  is  essentially  a  function  of  wear,  deflec¬ 
tion,  and  initial  shape.  For  a  moving  wave  or  tilted  seal,  the 
film  thickness  function  has  been  shown  to  be  [5]: 

h  =  hQ  +  w(r)  +  [vQ  +  (r  -  rc)  *0]  cos(n0  +  flt)  (4-29) 

where 

w(r)  =  -min{[vQ  +  (r  -  rQ)  $Q]  cos(n0  +  flt)}  (4-30) 

v0  and  <j)0  can  be  somewhat  arbitrarily  chosen,  and  their 
choice  influences  the  shape  of  the  film.  In  this  gas  seal  study, 
v0  and  <|>0  are  chosen  just  as  in  the  liquid  seal  work  [5]  such 
that  there  is  a  circumferential  wave,  a  radial  taper,  and  a  dis¬ 
tinct  sealing  dam.  n  represents  the  wave  motion  required  to 
maintain  the  wave  in  spite  of  wear  and  is  so  small  that  it  can  be 
neglected  in  the  pressure  computations.  Figure  4-9  shows  a 
three-dimensional  plot  of  the  wavy  film  thickness  for  a 
particular  +>  and  v. 


Table  4-1 


Gas  Seal  Parameters  and  Constants 
Face  Seal  Application  Parameters 
Viscosity  of  nitrogen  n  -  1.74  x  10"5  Pa»S  (2 


Outside  pressure 

Inside  pressure 

Angular  speed  of  shaft 
Seal  centroid  radius 

Shaft  diameter 

Face  Seal  Design  Parameters 

Ring  Stiffness  Ratio 

Compressive  strength  or 
carbon  ring 

Balance  ratio 

Surface  roughness 

Face  width:  inside  radius 


n  -  1.74  x  10-5  Pa»S  (2.52 
x  10"9  lbf-s/in.2) 

pQ  -  3.55  MPa  (514.7  psia)  (varies) 

p^  -  0.1  MPa  (14.7  psia) 

u  -  1800  rpm  (188  rps)  (varies) 
rc  -  49.18  mm  (1.9361  in.) 

D  -  88.90  mm  (3.50  in. ) 


A  -  EJX/GJQ  =  7.44 


pm  -  262  MPa  (38,000  psia) 

B  -  0.75 

c  -  0.51  pm  (20.0  pin.)  (varies) 
r-  -  48.26  mm  (1.9000  in.) 


outside  radius  rQ  -  53.04  mm  (2.0875  in.) 


Number  of  waves 

n 

-  9 

Waviness  amplitude 

ho 

-  0.14  pm  (5.4  pin.) 

Spring  pressure 

Psp 

-  0.196  MPa  (28.5  psi) 

Modulus  of  elasticy 

E 

-  20.7  GPa  (3.0  x  106  psi) 

Shear  modulus 

G 

-  9.6  GPa  (1.25  x  106  psi) 

Tilt 

♦ 

-  500  pm/m  (varies) 

Other  Constants 

Gr avity 

9 

-  9.81  m/s2  (32.174  ft/s2) 

Density  of  nitrogen 

P 

-  1.18  Kg/m3  (0.0023  slug/ft3) 

Spring  constant 

k 

-  3.36  KN/m  (19.159  lb/in.) 

Temperature 

T 

-  23 . 9°C  (75.08F) 

Universal  constant 

U 

-  287  J/Kg°K  (53.3  ft-lb/lbm/R0 ) 

i 


Narrow  Band 
of  contact 


Load  Support--Load  support  consists  of  the  fluid  pressure 
component  and  the  mechnical  component.  The  fluid  pressure  com¬ 
ponent  is  due  to  the  hydrostatic  and  hydrodynamic  effects  while 
the  mechanical  component  is  obtained  by  flattening  the  peaks  of 
the  asperities  in  contact.  As  before,  it  is  assumed  that  asper¬ 
ity  pressure  is  equal  to  the  compressive  strength  of  the  softer 
seal  face  material,  and  that  this  pressure  remains  constant. 

The  total  load  support  in  the  gas  seal  is  calculated  by 
integrating  the  hydrodynamic,  hydrostatic  and  asperity  pressure 
distributions  over  the  entire  seal  face  area.  That  is 

o 

[P  +  Pm(l  -  b)]  r  dr  de  ,  (4-31) 

where  b  is  the  fraction  of  the  face  not  in  mechanical  contact 
(given  in  Reference  [6]). 

The  load  that  must  be  supported  for  the  seal  to  be  in  equi¬ 
librium  is  a  function  of  the  seal  pressures  P0  and  P-j ,  the 
balance  ratio  B,  and  the  spring  pressures  Psp. 

W*  =  if (r q  -  r?)  [P03  +  P.(l  -  B)  +  P$p]  (4-32) 

As  before  [1-5],  the  total  load  support  W  increases  with 
decreasing  h0,  thus,  this  load  support  may  be  adjusted  to  be 
equal  to  the  applied  load  W*  by  selecting  a  suitable  h0.  Such 
a  procedure  is  used  for  the  solutions  presented  herein. 

For  comparison  purposes,  the  ratio  of  fluid  pressure  load 
support  to  the  total  load  (Wf/W)  expressed  as  a  percent  is  an 
index  to  indicate  the  relative  wear  rate.  If  the  ratio  is  equal 
to  one  (Wf/W  *  1.0)  theoretical ly,  adhesive  wear  will  become 
zero. 
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leakage--Gas  seal  leakage  is  computed  by  the  integration  of 
the  equations 


Q 


out 


h3 

hO  3P 
12n  3r 


de 


(4-33) 


and 


h?  3P 
12n  3r 


r.  de  , 


(4-34) 


Q-jn  and  Qout  represent  the  inside  and  outside  leakage,  res¬ 
pectively.  Both  the  inside  and  outside  leakages  are  used  as  a 
check  to  see  whether  accuracy  is  acceptable.  A  numerical  average 
of  the  inside  and  outside  leakages  is  used  to  represent  the  nom¬ 
inal  seal  leakage. 


Friction--From  the  previous  work  [6],  considering  the  fluid 
film  regime,  tangential  fluid  friction  shear  stress  is  given  by 

Tf  =  [noir  E(-jlj-)  +tp|^e(H)]  X  b  (4-35) 

where  E(  )  is  the  expectancy  operator  and  b  is  the  fraction  of 
the  face  not  in  contact  and  has  been  given  previously  [6]. 

Roughness  is  taken  as  being  given  by  a  parabolic  distribu¬ 
tion  as  before, 

f(M  =  -2L-  (c2  -  h2)3  ,  -c  <  h  <  c  .  (4-36) 

^  32c7  5 

The  expected  value  functions  are  given  in  Reference  [6], 

The  average  friction  shear  stress  due  to  asperity  contact  is 
given  by, 
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Ta  *  Ps*1  '  b> 


(4-37) 


The  mechanical  torque  is  calculated  by  the  integration  of 
the  tangential  fluid  shear  stress  plus  the  average  frictional 
shear  stress  due  to  the  asperity  contact  over  the  whole  seal 
face.  Thus, 


n  ^ 

rn  r  o  ? 

Torq  =  n  /  /  (t,  +  tJ  r  dr  de 

Jo  Jr.  f  a 


and  the  mechanical  driving  power  required  is. 


Power  =  a)  x  Torq  . 


The  coefficient  of  friction  is  defined  as 


Torq/r, 


where  r f  is  the  mean  seal  radius. 


(4-38) 


(4-39) 


(4-40) 


Stiffness--The  stiffness  in  the  gas  seal  is  computed  numeri¬ 
cally  by  varying  h0  by, 


-  dW 


(4-41) 


Example 

The  solution  technique  for  the  pressure  has  already  been 
discussed  in  previous  sections.  Considering  the  total  solution, 
h0  must  be  adjusted  to  satisfy  equilibrium.  Therefore,  h0 
becomes  a  dependent  variable  in  the  problem. 
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Specifically,  the  method  of  solution  is: 

1)  Select  a  system  of  seal  parameters  and  determine  the 
applied  load  W*. 

2)  Calculate  the  wear  term  w(r)  using  Equation  (4-30). 

3)  Assume  a  value  for  minimum  film  thickness  h0. 

4)  By  using  the  A.D.I.  technique  as  described,  the  pressure 
distribution  is  solved  for  the  given  boundary  conditions 
as  described.  A  typical  resultant  pressure  distribution 
is  shown  in  Figure  4-10. 

5)  Calculate  the  total  load  support  W  using  Equation 
(4-31). 

6)  Compare  W  to  W*  and  adjust  h0  using  a  secant  root 
finding  method.  Return  to  Step  4)  and  repeat  until  W  = 
W*  within  a  desired  accuracy. 

7)  Calculate  coefficient  of  friction,  leakage,  percent  load 
support  and  torque. 

Appendix  C  contains  a  computer  program  listing  for  the  solu¬ 
tion  to  this  problem. 

Using  the  parameter  values  in  Table  4-1,  an  example  solution 
was  obtained.  The  film  thickness  is  similar  to  that  shown  in 
Figure  4-9.  The  values  of  pressure  are  shown  in  Figure  4-10. 
Pressure  gradient  in  the  outer  edge  of  the  seal  is  low  and  at  the 
inner  edge  nigh.  This  suggests  that  accuracy  of  the  leakage  at 
the  outer  edge  (Qout)  1S  good  while  the  accuracy  at  the  inner 
edge  (Q-jn)  is  low.  It  is  much  more  difficult  numerically  to 
predict  the  steep  gradient  accurately.  The  load  support  required 
for  equilibrium  is  982.3  lbf  (4.4  KN).  The  total  load  support 
determined  by  the  fluid  load  and  asperity  load  is  found  to  be 
982.26  lbf.  Thus,  equilibrium  (W  =  W*)  is  attained.  The  minimum 
film  thickness  h0  is  30.2  yin.  (0.76  pm).  Since  h0  >  c  (0.5 
un),  no  contact  occurs.  NITA,  NITB  and  NITN  are  the  numbers  of 
iterations  in  the  A.D.I.  method  required  for  the  two  initially 
guessed  minimi-m  film  thicknesses  and  the  equilibrium  minimum  film 
thickness  respectively.  The  stiffness  of  the  seal  is 


Figure  4-10. 


high  (-1.15  x  106  Ibf/in.)  where  the  minus  sign  represents 
increasing  load  with  decreasing  gap.  The  leakage  rate  is  0.12 
SCFM  (3.4  x  10'3  SCMM) .  Negative  Qin  and  Qout  represent 
leakage  flow  radially  inward  as  expected  for  this  example.  The 
coefficient  of  friction  is  only  0.000045  which  implies  that  the 
frictional  force  is  small  (only  0.043  Ibf).  One  hundred  percent 
fluid  load  support  occurs.  Total  driving  torque  is  only  0.088 
in. -lb.  Thus,  both  the  leakage  and  torque  are  low.  Using  such 
results,  the  studies  below  were  completed. 

Parameter  Studies 

In  this  section  parameter  studies  are  made  using  the  compu¬ 
tational  model  developed.  Four  parameters,  tilt,  outside  pres¬ 
sure,  surface  roughness,  and  speed  have  been  considered  in  this 
initial  study. 

T i 1 t--T i 1 t  is  a  design  variable  over  which  nearly  complete 
control  is  possible.  Increased  tilt  increases  both  radial  taper 
and  waviness.  The  effects  of  changing  the  tilt  of  a  moving  wave 
gas  seal  from  negative  200  to  positive  800  ym/m  are  illustrated 
in  Figure  4-11.  The  former  value  denotes  a  profile  having  a 
somewhat  divergent  radial  taper  over  most  of  its  circumferential 
period,  while  the  latter  value  denotes  a  profile  which  has  a  pre¬ 
dominantly  convergent  radial  taper.  These  two  extreme  conditions 
are  illustrated  in  Figure  4-12. 

For  negative  and  low  values  of  tilt,  the  gas  film  stiffness, 
the  coefficient  of  friction  and  the  driving  torque  are  relatively 
high,  and  at  the  same  time  the  leakage  rate  is  low  indicating 
that  mechanical  or  asper ity-to-asperity  contact  has  occurred. 

This  is  verified  by  the  minimum  film  thickness  h0  being  less 
than  the  one-half  maximum  roughness  height  c  as  well  as  the  per¬ 
centage  of  the  fluid  film  load  support  being  less  than  100  per¬ 
cent  full  film  load  support.  However,  as  the  tilt  shifts  in  the 


convergent  direction,  higher  average  fluid  pressure  results.  The 
pressure  distribution  across  the  seal  face  has  already  been  shown 
in  Figure  4-10.  As  the  tilt  is  increased  the  seal  faces  even¬ 
tually  separate,  thus,  the  mean  film  thickness  increases  gradual¬ 
ly  and  the  leakage  is  higher.  Thus,  there  is  some  minimum  value 
of  tilt  increase  to  cause  lift  off. 

Outside  Pressure--The  results  of  changing  the  outside  pres¬ 
sure  P0  while  keeping  the  inside  pressure  P^  constant  are 
shown  in  Figure  4-13.  As  the  outside  pressure  increases  from  0.0 
to  1000  psig  (0  to  7  MPa),  the  percentage  of  fluid  load  support 
increases  rapidly  to  100  percent.  At  low  pressure  there  exists 
some  mechanical  asperity  contact  because  the  spring  pressure  is 
greater  than  the  total  lift  force  provided  by  hydrodynamic 
effects.  The  hydrostatic  effect  is  zero  at  zero  pressure. 

Full  film  operation  occurs  when  the  pressure  differential  is 
high.  As  p0  approaches  200  psia,  the  minimum  film  thickness 
h0  becomes  larger  than  the  one-half  maximum  roughness  height 
(c),  and  mechanical  asperity-to-asper ity  contact  no  longer 
occurs.  The  hydrodynamic  effect  is  not  sufficient  to  lift  off 
the  seal  against  the  spring  load  at  zero  sealed  pressure. 

In  the  gas  seal  analyzed  the  velocity  of  gas  discharge  is 
about  500  ft/sec  under  1000  psig  pressure  differential  (the  worst 
case  in  Figure  4-13  is  chosen)  which  is  about  0.5  Mach  number. 

The  method  of  calculating  discharge  velocity  is  shown  in  Refer¬ 
ence  [23].  Therefore,  sonic  discharge  does  not  occur.  However, 
if  the  pressure  differential  keeps  increasing,  sonic  flow  will 
occur  and  the  method  of  analysis  will  become  quite  complex 
according  to  [27]. 

Finally,  the  stiffness  of  the  gaseous  film  is  high  when  the 
outside  pressure  is  less  than  200  psig  (1.4  MPa)  because  of 
mechanical  contact. 
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Figure  4-13.  Effect  of  Outside  Pressure. 
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Surface  Roughness--The  effect  of  surface  roughness  is  shown 
in  Figure  4-14.  The  performance  parameters  are  quite  constant 
when  c  is  less  than  35  yin.  { 0.88  ym).  But  as  c  increases,  the 
friction  begins  to  rapidly  increase.  The  reason  is  that  surface 
roughness  c  is  now  higher  than  the  mean  film  thickness  and  thus 
asper ity-to-asperity  mechanical  contact  occurs. 

Figure  4-15  shows  the  effect  of  surface  roughness  when  there 
is  no  pressure  difference  between  the  two  sides  of  the  seal. 

Now,  as  the  surface  roughness  decreases,  the  minimum  film  thick¬ 
ness  decreases,  the  hydrodynamic  effect  becomes  larger  and  larger 
until  the  seal  lifts  off  (full  film  conditions).  This  is  an 
important  characteristic  of  the  gas  seal  because  if  the  seal  can¬ 
not  lift  off  at  zero  pressure  difference,  high  frictional  heat 
might  occur  at  the  contacting  surfaces  thus  destroying  the  seal. 
The  figure  shows  at  n  =  3  that  the  maximum  roughness  allowable  on 
the  seal  parts  in  order  to  provide  the  95  percent  lift  off  at 
zero  pressure  differential  is  in  the  order  of  c  =  0.10  ym  or 
4  yin.  This  corresponds  to  an  RMS  of  about  1.5  y-in.  In  prac¬ 
tice,  this  surface  roughness  is  difficult  to  achieve  and  main¬ 
tain.  Thus,  it  remains  to  be  seen  whether  or  not  such  a  seal 
could  survive  operation  at  zero  pressure  differential.  Clearly, 
n  =  9  would  work  better  according  to  the  figure.  Also,  a  smaller 
spring  load  would  help. 

Angular  Speed--The  effect  of  angular  velocity  on  performance 
parameters  is  illustrated  in  Figure  4-16.  With  the  tilt  at  5 00 
yn/m,  there  is  no  mechanical  asperity  contact  even  at  zero 
speed.  This  and  the  fact  that  h0  does  not  s ignif icantly  drop 
at  low  speed  means  that  there  is  no  significant  hydrodynamic  load 
support  component.  The  entire  load  is  supported  by  hydrostatic 
effects  due  to  the  radial  taper 

The  coefficient  of  friction  is  approximately  linearly  pro¬ 
portional  to  the  angular  rotational  speed  o>  because  the  viscous 
drag  force  is  proportional  to  angular  speed. 


Figure  4-16.  Effect  of  Angular  Speed 


Potential  Appl ication--The  preceding  results  show  that  a 
moving  wave  gas  seal  could  theoretical ly  operate  at  extremely  low 
friction  and  wear  with  very  low  gas  leakage.  Load  is  supported 
primarily  hydrostatically.  The  main  question  to  be  answered  is 
whether  or  not  sufficient  hydrodynamic  lift  will  develop  under 
zero  differential  pressure  conditions  so  as  to  avoid  serious 
thermal  distortion  and/or  wear  of  the  seal.  Theory  shows  that  a 
flow  surface  roughness  is  essential  to  develop  hydrodynamic  load 
support.  It  is  possible  that  a  carbon  surface  would  polish  to 
the  extent  needed.  Spring  load  could  also  be  lowered.  In  any 
case,  the  answer  to  this  question  must  be  obtained  experimen¬ 
tal  ly. 
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CHAPTER  5 

SEAL  RING  DEFLECTION 


Background 

It  has  long  been  recognized  that  seal  ring  faces  may  deflect 
out  of  their  original  lapped-in  plane  under  the  influence  of 
non-axisymmetric  loads,  non-axisymmetric  section  properties, 
non- ax i symmetric  swelling  of  the  carbon,  or  a  combination  of 
these  factors.  If  such  deflections  are  greater  than  those  which 
can  be  flattened  out  by  the  axial  load  applied  to  the  seal  faces, 
unwanted  gaps  and  leakage  will  result. 

It  is  thought  that  one  of  the  possible  problems  with  present 
submarine  seals  is  that  the  rings  are  too  stiff  to  be  able  to 
close  gaps  created  by  shifting  carbon  segments.  When  gaps  occur 
between  seal  faces,  erosion  continues  to  widen  the  gap  until 
leakage  becomes  excessive.  This  potential  problem  or  some  varia¬ 
tion  of  it  affects  the  design  of  all  seals  including  the  design 
proposed  and  tested  herein. 

Thus,  both  for  the  purpose  of  analyzing  problems  with  pre¬ 
sent  seals  and  for  designing  new  seals,  some  tool  is  needed  to 
predict  non-axisymmetric  deflection  for  seal  rings  and  also  to 
predict  potential  gaps  in  seal  ring  pairs  (the  distribution  of 
face  loading).  For  simple  rings  with  constant  cross  sectional 
properties,  no  splits,  and  simple  harmonic  deflections,  the  face 
loading  distribution  required  for  continuous  contact  is  rela¬ 
tively  easily  obtained  [41-42].  However,  for  more  realistic  con¬ 
ditions  relating  to  the  cross  section  properties  and  more  complex 
distortions,  finding  the  correct  distribution  of  face  loading 
becomes  a  very  difficult  problem. 

The  state  of  the  art  of  predicing  face  loading  in  complex 
cases  is  illustrated  by  a  report  by  Noell,  Rippel,  and 
Niemkiewicz  of  the  Franklin  Institute  [43].  An  evaluation  of 
out-of-plane  seal  distortion  caused  by  nonuniform  joints  was 
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made.  It  is  shown  how  a  nonuniform  cross  section  near  the  joints 
causes  out  of  flatness.  In  the  report  the  contact  between  the 
seal  faces  is  modeled  by  springs  where  tensile  stresses  across 
the  faces  are  allowed.  Since  such  stresses  cannot  occur  in 
reality,  the  computed  results  do  not  predict  how  much  and  where 
the  faces  separate.  Thus,  the  utility  of  the  work  to  date  is 
1 imited. 

To  serve  as  a  useful  tool  in  seal  design  a  more  exacting  but 
general  program  is  needed.  Specifically,  given  two  rings  of 
arbitrary  and  circumferentially  varying  cross  sections  and  load 
contacting  on  a  plane  and  pressed  together  by  axial  loads,  one 
must  find  the  resulting  face  loading  pressure  distribution  which 
arises  as  the  rings  are  loaded  up  and  distort.  At  regions  where 
the  faces  do  not  contact,  the  resulting  gap  must  be  found.  The 
same  program  would  be  capable  of  finding  the  resulting  pressure 
and  gap  distribution  given  faces  originally  distorted  out-of- 
plane  as  well. 

There  are  two  major  features  in  the  needed  model.  In  the 
first  place,  one  needs  a  very  general  tool  to  predict  deflections 
in  arbitrarily  loaded  variable  cross  section  rings.  Then,  the 
two  rings  are  brought  together  to  contact  and  the  pressure  dis¬ 
tribution  is  predicted.  This  second  feature  represents  a  specif¬ 
ic  type  of  contact  problem.  Suitable  models  for  rings  are  not 
available  in  the  literature  and  are  developed  herein.  Several 
special  ring  models  are  developed  and  each  one  is  tested  against 
known  solutions.  The  contact  problem  has  been  treated  previously 
using  finite  element  methods.  A  specific  example  which  illus¬ 
trates  the  technique  is  given  here.  Finally,  it  is  shown  how 
these  models  will  be  put  together  to  find  the  needed  solutions. 
All  of  the  work  is  based  on  the  finite  elements  techniques.  Such 
a  technique  is  essential  for  the  contact  problem. 


Friction  Force 


Figure  5-1.  Mechanical  Loading  on  a  Seal  Ring. 


deflections.  Based  on  the  previous  work  [41-42]  and  neglecting 
all  distributed  loads  (these  will  be  added  later  at  the  nodes), 
the  equations  of  equilibrium  are: 


v; 

+  Ne  = 

0* 

(5-1) 

V 

0 

(5-2) 

N6 

-  Vx  = 

0 

(5-3) 

+  Me  - 

v- 

0 

(5-4) 

V 

RVx  = 

0 

(5-5) 

Me 

-  Mx  = 

0 

(5-6) 

The  above  equations  assume  the  distance  between  the  shear  center 
and  the  centroid  is  negligibly  small.  Combining  the  above,  the 


equilibrium  equations  become 

N“  =  N9  =  0  (5-7) 

M“  +  =  Q  (5-8) 

M' "  +  =  0  (5-9) 

1%  -  Mx  =  0  (5-10) 

Stress  resultant-displacement  equations  are  given  by 

=  Jr  (5-n) 


Mx  *  ^  [UR  -  v")  Jx  -  (u"  +  u)  J  ] 
R  y 


(5-12) 


♦Primes  denote  derivatives  with  respect  to  9. 


I 


1 


( 


I 

I 


[ 

l 

[ 


L 


M  =  E2  [(u  +  u“)  J  -  UR  -  V)  J  1  (5-13) 

y  y  *y 

M0  =  ^  J0[<J>'R  +  V]  (5-14) 

R 

where 

JX=  frAm*1  (5-15) 

Jy-/ar^4rdj  '5-16> 

Jxy  -f  da  (5-17> 


J0  =  torsional  stiffness  constant  such  as  in  Oden  [45]  (5-18) 


Substituting  Equations  (5-11)  through  (5-14)  into  Equations  (5-7) 
through  (5-10)  results  in  the  following  differential  equations 
describing  deflections: 


w’  +  w' ' ' 


J  J 

(u  +  u")  (1  +  -K)  -  (U"  +  U"") 

-,D  C 


aR 


aR 


+  (<j,R  +  <j> '  ' R  -  v'  '  -  v'  '  '  ')  -M  =  0 

aR^ 


(5-20) 


/  GJ9 

t"R  (l  ^ 


-  v  '  '  '  '  +  V  '  ' 


GJ, 


-  ( U "  "  +  u 


. .  \  M  =  0 
x 


(5-21) 
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( u *  +  2u"‘  +  u . )  J  -  U'R  +  4) "  1 R  -  v'"  -  v""1)  J  =0 

y  xy 

(5-22) 

GJ.  J 

U"R  +  v")  -  (*R  -  v”)  +  (u"  +  u)  -j^=  0 

x  x 

(5-23) 

It  can  be  seen  i  f  Jxy  =  0,  then  the  out-of-plane  displacements 
v,  4>  are  decoupled  from  the  in-plane  displacements  u,  w.  The 
above  equations  consist  of  four  differential  equations  in  the 
unknowns  u,  v,  w,  <j>. 

Finite  Element 

Given  the  governing  deflection  equations,  one  can  now  derive 
a  finite  element  stiffness  matrix  for  a  general  ring  element. 

The  stiffness  matrix  is  defined  by 


[F]  -  [K]  [5]  . 


(5-24) 


That  is,  given  an  arbitrary  set  of  displ acmeents  at  the  end  of 
the  element  <5,  what  are  the  forces  and  moments  necessary.  [K], 
the  stiffness  matrix,  defines  this  characteristic  and  is  all  that 
is  required  to  define  a  finite  element. 

There  are  various  coordinate  systems  in  which  the  stiffness 
matrix  can  be  defined.  In  this  work,  a  polar  coordinate  system 
was  chosen.  This  makes  the  stiffness  matrix  invariant  around  the 
ring  assuming  constant  section  properties.  Thus  the  stiffness 
matrix  only  needs  to  be  computed  one  time  for  a  given  element 
size.  It  can  then  be  used  at  any  location  around  the  ring  with¬ 
out  using  a  coordinate  transformation  as  is  required  for  other 
stiffness  matrices  such  as  those  presented  by  Morris  [49], 
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Three  different  stiffness  matrices  will  now  be  derived,  one 
for  in-plane  bending,  one  for  out-of-plane  bending,  and  one  for 
the  coupled  problem.  The  two  simpler  cases  are  used  to  provide  a 
check  on  the  limiting  cases  of  the  coupled  problem  of  interest. 
The  solution  of  the  two  simple  cases  also  helps  in  understanding 
the  solution  to  the  coupled  problem. 

In-Plane  Stiffness--Figure  5-3  shows  an  element  defined  for 
in-plane  stiffness.  For  the  in-plane  case.  Equations  (5-20)  and 
(5-22)  for  JXy  =  0  reduce  to  the  following 

w'  +  w"  1  -  (u  +  u1  ')  (l  +  -  (u‘  '  +  u'  1  ")  i  *  0  ,  (5-25) 

\  ar  /  aR^ 

u'  +  2  u'  "  +  u .  =0  .  (5-26) 

The  objective  is  to  solve  the  above  two  equations  subject  to  the 
displacement  boundary  conditions  shown  in  Figure  5-3. 

The  general  solution  to  Equation  (5-26)  is 

u  =  Cj  +  C2  sin  9  +  C3  cos  9  +  C4  e  sin  8  +  Cg  e  cos  8  .  (5-27) 

Substitution  of  (5-27)  into  (5-25)  and  solving  for  w  gives 


w  =  Cg  +  Zj  sin  9  +  Cg  cos  9  +  +  1  j  9 


C4  9  cos  9  +  Cg  9  Sin  9  . 


(5-28) 


There  are  six  displacement  boundary  conditions  shown  in  Figure 
5-3  and  eight  arbitrary  constants  above.  It  follows  that  there 
are  too  many  constants  and  that  substitution  back  into  the  ori¬ 
ginal  differential  equation  will  eliminate  two  of  the  constants. 

Using  Equations  (5-11),  (5-12),  (5-27),  and  (5-28),  one  may 
show  that 
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,C 

s 

f 


Ea 


Nq  =  ip  [C7  cos  8  -  Cg  sin  8  -  C4  cos  8  +  Cg  sin  9 


-  C2  sin  e 


EJv  , 

— U —  (2C4  cos  9  -  2Cg  sin  8) 
R 


cos  8] 
(5-29) 


Simi 1 arly. 


EJv 

My  =  ^  <C1  + 


2C4  cos  9  -  2C 5 


sin  e) 


(5-30) 


Going  back  to  the  original  equilibrium  equations.  Equations  (5-3) 
and  (5-5)  show  that 

=  -RNg  .  (5-31) 


Substituting  (5-29)  and  (5-30)  into  (5-31)  shows  that 

C8  *  C5  -  C2  (5-32) 

C7  =  C4  +  C3  (5-33) 


so  that  two  of  the  constants  are  eliminated. 

The  boundary  conditions  are  shown  in  Figure  5-3.  Rotation 
represents  the  total  in-plane  rotation  which  is  compatible  with 


_  1  du  , 

“  TTSd  v;  • 


(5-34) 


The  negative  forces  and  moments  shown  at  end  (T)  of  the  beam 
elements  are  in  the  same  direction  as  the  displacements  as 
required  for  the  stiffness  matrix. 
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Using  Equations  (5-27),  (5-28),  (5-32)  and  (5-33)  and 
summarizing 

u  =  C1  +  C2  sin  0  +  cos  9  +  C4  9  sin  9  +  Cj  9  cos  9  , 


♦  =  * 


:i(i 


+1)9+  2C4  sin  9  +  2C5  cos  9  +  Cg 


w  =  Ci(i  * x)  9  -  C2  cos  9  +  C3  sin  9  +  C4(sin  9  -  9  cos 


+  C5(cos  9  +  9  sin  9)  +  Cg 


(5-35) 

(5-36) 


9) 

(5-37) 


Thus,  to  solve  for  the  constants  Ci  boundary  conditions  as 
shown  are  applied  to  Equations  (5-35),  (5-36)  and  (5-37)  as 
follows. 


1 

Sin  9j  COS  9j 

91  sin  9| 

9^  COS  9^ 

— 

0 

r—  -s 

ci 

(SH* 

0  0 

2  sin  9| 

- R - 

2  cos  «j 
- R - 

1 

R" 

-  COS  9j  sin 

sin  9j  -  9^  cos  9j 

COS  9j  +  9  j  Sin  9j 

l 

C3 

sin  cos  9j 

9j  sin  9g 

6j  COS  9j 

0 

C4 

&•■)- 

0  0 

2  sin  »2 
- R — 

2  COS  92 

- R — 

1 

IT 

C5 

L. 

-  cos  9j  sin 

sin  9j  -  cos  9j 

cos  9j  ♦  »2  sin  a2 

1 

C6 

- 

(5-38) 

or 

CAinJ 

[C]  •  C5in]  . 

(5-39) 
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Now,  using  Equations  (5-5),  (5-29),  and  (5-30)  one  may  summarize 
the  forces  as  follows. 


EJ. 


Vx  =  (2C^  sin  9  +  2Cg  cos  9)  , 

R 


EJ. 


My  =  — ^ —  (C^  +  2C^  cos  9  -  2Cj  sin  9) 


EJV 

Ng  =  -ji-  (-2C4  cos  9  +  2C5  sin  9) 
R 


The  generalized  force  vector  is  given  by 


or 


EJ 


^  =  cFin] 


(5-40) 


(5-41) 


(5-42) 


(— 

0 

0 

0 

-Ism 

91 

-|cos 

91 

0 

r- 

C1 

-1 

0 

0 

-  2  cos 

91 

2  sin 

91 

0 

C2 

0 

0 

0 

1  cos 

91 

*  1  s1n 

91 

0 

C3 

a 

0 

0 

0 

1  sln 

92 

|cos 

92 

0 

C4 

1 

0 

0 

2  cos 

92 

-  2  sin 

02 

0 

CS 

0 

L. 

0 

0 

-  Jcos 

92 

1  s1n 

92 

0 

C6 

L  J 

-M 


(5-43) 


(5-44) 


From  (5-39)  and  (5-44) 


■  /  [D1n]  [A^r1  ctn] 


(5-45) 


[K,„]  -  ^  [D,„]  C A | n 3 - ^  . 


(5-46) 


To  provide  an  example  for  future  reference,  the  following 
parameter  values  are  selected 

R  =  3.0  in 


■X**  0.01 

aR^ 


e2  =  o 


02  =  ir/6 


( 

For  this  example  the 

stiffness 

matrix 

[K]  is: 

■"30.72 

21.84 

9.21 

22.00 

20.54 

X 

21.89 

22.80 

-3.24 

-20.54 

11.26 

EJtf  9-21 

-3.24 

62.62 

23.34 

-8.11 

C’Kin^  =  ~Z~ 

n  R^  -22.00 

-20.54 

23.34 

30.72 

-21.84 

• 

> 

20.54 

11.26 

-8.11 

-21.84 

22.80 

^23.34 

-8.11 

-58.84 

-9.21 

-3.24 

(5-47) 


23.34 

-8.11 

-58.84 

-9.21 

-3.24 

62.62 


(5-48) 


The  above  equations  have  been  checked  against  those  given  by 
Morris  [49].  The  stiffness  matrix  is  invariant  with  orientation 
given  a  constant  a2  -  3i  as  it  should  be. 
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Out-of-Plane  Stiffness--Figure  5-4  shows  an  element  defined 
for  out-of-plane  stiffness.  The  moments  and  forces  are  shown 
which  are  compatible  with  the  displacements.  The  stiffness  for 
this  element  is  determined  in  a  manner  similar  to  that  for  the 
in-plane  case.  Letting  JXy  -►  0,  the  out-of-plane  equations  are 
taken  from  Equation  (5-21)  and  (5-23). 


4>'  1  R(1  +  A)  -  v' ' 1  'A  +  v"  =0  , 
( <j>" R  +  v "  )  -  A(*R  -v")  =0  , 


where 


A  = 


EJx 


(5-49) 

(5-50) 

(5-51) 


Elimination  of  <j>  between  (5-49)  and  (5-50)  gives 


+  2v '  "  '  +  v 1  1  =  0 


(5-52) 


to  which  the  general  solution  is 


v  =  +  $2  9  +  B3  sin  9  +  B^  cos  9  +  Bg  0  sin  9  +  Bg  9  cos  9  . 

(5-53) 

Using  (5-49)  to  solve  for  $  gives 

4>R  =  -B^  sin  9  -  B^  cos  9  -  Bg  9  sin  9  -  Bg  9  cos  9 

+  2B5  TT-*  cos  0  '  2B6  sin  9  +  B78  +  Bg  *  (5'54) 


Substitution  into  (5-50)  shows  that 


8-7=0,  8q  =  0  . 


Thus,  the  solutions  for  $  and  v  are  now  given. 


(5-55) 
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Applying  the  boundary  conditions  of  Figure  5-4  to  Equations 
(5-53)  and  (5-54)  gives  the  following  matrix  equations 


9 


-n 

l 

•l 

din  9j 

COS 

»j  Sin  *J 

dj  cos  dj 

0 

1 

T 

ff  coi  8j 

-  ff  s'"  »1 

j  sin  8j  *  j  dj  cos  dj 

|  cos  dj  -  J  dj  sin  dj 

B? 

0 

0 

-  R  5|"  9l 

-  R  c"  91 

'  IT  $'n  *1  *  R  c',,  3 

-  jfi  COS  dj  -  ly-rTT  s'"  \ 

B3 

l 

9? 

sin  »2 

COS  0^ 

d?  sin  dj 

«2  cos  9 2 

\ 

• 

0 

1 

If 

1  cos  oz 

-  V  s,n  ”2 

|  sin  d2  »  ^  o2  cos  $2 

‘cos  »2  •  ^  d2  sin  o2 

i 

B5  ! 

-  | 

0 

a 

-  {  sin  02 

•  ff  co’  "2 

-  £  sin  d2  ♦  |  ^-4-n  cos  d2 

■  ir CM  ’?  •  1  r^nt s,n  92 

BS 

w  J  l 

1  dv 

I  3o  , 


1  dv  i 
S?3o  . 


L  J 


( 


or 


[AQut]  [8]  =  [6nilt]  . 


out- 


(5-56) 


(5-57) 


Using  Equations  (5-4),  (5-13),  and  (5-14)  and  substituting 
for  v  and  $,  one  can  show  that 


V  = 


EJx  82 


R* 

El 


Mx  =  jr  (rri)  ("s5 cos  9  +  86  sin  a) 


(5-58) 


(5-59) 


EJy  /  Bn  28s-  2Bfi  ) 

m9  =  [t  '  r+“A sin  9  ‘  rrr cos  9y 


(5-60) 


Thus  the  generalized  force  vector  is  given  by 
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rhrcos  9 

o 

-  TT7  sin 

‘  i~rT  cos 

[Dout]  [B]  *  [Fout]  • 

From  Equations  (5-57)  and  (5-62) 

tFout]  =  ~^r  tDout^  CAout]  1  £5out^ 
so 

CKout3  =  JT  [Dout]  ^Aout3'1  . 


Using  as  an  example 


R  =  3.0  in. 


A  =  5.0 


9^  =  0 

02  -  ir/6 

The  stiffness 

» 

matrix  for  the  out-of-plane  case  is: 

(5-65 

27.12  -21.20 

2.24  -27.12 

21.20 

-2.24 

21.20  21.55 

2.77  -21.20 

11.75 

-0.15 

2.24  2.77 

1.63  -2.24 

0.15 

-0.92 

[Kout]  = 

-27.12  -21.20 

-2.24  27.12 

-21.20 

2.24 

21.20  11.75 

0.15  -21.20 

21.55 

-2.77 

-2.24  -0.15 

-0.92  2.24 

-2.77 

1.63 

(5-66 

Coupled  Stiffness--Figure  5-5  shows  an  element  defined  for 
both  in-  and  out-of-plane  stiffness.  The  forces  and  displace¬ 
ments  defined  for  the  two  previous  cases  are  shown  together  in 
Figure  5-5. 

To  find  the  stiffness  for  this  element,  the  set  of  coupled 
Equations  (5-20)  through  (5-23)  must  be  solved  subject  to  the 
general  displacement  boundary  conditions  shown  on  Figure  5-5. 
Eliminating  u  and  <j>  between  Equations  (5-21),  (5-22),  and  (5-23) 
an  equation  in  v  only  is  found. 

v .  +  2v' 1 1  '  +  v"  =0  .  (5-67 

As  before,  the  general  solution  is 

v  =  Gi  +  Go  9  +  G,  sin  3  +  G„  cos  9  +  Gr  9  sin  9  +  Gr  9  cos  9  . 


- 

» 

4 

hr 

' 

f 

Eliminating  u  between  Equations  (5-21) 

and  (5-23) 

A 

> 

1 

> 

1 

II 

+ 

(5-69) 

; 

Substituting  for  v, 

1 

<j>' '  +  4>' '  "  =  -2Gg  sin  0  +  2Gg  cos  9  . 

(5-70) 

- 

The  homogeneous  and  particular  solutions  to 

this  equation 

give 

f 

4>R  =  Gy  +  G^  0  +  Gg  sin  9  +  Gg  cos  0 

: 

-  Gg  9  sin  e  -  Gg  0  cos  0  . 

(5-71) 

L« 

Rewriting  Equation  (5-23), 

(u  +  u")  =  4.R  -  v"  -  (^"  R  +  v" 

}  J  » 

(5-72) 

i 

where 

1 

EJx 

(5-73) 

Substituting  for  $  and  v,  the  equation  becomes 

< 

(U  +  U")  ^  =  Gy  +  G15  9  +  [(Gg  +  Gg)  (I  + 

i>  *  2Ge]  sin 

9 

+  C(2g  +  G^)  (1  +  j)  -  2Gg] 

COS  9  . 

(5-74) 

• 

The  homogeneous  and  particular  solutions  for  u  give 

i 

i 
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c: 


( 


to 


•j1  ^  U  =  Gy  +  G^r  9  +  G^q  Sin  9  G^^  COS  9 

°x 


+  ^7  ^G9  +  G4^  ^  +  7^  "  G5^  9  S1’n  9 


i\  (Gg  +  G3)  (1  +  -J)  +  Gg]  0  cos  9 


(5-75) 


Substitution  for  v,  <j>,  and  u  into  Equation  (5-22)  shows  that 


G151 


J2  ' 

xy 

IT 

x  yi 


=  0 


(5-76) 


So  Gls  is  an  extra  constant  and 


G15  =  0 


(5-77) 


Now  substituting  for  u,  v,  and  $  into  Equation  (5-20), 


w'  +  w' 1 


/I  *  iy)  ^ 
'  aR2'  Jxy  aR2 


+  2  [i  (Gg  +  Gg)  (1  +  ^)  +  Gg]  sin  0 
xy 

+  2  -j —  [■£•  (Gg  +  G^)  (1  +  ■£•)  -  Gg]  cos  0  .  (5-78) 
fhe  particular  and  homogeneous  solutions  for  w  give 
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I 


< 


w  = 


G7  +  G13  sin  9  +  ^14  cos  9 


P 


M 


P 


[7  (Gg  +  Gj)  (1  +  -j)  +  Gg]  8  sin  6 

[7  (G9  +  G4)  (1  +  j)  -  Gg]  8  cos  9  .  (5-79) 

We  now  have  the  solution  for  u,  v,  w,  and  <f>  in  terms  of  14 

arbitrary  constants,  two  more  than  dictated  by  the  boundary  con¬ 

ditions  of  the  problem.  To  find  the  extra  constants,  one  must 
check  to  see  that  the  original  differential  equations  of  equili¬ 
brium  (5-1)  through  (5-6)  are  satisfied.  A  combination  of  (5-3) 
and  (5-5)  give 

M'  +  RN'  =  0  .  (5-80) 

y  o 

But  from  (5-11) 

+  RNg  =  Ea[w‘ 1  -  u’]  (5-81) 

so 

w"  -  u'  =  0  .  (5-82) 

Substituting  for  w  and  u  from  Equations  (5-79)  and  (5-75)  into 
the  above  shows  that  G 1 3  and  are  extra  constants  and  can  be 
expressed  in  terms  of  the  other  constants. 

G13  =  [G11  +  7  (G9  +  G4>  ^  +  j)  ~  G5^  ’  (5-83) 

xy 

G14  *  T-  C_G10  "  I  (G8  *■  G3)  <1+^>  -V  '  (5-84) 

xy 
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Thus,  in  summary,  in  terms  of  the  12  remaining  constants, 
v  =  Gi  +  G2  0  +  G3  sin  e  +  G^  cos  e  +  Gg  9  sin  9  +  Gg  0  cos  0 

(5-85) 

<j>  =  ^  {G7  +  Gg  sin  9  +  Gg  cos  9  -  Gg  e  sin  9  -  Gg  9  cos  0}  (5-86) 

Jx  1  1 

u  =  -3 —  {G?  +  G1q  sin  0  +  Gn  cos  0  +  [|  (Gg  +  G4)  (1  +  -J) 

xy 

-  Gg]  0  sin  9  -  [y  (Gg  +  Gg)  (1  +  +  Gg]  0  cos  0}  (5-87) 


w  =  G12  + 


Jxy  aRTJ 


aR 


G7  9 


+  T~  { ^Gll  +  \  (Gg  +  64)  +  ~  Gg]  sin  9 

xy  l 

+  [-  G10-i(G3+G3>  '  G6^  cos9 


[7  (Gg  +  Gg)  (1  +  -^)  +  Gg]  G  sin  0 


[7  (G9  +  G4)  U  +  7)  -  Gg]  8  cos  0  .  (5-88) 


To  satisfy  the  boundary  conditions  one  needs  in  addition 
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J 

J  1 

[( 

L 

aR‘ 

?)  T~ 

’  xy 

xy 

*?\ 

Gy  9 

.  _  u'  w  _  1  L 

*  "  IT  +  IT  '  TT  |  G12 

+  T-  j2^  (G9  +  V  (1  +  ^  ■  G5^  sin  9 


xy 


-  2[j  (Gg  +  G3)  (1  +  +  Gg]  cos  9  j  ,  (5-89) 

=  i-  ^  G~  +  G,  cos  9  -  G,  sin  e  +  Gc(9  cos  0  +  sin  e) 


■R  TT  |  2  3 


^  J  Ml  V 

+  Gg(-  9  sin  9  +  cos  9)|  ,  (5-90) 


Applying  the  boundary  conditions  shown  in  Figure  5-4  gives 
the  matrix  Equation  (5-91)  (next  page)  which  is  equivalent  to 


^coupled-^  ^  ^coupled^ 


(5-92) 


To  get  the  corresponding  forces,  substituting  for  u,  v, 
and  w  into-Equations  (5-11)  through  (5-14)  and  (5-3)  and  (5-4) 
gives 


y— j  (2Gg  cos  0  +  2Gg  sin  9) 


+ 


-  (Gg  +  G^)  sin  9  +  (Gg  +  Gg)  cos  9 


I 

< » 


(5-93) 


V 


y 


EJxG2 


(5-94) 
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Equation  (5-91) 


EJ  (/J  J  \ 

~£T  | V3^  "  /  (2Ge  Sln  9  "  2Gs  C0S  9) 

i}  +  [(69  +  S^)  cos  0  +  (Gg  +  Gg)  sin  0]  j  , 

xy  x  j  1 


Mx  =  — y- J  [(S3  +  Sg)  sin  0  +  (G^  +  Gg)  cos  0] 
R 


My  =  jjlT  j  ("3^  "  T~j  {G7  "  2G5  C0S  9  +  2G6  sin  0) 


(5-95) 


(5-96) 


♦  i 


M)- 


Jxvl  ) 

X1  [(Gg  +  G4)  cos  0  +  (Gg  +  G3)  sin  0]  j  . 

(5-97) 


Mq  =  — j-  [^2  +  (Gg  +  Gg)  cos  0  -  (G4  +  Gg)  sin  0] 

R 


(5-98) 


Using  the  above,  the  generalized  force  vector  is  given  by 
Equation  (5-99)  (next  page).  This  equation  has  the  form 


[Gcoupled^  [G-^  ~  [^coupled-^ 


(5-100) 


Using  Equations  (5-92)  and  (5-100) 


[^coupled-^  r2  [^coupled-^  [\oupied^  [5coupled^  -(5-101) 


Thus , 
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r 

i 


>s> 


C 


b 


< . 


b 


£JX  _! 

^coupled^  "  02  ^coupled^  ^coupled-^ 


As  an  example,  let 
R  =  3.0 

91  =  0 

@2  =  it/  6 

A  =  5 


(5-102) 


(5-103) 


-X-  =  0.01 

aR^ 


-jX  =  0.1 

Jx 

^■o-1 

For  this  case  the  stiffness  matrix  [KC0Upiecj]  Is  calculated  as 


t*, 


coupled 


.  “« 
]'7T 
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30.71 

2.64 

9.21 

2.03 

21.84 

0.17 

-21  99 

-2.64 

-23.33 

2.12 

20.53 

-0.19 

2.64 

27.12 

-0.71 

21.20 

2.12 

2.24 

-2.64 

-27.12 

-0.71 

21.20 

2.12 

-2.24 

9.21 

-0.71 

62.6  2 

-0.72 

-3.24 

-0.09 

23.33 

0.71 

-58.83 

-0.39 

-8.11 

0.00 

2.03 

21.20 

-0.72 

21.55 

2.16 

2.77 

-2.12 

-21.20 

-0.39 

11.75 

1.17 

-0.15 

21.84 

2.12 

-3.24 

2.16 

22.79 

0.28 

-20.53 

-2.12 

-8.11 

1.17 

11.25 

-0.01 

0.17 

2.24 

-0.09 

2.77 

0.28 

1.63 

-0.19 

-2.24 

-0.00 

0.15 

0.01 

-0.92 

-21.99 

-2.64 

23.33 

-2.12 

-20.53 

-0.19 

30.71 

?.  64 

-9.21 

-2.03 

-21.84 

0.17 

-2.64 

-27.12 

0.71 

-21.20 

-2.12 

-2.24 

2.64 

27.12 

0.71 

-21.20 

-2.12 

2.24 

-23.33 

-0.71 

-58.83 

-0.39 

-8. 11 

-0.00 

-9.21 

0.71 

62.62 

■  0.72 

-3.24 

0.09 

2.12 

21.20 

-0.39 

11.75 

1.17 

0.15 

-2.03 

-21.20 

-0.72 

21.55 

2.16 

-2.77 

20.53 

2.12 

-8.31 

1.17 

11.26 

0.01 

-21.34 

-2.12 

-3.21 

2.16 

22.79 

-0.29 

-0.19 

-2.24 

0.00 

-0.15 

-0.01 

-0.92 

0.17 

2.24 

0.09 

-2.77 

-0.28 

1.63 

_/ 

(5-104) 

k_ 


Verification  of  Stiffness  Matrices--The  correctness  of  the 
preceding  stiffness  matrices  has  been  verified  in  detail  by 
Knowlton  [24].  He  assembled  the  elements  into  rings,  performed 
solutions  and  compared  to  known  results.  The  in-plane  results 
were  compared  to  those  of  Morris  [49]  and  straight  beam  theory. 
There  is  some  slight  difference  when  compared  to  Morris  because 
he  uses  straight  beam  theory.  Comparison  is  perfect  of  course  as 
R  -*■  ®. 

The  out-of-plane  case  was  verified  by  comparing  to  a  problem 
solution  for  out-of-plane  deflection  presented  by  Roark  [44]. 
Agreement  was  perfect. 

Finally  for  the  coupled  problem,  the  only  completely  inde¬ 
pendent  solutions  formed  are  those  by  Meek  [47].  Comparison  was 
made  to  a  cantilever  curved  beam  section  loaded  in  two  planes  and 
with  JXy  *  0.  Again,  agreement  was  perfect.  Coupled  numerical 
solutions  were  also  compared  to  those  analytical  solutions  pre¬ 
sented  by  Lebeck  [41].  Again  agreement  was  perfect. 

Thus,  based  on  these  limited  checks  the  preceding  results 
appear  to  be  correct.  The  reader  is  again  cautioned  that  warping 
is  not  included  here  and  this  is  valid  for  compact  sections. 
However,  warping  does  play  a  major  role  in  many  practical 
problems  and  must  be  included  where  necessary. 

General  Contact  Problem 

Given  the  stiffness  matrix  results  for  a  ring,  the  next  step 
in  the  solution  of  the  general  seal  ring  contact  problem  is  to 
create  a  model  of  contact  and  incorporate  it  into  problem  solu¬ 
tions.  Before  illustrating  how  this  is  to  be  done  for  the  actual 
seal  assembly,  a  simple  example  will  be  used  to  show  how  such  a 
contact  problem  can  be  solved. 

The  example  problem  will  start  with  a  beam  on  an  elastic 
foundation  (Figure  5-6).  This  problem  can  be  modeled  using  beam 
finite  elements  (Figure  5-7).  The  stiffness  matrix  for  a  beam 
element  as  given  by  Irons  [51]  is 


\ 


Figure  5-6.  Beam  on  Elastic  Foundation. 


Figure  5-7.  Beam  Element. 
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where  the  displacement  and  force  vectors  are 


[V 

Pa] 

[6]  = 

9A 

(5-106) 

[F]  = 

ma 

WB 

PB 

^9B> 

AJ 

(5-105) 


(5-107) 


For  a  solution  by  the  finite  element  method,  one  requires  an 
expression  for  the  total  strain  energy  of  the  system.  The  strain 
energy  for  one  beam  element  and  the  springs  under  the  element  is 


i(e)  =|[s]T  oy  [si 


1  kA 

1  T 


(wA  '  A0A)2  +  \~r  (wB  ■  A0B)2  • 


(5-108) 


The  terms  k/\/2  and  kg/2  assume  the  spring  under  each  node  is 
divided  into  two  equal  parts,  one  part  associated  with  each  adja¬ 
cent  element. 

Total  potential  energy  is  the  strain  energy  over  all  of  the 
elements  plus  that  due  to  the  potential  of  applied  loads  or 

U  =  £u(e)  -  [F]T  [5]  .  (5-109) 

e 


According  to  the  principle  of  minimum  potential  energy 

SU  =  0  .  (5-110) 

Applying  this  principle  to  Equations  (5-108)  and  (5-109)  combined 
gives 
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jLC<J  +^|7“  (wA  "  &0AJ  +  T  (wB 


-  An.)  .  C«3  -L CFJ 


(5-111) 


In  matrix  form,  using  the  nunbering  system  shown  in  Figure 
(5-6),  the  system  of  equations  to  be  solved  is: 


^‘r-rr 


St  -12 


-St  24  *  k. 


-12  St 


-St  24  ♦  k. 


-12  -St  24  ♦  k 

St  2t*  0 


-12 

St 

-St 

a2 

,3 

0 

0 

8t2 

(5-112) 

System  of  Equations  (5-112)  can  be  readily  solved  for  the  equili¬ 
brium  displacements  w-j  for  arbitrary  kj  and  .  The  a, 
allow  for  any  preload  of  the  springs.  An  algorithm  solving  Equa¬ 
tion  (5-112)  was  written  and  solutions  to  various  beam  on  elastic 
foundation  problems  were  carried  out.  For  the  case  of  a  semi¬ 
infinite  beam  on  elastic  foundation  loaded  at  the  end  only,  the 
displacement  at  the  end  was  found  to  be  within  9.2  percent  with 
20  elements,  2.5  percent  with  40  elements,  and  0.4  percent  with 
30  elements  of  the  exact  solution  results.  Thus  the  method  was 
verified. 


Now  to  solve  a  contact  problem  one  need  only  make  the 
springs  stiff  enough  to  approximate  contact  stiffness,  allow  for 
preload  by  adjusting  the  A-j,  and  allow  for  noncontact  by  caus¬ 
ing  the  k-j  -*•  0  where  appropriate.  The  example  problem  to  be 
solved  is  shown  in  Figure  5-8.  Here  a  beam  is  brought  into  con¬ 
tact  with  two  circularly  curved  surfaces  of  high  rigidity.  Loads 
Pi  and  P2  act  to  deflect  the  beam  making  it  contact  the  two  sur¬ 
faces  with  some  knowable  contact  pressure  distribution.  The 
problem  to  be  solved  is  to  find  this  pressure  distribution,  or 
more  precisely  to  find  which  points  in  the  discretized  problem 
make  contact  and  at  what  load. 

The  procedure  used  to  solve  the  stated  problem  is  as 
follows: 

1)  Pick  initial  contact  points  A  and  B  and  set  k/\  and 
kg  equal  to  contact  stiffness.  Set  all  other  ki  -*• 

0. 

I  I 

2)  Set  P2  =  nF  x  Pi  and  P2  =  nF  x  P2  where  F  is  some  frac¬ 
tion  of  1.0  and  n  =  1  (F  =  0.1  appears  to  work).  Solve 
for  the  deflection  using  this  reduced  load. 

3)  Compare  the  deflection  to  the  known  gap.  If  deflection 
is  greater  than  gap  at  i,  set  k,  equal  to  contact 
stiffness  and  set  Ai  also. 

4)  If  deflection  is  less  than  A  at  i,  the  set  kj  *  0. 

5)  Repeat  steps  2),  3),  and  4)  with  n  =  2,3  •••  1/F. 

6)  If  assumed  contact  pattern  agrees  with  solved  for  con¬ 
tact  pattern  at  nF  =  1,  the  solution  is  consistent  and 
in  hand.  If  it  does  not,  then  set  F  equal  to  a  finer 
increment  at  least  near  full  load. 

The  above  procedure  was  found  to  give  a  unique  solution  to 
the  example  problem.  It  was  found  that  only  by  slowly  increasing 
the  load  was  it  possible  to  arrive  at  a  consistent  solution  (see 
6)  above).  Applying  full  load  resulted  in  continual  changes  in 
contact  pattern  with  each  iteration--some  of  the  time. 


For  the  example  problem. 


i  =  1  in.  (total  length  20  in.) 

E  =  30  x  106  lb/in.2 

k.  =  100  x  106  lb/in. 

Px  =  20000  lb  (5-113) 

P2  =  10000  lb 
R  =  100  in. 

I  =  0.085  in.4 

A.j  -  based  on  circular  support 
The  computed  result  is  shown  in  Table  5-1. 

Ring  Contact  Problem 

The  contact  problem  for  two  rings  is  considerably  more  dif¬ 
ficult  than  that  for  the  beam.  Figure  5-9  shows  how  the  problem 
can  be  set  up.  The  two  seal  rings  are  pushed  together  by  arbi¬ 
trary  loads.  ks  represents  the  stiffness  of  the  support  and 
kf  is  the  stiffness  of  the  contact  at  the  face.  The  contact  is 
represented  only  by  one  set  of  springs  at  the  mean  contact  rad¬ 
ius.  Additional  sets  could  be  included  to  better  model  face 
contact,  but  it  is  thought  that  for  initial  analysis  at  least, 
one  set  is  sufficient. 

The  solution  to  the  problem  represented  in  Figure  5-9  is 
accomplished  in  much  the  same  way  as  that  for  the  beam.  First, 
using  the  stiffness  matrices  derived  for  the  ring  elements,  the 
strain  energy  of  the  entire  assembly  including  the  springs  is 
found.  Then,  derivatives  are  taken  and  the  set  of  resulting 
equations  solved. 

These  final  steps  remain  to  be  completed. 


Table  5-1 


Computed  Results  -  Beam  Contact  Problem 


Node 

w 

i 

A 

i 

Contact 

Force 

-0.1141 

-0.1251 

0 

-0.07761 

-0.08003 

0 

-0.04505 

-0.04501 

3459 

-0.02007 

-0.02000 

7014 

-0.005011 

-0.005005 

602 

-0.00008922 

-0.0 

8922 

-0.004347 

-0.005005 

0 

-0.01446 

-0.02000 

0 

-0.02652 

-0.04501 

0 

10 

-0.03662 

-0.08003 

0 

11 

-0.04085 

-0.1251 

0 

12 

-0.03662 

-0.08003 

0 

13 

-0.02652 

-0.04501 

0 

14 

-0.01446 

-0.02000 

0 

15 

-0.004347 

-0.005005 

0 

16 

-0.00008922 

-0.0 

8922 

17 

-0.005011 

-0.005005 

602 

18 

-0.02007 

-0.02000 

7014 

-0.04505 

-0.04501 

3459 

-0.07761 

-0.08003 

0 

-0.1141 

-0.1251 

0 
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Two  Ring  Contact  Problem. 


CHAPTER  6 

SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


Nine-Wave  Seal  Design 

As  of  Test  117,  nine-wave  seal  performance  began  to  approach 
expected  performance  and  the  various  design  problems  appear  to 
have  been  solved.  Conclusions  concerning  the  long-term  perform¬ 
ance  must  await  the  completion  of  the  scheduled  2000  hour  test. 
However,  several  conclusions  can  be  reached  concerning  the  design 
itself. 

1)  As  designed,  54  thin  arms  must  be  machined  out  of  a 
solid  Inconel  625  piece.  While  there  was  success  in  accomplish¬ 
ing  this  difficult  machining  task,  a  great  amount  of  NC  machine 
time  was  required.  A  simpler  design  or  a  more  easily  worked 
material  must  be  fc-nd  for  practical  application  of  the  concept. 

2)  As  designed,  many  1/32-in.  holes  were  needed  in  the 
Inconel  625.  Again,  while  there  was  ultimately  success  in  dril¬ 
ling  these  holes,  a  great  amount  of  machine  time  was  required. 
Special  tech'iiques  had  to  be  developed  to  remove  broken  drills. 
This  problem  could  probably  be  overcome  by  using  an  EDM  machine 
to  make  these  holes.  Alternatively,  more  easily  worked  material 
or  a  design  which  does  not  require  so  many  small  holes  in  a  hard 
to  work  metal  is  to  be  preferred. 

3)  Larger  than  expected  arm  deflections  led  to  hydraulic 
leakage  problems.  A  3-D  stress  analysis  of  the  configuration 
would  have  shown  this  problem  or  a  simpler  more  readily  analyzed 
configuration  should  be  used. 

4)  Ring  stiffness  with  regard  to  waviness  turned  out  to  be 
larger  than  expected  because  significant  torsional  warping 
stresses  developed.  Again,  a  3-D  stress  analysis  would  have 
shown  this  problem  earlier  and  should  have  been  used.  This  prob¬ 
lem  was  overcome  by  reducing  the  cross  section. 

5)  A  method  was  developed  to  retain  all  0-rings  on  assem¬ 
bly.  This  should  have  been  incorporated  in  the  original  design. 
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6)  A  large  thickness  epoxy  bond  of  the  carbon  into  the 
Inconel  allows  the  carbon  to  creep  such  that  it  need  not  conform 
to  the  desired  waviness.  This  problem  was  eliminated  by  reducing 
the  bond  thickness  to  0.001  in.  or  less. 

7)  To  eliminate  hydraulic  leakage  in  spite  of  the  large  arm 
deflection,  a  piston  arrangement  was  used  which  gives  the  0-ring 

a  zero  gap  through  which  it  can  extrude.  After  trying  several 
other  approaches  it  was  concluded  that  if  there  exists  any  seal 
ring  gap  which  must  be  sealed  by  the  0-ring,  the  periodic  motion 
of  the  arms  will  encourage  the  0-ring  to  extrude  through  the 
gap--even  though  the  gap  size  is  perfectly  acceptable  for  a 
static  0-ring  application. 

Test  Machine  Modification 

The  test  machine  was  modified  so  that  it  can  automatically 
run  a  seal  test  consisting  of  variable  speed,  pressure,  and  mis¬ 
alignment.  The  test  becomes  completely  computer  controlled. 

These  changes  required  considerable  modification  of  the  test 
apparatus  and  its  associated  electronics.  While  existing  soft¬ 
ware  was  modified  to  run  the  2000  hour  test,  it  became  clear  that 
a  new  software  package  incorporating  greater  flexibility  and 
refined  logic  should  be  written.  The  computer  controlled  test 
cycle  works  as  planned. 

Seventy-Five  Percent  Balance  Ratio  Test 

A  500  hour  75  percent  balance  ratio  flat  face  test  was  run 
to  try  to  duplicate  and  explain  some  previous  data.  It  was  con¬ 
cluded  that  the  seal  can  start  up  in  one  of  at  least  two  modes  of 
operation  and  that  the  second  mode  of  high  friction  zero  leakage 
is  not  stable. 


Nine-Wave  Seal  Design  Methodology 

Chapter  3  points  out  the  design  criteria  for  the  nine-wave 
seal.  Many  factors  were  considered  including  zero  moment  design, 
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location  of  the  centroid,  low  stiffness,  and  corrosion  resis¬ 
tance.  A  good  part  of  the  design  was  accomplished  using  a  compu¬ 
ter  program  which  relates  the  various  geometrical  variables  to 
the  design  criteria.  Alternative  conf igurations  were  evaluated. 
Finite  element  and  finite  difference  stress  and  deflection  calcu¬ 
lations  were  made.  Performance  was  predicted  using  previously 
developed  models.  A  clear  design  methodology  has  been  set  forth 
which  should  be  generally  applicable  to  wavy  seal  design  for 
other  purposes. 

Misal ignment 

Since  submarine  seals  undergo  misalignment  as  a  submarine 
decends,  a  condition  for  the  simulated  submarine  operation  test 
is  to  create  a  similar  variable  misalignment.  This  was  achieved 
in  a  unique  way  for  the  wavy  seal  tests.  Using  the  worm  gear 
drive  seal  support  originally  used  to  move  waviness  and  machining 
angles  into  the  parts,  the  amount  of  misalignment  can  be  adjusted 
*1  _  as  needed.  This  saves  a  lot  of  time  in  not  having  to  remisalign 

the  machine  by  hand  for  each  test  condition.  The  misalignment  is 
computer  controlled. 

\  Moving  Wave  Gas  Seal 

The  complete  analysis  of  a  moving  wave  gas  seal  is  included 
in  Chapter  4.  Certain  advantages  of  a  moving  wave  gas  seal  com¬ 
pared  to  other  types  are  pointed  out.  The  analysis  uses  the 
alternating  direction  implicit  method  to  solve  the  compressible 
flow  equation.  Other  aspects  of  the  analysis  are  similar  to 
those  already  used  for  liquid  seals.  Parameter  studies  are  made 
which  show  that  at  low  speed  operation  the  fluid  pressure  left  is 
predominately  due  to  hydrostatic  effects.  At  zero  pressure  dif¬ 
ferential,  a  very  smooth  surface  is  required  to  allow  the  seal  to 
hydrodynamical ly  lift  off.  The  design  theoretically  appears  to 
,  work  and  offers  a  gas  seal  with  very  low  friction  and  leakage. 

Experiments  are  needed  to  confirm  these  predictions. 
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Predicting  Seal  Ring  Gap 

To  address  the  difficult  problem  of  predicting  seal  ring  gap 
(or  zero  gap  as  is  desired)  when  rings  are  not  axisymmetric  or 
when  seals  have  complex  out  of  flatness,  two  significant  develop¬ 
ments  were  made  which  will  lead  to  the  development  of  a  general¬ 
ized  computer  program  for  predicting  seal  gap. 

1)  A  new  ring  finite  element  was  developed  for  the  case  of 
coupled  in-plane  and  out-of-plane  loading  and  deflection.  The 
element  has  been  checked.  It  provides  a  relatively  simple  method 
to  solve  the  contact  problem  compared  to  3-D  solid  elements. 

2)  Contact  problems  for  beams  otherwise  similar  to  the  one 
of  interest  have  been  solved  using  an  iteration  method. 

After  1)  and  2)  are  combined,  the  analysis  of  actual  seal 


rings  can  be  initiated. 
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Figure  A-l.  Test  #111,  75%  Balance,  Parallel  Faces,  500  psi,  1800  rpm. 
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*IY*AlY*l./l2,»X6»*i»IY5**9-Y4J*A4*Ul*X7*Xb/2.-X3*4)«*2 

Alf*AlY*l./l2.»lX5-<4)»nJ*(YS-Y4} 
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PREVIOUS  PAGE 
IS  BLANK 


30660 
30673 
30663 
30693 
30703 
30713 
30733 
30733 
307  4Q 
33753 

30  763 
33773 
30763 
30793 
30603 
33313 
33833 
30833 
33643 
33853 
33363 
30973 
30633 
33393 
30903 

30913 

30923 

30930 

33940 

30953 

30963 

33973 

30980 

30990 

31300 

31310 

31323 

31333 

31350 

31360 

31370  3 

31393 

31390 

31100 

31113 

31123 

31133 

31140 

31153 

31160 

31173 

31190 

3H9Q 

31203 

31213 

31220 

31230 

31243 

31250 

31263 

31273 

31263 

31293 

31300  : 

31313 

31  323 


AIY»AIY*A5»(X1»X6*X7»CXS*X4)/2,-XBAR}a*2 

AiriAlY*l(/l2.*(X4»X66}**2*f5*A6«U2»(3«'((4»X66)/2fA9A4Jt*2 
AlY*AlYM./12.»X65»»J*{r5-rS5}*A7*CX2»X3*<U»X66/2,»XSM)»*2 
110  f»i.*eCC-YaAR 
AJR»AXX*AXY 

Fl»2,*PI»(RLFX5*Xb) »Y  7 
F2*PI»  CR30*»2-R3<*»2) 

F3»2,*Pi»(Raa-xn«Yi 
F5*2,»PI«(R33) *72 
Fb»PI*( (R00)«»2»(RI*X*3)*«2J 
F7«2,«PI*IRI*XX3) «YY3 
F3«FJ.F».f6 
A«U»-Fl*(Y9AR-r  7/2. J 

R2*2,/3,»(R33**3«R3»*3J /  (  R0a**2-Rd* • 2) 

Al2*-f2*(R2.RC3 

A'l3«F3*(Y6AR«Yl/2,J 

Ai5»-f5«(Yl*Y2/2,«Y3ARl 

R6«2,/3.*(CRa35»»3"(Rl*XX3)»*3)/tCR66)**2»lRI*XX3j»»2J 

A1b*Fb«(Rb-RC) 

A97»«F7»(Y1*Y2*YY3/2**Y34RJ 

R9a2,/3.»XX3*RI 

A48*F8*(R8»RC) 

T4«XAI,11*AM2*Art3*4'H*4'15*A*b*A!i7*AI961  7  (2,*»1»RC) 
RaTAT£*TR»RC*»2/E/AlX*R4 

140  WRITE  (6,  ISO'  XBaR|  Y  3  A  9  ,  AI X , AIT , AJR, A J 3 , A, 4*3 , V , ECO ,  RC,SIX,R3rATE 
1  30  F3RMATC/,5X,  'X3AR»*,F7,4,/,3X,  * Y BAR* • , F7 , 4 , / , 5X ,  «  I X  *»  ,  £  J  2, 3,  /  ,  5*  , 

•  ,IYi»lei2t3,/,Sx,*JR«,.cl2.5,/,5X,»JrHeTA«*,ei2.S,/,9X, *A«',E12,S 

•  , /,5X, 'MTHETAa',E12.5,/,5X,'V»*,ei2.3,/,5X, *ECC**,:i2.5,/f 3X, 

•  *RC*»,F7,4,/,5X,'iIX*,,*12.i,/,SX,»RaTATI3.'l»»,ei2,a) 
wRITE(6,160)  Ya»MI5,vaPiI,AIXY,Xl,X2,X3,X6,X5,X6»X7 
RRITE16,170J  ri,Y2,Y3,r4,Y6,Yb,T7,R3a,RPX,Pl.,RI,R|.. 

160  F0R4AT(5x,'V/PRI3*',E12,5,6x,  * V/PRI •• , £  1 2« 5 , / , 5X , ' I X Y* ' , Ei 2, 5, 

•  /,5Xr 'Xt«»lF6.4,2X#*X2«<',F6,4,2X,*Xi««,F6t4,2x#  *X4*»,F6t«,2X, 

«  'XS«*,Fb,4,2X,*Xb**,F6,4,2X,'X7»*,Fb.4J 

1 7  o  igRMAT(5X,*YI**,F6,a,2X,  '  Y 2 » •  ,  F 6 , 4 , 2 X , » Y 3  * • , F 6 , 4 , 2X • ,Y»«»,F6,» 

•  , 2X, 'Y5«*,Fb.4,2X, *Yb**,Fb,4,2X, *Y7*',Fb,4,/, 5X, * R33* * , F6, 4 , 2* , 

•  *RPI»Y,F6,4,2X,*PU»'#^5,a,2X,‘RI»'f F6,4,2X,  'RL»*iF6,4,///j 
G3  T 3  1 3QU 

£90 

03UdL£  PRECISI39  FUNCTI3N  FliCRU) 

C3MH0N  RC,X4,X5,X7,Y5,Xl,X2,X3,A6,l66,Yl,r2,r3,Y4,rSS 
Ca.-'naN  Yb,Y7,3,3HI,AK,RI,R9,XXJ,YY3,ri,T2,AJR 
C3MM3N  Al, A?, A3,A4,A3,A6, A/, Xd Art ,Y3AR, *33, RPX,P3 
C3*M0N  PI,  AMg,Pw,  e«  J,  V3PRIG,  V3PHI,£CC,  A,  A  13,  AJ3,  Y,  CU,  A  I  X  ,  Y 

oauaus  ?reczsi3n  x4,xs, xt,y5,rc 

OIRE»SIaN  *(30) ,Y (53) ,P  (53# S3) ,41*1  133,50) ,AIP1 (53#  53 J 

0 1 *£NS 1 3N  A jPi (50,53) , A J*1  t 33,60 J ,3915150, 53 J 

X4»RI-TI-«l 

X5*R8*T2-RU 

X7«»X1+X2YX3*X4.X3«<6 

Y9»Y1*y2»Y3*Y4»Y6*y7 

Al»*l»Y2 

A2»(X2-*n«(Yi*Y2J 

Ai«(x7*xi-X2)»(Yl*Y2-rJ) 

A4*X6» (Y5*»8«raj 
A5»IX5-*4J  »(Y5.Y4J 
AbHX4-*boJ  » y 5 
a7»X66»(Y5-y55J 

X3AR«Al*Xl/2,*A2*(Xl*tX2-XlJ/2,J*A3«U2YCX7*Xl-*2J/2,) 

X3AR«X3AR*A<»»(Xl*x7*Xb/2.)*A8*(Xl*Xb*xf*(X5-X4)/2tJ 

XiArt*xaArt*Ab*(X2«X3*(X4-»bb)/2,J*47»U2*Xi*«4-Xb6/<,J 

X3AW«X9AR/(Al*42*A3*4«w»5»Ab»A7J 


Y3AR*Al«CYl*Y2/3.J*A2»vrl*Y2)/2,»Ai»lM*»2-Y3)/2. 

Y3AR»T9AR*A<»»lY7*CYb*Y5*Y4J/2.j*A5«Xrb*Y7f(Y5«*4J/2,l 


31333 

Jl  383 

31353  C 

31363 

31373 

31383 

31393 

3U03 

31813 

31  823 

31833 

31840  C 

31850 

31860 

31870 

31880  3 

31890  3 

31500 

31513 

31523 

31530 

01540 

01550 

01563 

31573 

01580 

01590 

31500 

31610 

01520 

01530 

31640 

01550 

01560 

31670 

01680 

31590 

31700 

31713 

01720 

31730 

01740 

31750 

01760 

01773 

31  790 

31790 

01803 

31810 

31620 

31530 

31880 

31550 
31860 
31870 
31883 
31590 
01900 
31913 
31  823 
31930 
31  98  3 
31853 
3H60 
318  7  3 
31  880 


73  AS  *'?8A9  +  *6»  (7  5*7  7*75/5.  j*A/»(Y55*r6*77*(7  5- 7551/2,1 
7  3AR*YSAS/ (Al*A2*A3*A8*45*A6*A7) 

AIX«t ,/ia,»Xi*72*»3*Al«(ri*72/a,-73A8) **a 

AIX»AIX*l,/l2,*CX2-XU*l71*72}«*3*A3.(t7l*72)/g,-7SA<l*«2 
AIX*AlX*l,/ia,»(X7*Xl-X21 *(71*72*73) •*3*A3»((71*73«r 5)/ 

*2 

AIX*AIX*l,/12,*Xb* (75*76-781 •• 3* A«*U7*(Y5*Y5-Y8! /2, •rttAS}**a 
AIX*AIX*t,/ia,*(X5-X4)«(75-78)*«i*A5*(7b*T7*(75-781/a,-f3ASl»*2 
AjX»AlX*l,/ia,»(X8-X661*75**5*A6*(76*f7*7i/a,-73A81»«a 
AIX»AXX*t,/l2,»*65*(75-7551**3*A7«(755*76*77*(75-755l/a,*7iASl*»a 

saa«x2*x3*X8*su 
—  S3I*RU*X5*Xb*x7/2,  «4— 

R0*saa-x8AS 


EIX«e*AIX 
£89*0,000001 
AAl 5X66 
AA2»X8 

A  a3*X5 
A A4*  X  5*X  b 
AaS*X5*x6*X  7 
AA6*Xl«XS*Xb*X7 
81»77 
92*7  1 
93*76*77 
98*71*72-73 
95*755*76*77 
95*71*72 
11*8 
12*13 

13*17  n 

18*19 

15*20 

15*22 

31*17 

32*23 

33*24 

38*26 

35*30 

36*37 

X  X  *0, Q 

X(il *0,0 

OX *61/(11-1) 

03  IQ  1*2,16 

:*ci»cT,ii)  ox*(B2-ai)/(ia-u) 

17(1, 5T, 12)  OX* (33*92) /(I  3-12) 

17(1,  GT, 13)  0X*(84-93)/C4-t3) 

17(1, 37,18)  OX* (95-98) / (13-18) 

17(1,37,15)  0X»(3b»35)/(I5-15) 

XX*XX*0X 
X (I) *XX 
13  CONTINUE 
77*0,0 
7 (1) *0,0 
07«*A1/ (31-11 
03  20  3*2,35 

IF(j,Gr,3l)  07 *(A a2-a*U / (32-311 
17(3,07,321  07*(4A3-A*2l / (33-321 
17(3,37,331  07* (AX8-AA3) / ( J  8-331 
17(3,37,38)  0  7 , l 4 45-4 A4) / (35-J8) 
:rCJ.3T,J5J  07*(4*b-445)/(J5»J5) 

77*77*07 
y  f  J  J  ITY 


‘ 

t 

• 

r 

31990 

20 

continue 

32000 

l4AX*I6-l 

32010 

J84X»J6«1. 

32320 

00  30  I*2,IMAX 

32330 

03  30  Ja2,JMAX 

32340 

AiJ*42t/(x(x-n-xu))/cxci*u-xcin*a,/(UJ-i)-r<jn/(TUM  wcjj 

32350 

ft 

) 

$ 

32360 

AiiiUiJW./(xci.n-xcin/cxci-u*xii*i))/AU 

32370 

*121(1,  J)»2./(X  (1*11  -X  (1)1/  (XU-1) -<<m  J3/AIJ 

32360 

*JKl(I.J)*2./CT(J-l)-T(J))/CUJ-l)-»(JM5  J/4IJ 

32390 

*jPici,j)»-2,/(rcj*n-TU)  j/cru-n-f  (j*iij/Au 

32100 

89MSCX, J)#2./AIJ 

32113 

30 

C3NTINUE 

32120 

oa  40  i a i « x 6 

32130 

03  40  3*1,36 

. 

32140 

40 

2(1,3) aO. 0000003030001 

i 

32150 

aieaai ,5 

32160 

a^aai,-aiEG 

32170 

oa  60  IT»1,200 

32190 

A  4AX  *Q  ,  Q 

32190 

oa  SO  J  a 2 , J  4 A X 

32203 

IMAX*X6-1 

32213 

HIN»I5*1 

32220 

IF(3,GT,3t)  1 1 1 N*  X  3* 1 

• 

32233 

IF (3.GT.32)  I4AX  a  1 4» 1 

k 

32240 

IFCJ,GT,33)  XHNallM 

* 

32250 

IffCJ.GT,j4)  I4I4«2 

32263 

IF(3,GT.35]  I-inaia*! 

. 

32273 

IF  CJ.GT,  J5]  I4AX*X6«l 

• 

32250 

03  50  lalHlN.IMAX 

- 

02290 

23L0«P(I*J1 

32303 

PNEa* AlMl (I, 3) »P ( I- l,J)**IPl(X, 3) *2 U *1,3) *A3Nl (I, 3) «?( I, 3«1J 

i 

32310 

2NE»l»PN6W*43Pl(I*3)*2(I»3*l)*89H!JtI,3) 

32520 

2(I#J)aPNEwa3'lC3*aH3aP0l.O 

02333 

V  ALa  A5S  UP  C I,  J] -P3L03  7  2  31.31 

02340 

IFCVAL.GT.AMAX)  A4AXa¥Aw, 

32350 

50 

C3NTINUE 

32360 

if(val,ut,e99)  30  ra  ?o 

32370 

60 

caNTiNUE 

a 

32380 

*4ITE(6,651 

32590 

65 

F 3 PM AT  tax, 'EXCEEDS  ITERATIONS' J 

32400 

70 

3  4 AX  *  36« 1 

32410 

1 4AX ■ I 6«  1 

3242Q 

AjaO.O 

32430 

0  3  80  I ■ 1 ,  I M  A  X 

32440 

03  80  3*1, 3 MAX 

32450 

Z*2(I,3)*P(X*l*3)*2(IM,3*l)*P(I,3*tJ 

• 

32460 

80 

AJ2AJ*,25*Z*(X(I*U*X(I)J  *(r(3*l)-r(3)) 

-  • 

.  >470 

3*36*1 

32483 

*3*2, a*3 

3244Q 

A  J  3a  A3 

32300 
32310  : 

a*c»aXX/G/*33 

32520 

A'laaCPHXataAlX) /9C**2*(N**2-1,1 «*a/(»*N»*2*l,) 

32533 

E;C*AM0/O*9C/«2J./3G/3INCNaU/2,/*?n 

32540 

7*9C**3*AH0/£/AXX*(i,*AJ/(Sa*a-l,Ja#a 

• 

32550 

VaPHICi*V/P4I 

32560 

VaPHla»C-AA»4I 

32570 

F 1 lav3PMlG-y3PHI 

32580 

9£TUSN 

32590 

EnO 

32600 

0  0  J  B  L  £  24ECI5I3N  FJNCTI3S  X *  4 0 1 N ( AX , iX , f , T 3U 

32#  1  3 

3 3 jBLE  99ECISIJN  a  j  ,  3  x  ,  F  ,  r  3U 

3  2  620 

0  3  JdUE  0uECIoI34  A  ,  i  ,  Z .  3  .  :  «  EPS  ,  r  A  ,  F  i  ,  FC  ,  r  3(.  1 ,  X  4  ,  P  ,  3  ,  •» ,  3 

32  #  30 

j  a  j  3le  >4£c:s:j4  )a3s 

32640  : 

173 

• 

.  j 

32sso 

32663 

32673 

32683 

32690 

32703 

32713 

32723 

32733 

32793 

32750 

32760 

32773 

32783 

32793 

32900 

32813 

32823 

32933 

32943 

32850 

32963 

32873 

32883 

3299Q 

32903 

32913 

32923 

32933 

32943 

32950 

02960 

02973 

32980 

32990 

33300 

32313 

33320 

33330 

33340 

33330 

33360 

3337Q 

3338Q 

33390 

33100 

33113 

33123 

33133 

J3l«3 

33150 

33U0 

33173 
33133 
33193 
33203 
33213 
33223 
33233 
33240 
3325Q 
33263 
33273 
J  3  2  9  3 
33293 
33303 


3<  COMPUTE  EPS,  T*£  RELATIVE  MaC-IINE  PRECISION 

£PS  «  1,000 

10  EPS  ■  EPS/2.303 

ran  «  l.oDo  ♦  eps 

17  (T3L1  ,GT,  1,303)  G8  18  10 

c 

C-  INITjAUlZATlaN 
C. 

A  a  AX 
3  a  3x 
7*  «  F  C A3 
FB  a  P i 9) 

c 

:  begin  step 

c 

20  C  a  a 
PC  a  Pa 
3  a  a  »  A 
S  «  3 

30  IP  (QA63(FC3  tSE,  3AB3CF35)  ia  TO  40 
A  a  3 
3  a  C 
C  «  A 
PA  a  P9 
F8  a  PC 
PC  a  pA 
C 

3:  caNVE.RGENCE  TEST 
C 

40  T0L1  «  2,OOOa£?SaOAB3{3)  ♦  3,500arat 
X*  «  , 5*  CC  «  8) 

ip  (OAastxn)  , uE,  Tauu  u  u  90 

IP  CP3  ,EQ,  0,0035  Ga  T3  9Q 


•  IS  BISECViaN  NECESSARY 

IP  (0A8SCEJ  ,IT.  T 3L 1  3  33  T0  70 
IP  C0A8SCPAJ  ,l£,  3A3S  (F3) )  Ga  TO  70 

IS  3JA3RATIC  INTERPOLATES  PaSSIdLE 

IP  (A  ,NE,  C]  33  r3  50 

LINEAR  lNT£RPai,ATI3N 

S  a  P9/P A 
*  a  2,000»XN*3 
3  *  1,000  •  S 
G  3  T  3  6  0 

INVERSE  30 AOR  AT  I C  INTERP3LATI JN 

30  3  a  FA/fc 
R  a  p  3/P  C 
S  a  F9/FA 

P  »  Sa(2,3Q0aX4al*!3  .  Rj  .  jg  .  a)*;r  - 

3  •  {3  ■  l  ,300) * (R  •  l,000)»(S  -  l.OJJ) 
ADJUST  SIGNS 

80  IP  CP  ,GT,  0,300)  3  •  -0 
P  *  3AaS(P) 

IS  Interpolation  accept able 


,300)  ) 


174 


4 


X 


L 


[■ 

j: 

t 

, 


33313  '  if  aa,ooa»F)  ,se,  C3,300*x**a  -  oabsu am *aj ) )  oa  fa  fa 

33323  IF  CP  ,G£,  OABS (3. 500»£*3) )  fia  ra  73 

33333  £  »  0 

333*0  3  «  P/U 

33333  33  T2  SO 

33360  C 

03373  C  BISECTION 
33380.  w 

33393  ‘  70  0  ■  X« 

33403  t*0 

33413  C 

33420  C.  COMPLETE  STEP 
33433  C 

33440  80  A  i  8 

33450  FA  •  FB 

33463  IF  (DABS(O)  ,3T,  T3L1)  3*8*0 

33473  IF  (QABSCO)  _,l£,  T JL i J  3*8*  OSIGN(raH»  XI) 

33480  F8  *  F (3) 

33493  IF  C(FB*(FC/0A3S(FC) ) )  ,3T,  0,000)  S3  Ta  23 

33500  02  T3  30 

33510  C 
33520  C  03n£ 

33530  C 

33540  '  90  zesaiN  ■  a 

33553  RETURN 

33560  ENO 


V 


t 
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10 

C 

PROGRAM  NAVSEAL1 

20 

C 

DESIGN  MODIFICATION  PROGRAM  TO  LOWER  THE  TORSIONAL  STIFFNESS 

30 

c 

THE 

EXISTING  SEAL  RING, 

40 

c 

50 

DIMENSION  X<50>, V<30), P<50,  30),  AIM1C50,  30),  AIP1C50,  50) 

£0 

DIMENSION  AJP1C50,  50),  AJMK50,  50),  BRHSC50,  50) 

70 

REAL  K 

30 

5 

WRI TE<6.  10) 

90 

10 

FORMAT (./ ,  5X,  "  INPUT  XI,  VI'  > 

100 

READ <5,  *)  XI,  VI 

110 

WRI TE<6,  20) 

120 

20 

FORMAT </,  5X, ' INPUT  X2,  X3,  V4' ) 

120 

READ <5, *)X2, X3,  V4 

140 

WRITE  <6, 30) 

150 

30 

FORMATC/,  5X,  '  INPUT  D, V6' ) 

160 

READ<  5, *)D, V6 

170 

K*l.  015 

180 

RI-1.  9 

190 

fir=i.  9 

200 

RL*1.  339 

210 

XX3*.  1875 

220 

VV3*.  125 

220 

Tl=.  010 

240 

T2=.  003 

250 

X7*.  36 

260 

X8-.  154 

270 

X9*  051 

230 

X10=. 1435-X1 

290 

X4-X2-X3 

200 

X3*XX3+2.  *T  1 

310 

31 

X6=RI-T1-RL 

320 

V2*.  38-V1 

330 

V3*.  22-V4 

340 

VS*.  1 

350 

V7-.  09-< .  1-V6) 

360 

V8=.  08 

370 

V9=.  35 

330 

V10*.  08 

330 

P I *3.  141592634 

400 

AMU*.  278 

410 

PU*S00. 

420 

PL*.  5602 

420 

N=9 

440 

E-31.  E+06 

450 

G=E/2.  /Cl.  +AMU) 

460 

PG=750.  0 

470 

100 

A1=X3*V3 

480 

A2=< V3*V4 j*X4 

490 

A3*X1*(V2*V3+V4-V5) 

500 

A4=X10*(VH-V2*V3*V4-V5> 

510 

A5=X9*(  Vl-*-V2-*-V3+V4-V5-V10) 

520 

A6=<X7+X3-X6)*<V7+V8-V6) 

330 

A7*tX6-X7>*< V7+VS) 

540 

A3=X7*V7 

550 

I F < Al.  LT.  0.  0)  GO  TO  26 

560 

IFCA2.  LT.  0.  9)  GO  TO  26 
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570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710  C 
720 
730 
740 
750 
760 
770 
780 
790 
300  C 
310 
320 
830 
340  C 
350 
860 
870 
830  C 
890 
900 
910 
920 
930 
940 
950 
960 
970  C 
980 
990 
1000 
1010 
1O20 
1030 
1040 
1050 
1060  C 
1070 
1080 
1090 
1100  C 
1110  C 
1120 
1130 
1140 
1150  C 
1160 
1170 
1120 
1190 
1200 
1210 
1220 


IF<A3.  LT.  0.  0)  GO  TO  36 
IF<A4.  LT.  0.  0)  GO  TO  36 
IF<A5.  LT.  0.  0>  GO  TO  36 
I F<A6.  LT.  8.  0>  GO  TO  36 
IF<A7.  LT.  0.  0)  GO  TO  36 
IF<A8.  LT.  0.  0)  GO  TO  36 
FXAl=*X3/2. 

FXA2»X3+X4,'2. 

FXA3*X2+Xl/2. 

FXA4=*Xl+X2+X10/2. 

FXA5*Xl+X2+X104-X9/2. 

FXfl6=»Xl+X2-t-X10+X9+CX3-t-X4+X5-X2-Xl-X10-X9)/2. 

FXA7=X3*X4+X5-KX6-X?>/2. 

FXA8=X3't-X4+X5+X6-X7/2. 

FYRl=YH-Y2+Y3/2. 

FVA2*Vl-t-V2'KV3-t-V4>/2. 

FYA3*Yl+<Y2+V3+V4-V5>/2. 

FYA4=CYl*Y2+Y3+-V4-V5)/2. 

FYA5*Y10>< VI +Y2+V3+V4-Y5-Y10) /2. 

FYA6«Y10+Y9-KY8*V7-Y6>/2. 

FYA7=Y10+V9+CV8+Y7>/2. 

FYA8*Y10-*-V9-fY8+Y7/2. 

X6AR»Al*FXAl+A2*FXA2-t-A3*FXA3«-A4*FXA4 
XBAR»XSAR-t-A5*FXA5+R6*FXfl6+A7*FXA7-t-A8*FXfiS 
XBAR^BAR/C  Al-t-A2*A3+'R4  +  R5'*-R6'<-A?+R8> 

Y6AR»Al*FYAH-A2*FVR2-»-A3»FYR2  +  A4*FYA4 
YBRR=VBAR+A5*FVft5+A6*FYB6+ft7*FVA7+A6*FVA3 
VBAR»VeARX<flH-A2t-R2-t-A4  +  B5-i-fl6«-R7*B8> 

RIX*1.  712.  *X3*V3**3+fll*(FVAl-YeAR>**2 
flIX»flIX*i.  /12.  *X4*< V' 3*  Y4 > **2  +R2* C F VB2- VBRR>  **2 
BIX-flIX+1.  /12.  *X1*<Y2*Y2«-Y4-Y5>**3+A3*<FYA2-YBRR>**2 
RIX*flIX+i.  X12.  *X10*C  Vl+Y2+Y2+,Y4-Y5>**2'*-H4*<FyR4-YBflR>  *»2 
BIX-RIX+1.  X12.  *XS*<Yl-t-y2+Y3*Y4-Y5-V10>**3+A5*<FYA5-YBRR>** 
BIX=»BIX+1.  X12.  *(.X8-<X6-X?))*(Y?+Y$-Y6>**J+?I6*(FYR6-VSFIR>** 
BIX»BIX+1.  X12.  *<  X8>  *( V7+Y8) **2+B7*<  FYA7-YBRR>  **2 
flIX=RIX+l.  X12.  *X7*(Y7 >**2*R8*<FYB8-VBBR>**2 

BIV=1.  /12.  *X3**2*Y2*Rl*<FXfll-XBAR>**2 
AIV«=BI  Y+l.  /12.  *X4**3*( Y3+V4>  +  R2*(FXB2-X£RR>  **2 
BI  Y=BI  Y+l.  X12.  *Xl»+3*t  V2+Y3+V4-Y5>-»-A2*(FXfl2-XBAR>**2 
A I  Y»A  I  Y  +  l.  /12.  t‘X10**3*<:  Yl-,rV2+Y3'*-Y4-Y5>  +R4*<  FXR4-XBAR >  **2 
BIY=BIV*1.  X12.  *X9**3*<Yl+Y2+V3*V4-V5-Y10>-t-A5*<FXA5-XBAR>*» 
BIY=BIY+1.  /12.  *  <X8-<X6-X7>  >**3*(  V7*Y8-Y6> *A6*C FXA6-XBAR > ** 
AIV=BIY*1.  X12.  *XS*,*3*<Y7+Y8>'*-B7*(FXA7-XBAR.''**2 
AIV=*AIY+1.  712.  *X?*»<3*V7 *AS*<  FXA8-X8BR>  **2 

R00*RL+X6+X5+X4-fX3 

RC=ROO-XBBR 

RP I  =2.  484 


BJR=AIX+BIY 

£IX=E*AIX 

AJO=(.  025*(  Al't-R2-**A3,i-A4't-A5+A6+R7+A8>  **4  JXAJR 
TORSIONAL  CONSTANT  CALCULATION  USING  FINITE  DIFFERENCE 
ERR=0.  0000001 
AA1=X7 
AA2=X6 
AA2  =  X7+X'S 
hA4*X7-‘-X3+-X9 
AA5  =  X7 +XS+X9+X13 
mA63X5*X6 
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1230 

R87=X4'*-X5+X6 

1240 

Afl8“X3+X4+X5+X6 

1250 

BBi*Yie 

1260 

SB2-Y1 

1270 

8B3=Y1+Y2 

12S0 

BB4-Y9+Y10 

1290 

665-Y1+Y2+V3+V4-Y5 

1300 

686=Y8+Y9+Y10 

1310 

BB7=Y1+Y2+Y3 

1320 

BBS*V7+Y8-t-Y9+Y10 

1330 

BB9=V1+Y2+Y3+V4 

1340 

11*4 

1350 

12=5 

1360 

13=13 

1370 

14=17 

1330 

15=19 

1390 

16=20 

1400 

17=21 

1410 

18=22 

1420 

19=23 

1430 

Ji=17 

1440 

J2=23 

1450 

J3=24 

1460 

J4=26 

1470 

J5*29 

1480 

J6=30 

1490 

J7=36 

1500 

J8=37 

1505  C 

MESH 

GENER8TI0N  ROUTINE 

1510 

XX  =  0.  0 

1520 

Xil)=0.  0 

1530 

DX=6Bi/( Ii-1) 

1540 

00  60  1  =  2,  19 

1550 

IFCI.GT.  ID  0X=CBB2-BB1)XCI2-I1> 

1560 

I F < I .  GT.  12)  0X*<BB3-BB2)X< 13-12) 

1570 

IFC I.  GT.  13)  DX*CBB4-BB3> X ( 1 4- 1 3 > 

1580 

IFCI  GT  14)  0X=(B65-BB4)X(I5-I4> 

1590 

IFCI  GT.  15)  OX=t  BB6-BB5>/( 16-15) 

1600 

IFC I.  GT  16)  DX=C6B7-6B6)XC 17-16) 

1610 

IFC I. GT.  17)  DX=CBB8-B67>/< 18-17) 

1620 

IFCI.GT.  18)  0X=CB89-BB8)XC 19-18) 

1630 

XX-XX+DX 

1640 

XCI)=XX 

1650 

60 

CONTINUE 

1660 

YV=0.  0 

1670 

YC 1 )  =0.  0 

1680 

0V=RR1/C  Jl-i) 

1690 

00  70  J=2,  J8 

1700 

IFCJ.  GT.  Jl)  0V=CRR2-RR1 >XCJ2-J1) 

1710 

IFCJ.  GT.  J2)  0Y=CRR3-RR2)XC J3-J2) 

1720 

IFCJ.  GT.  J3)  DY= C 884-883 >XCJ4-J3> 

1730 

IFCJ  GT.  J4>  OY=Cflfl5-0ft4)XCJ5-J4> 

1740 

IFCJ  GT.  J5>  0Y=CRR6-RR5)/CJ6-J5> 

1750 

IFCJ.  GT.  J6>  DY=CRR7-Rfi6)/C J7-J6) 

1760 

IFCJ.  GT.  J7>  0Y=  (.888-887) /CJS-J7) 

1770 

YY=YV+OY 

1730 

V  C  J ) =VY 

1790 

70 

CONTINUE 

1300 

IM8X=I9-i 

1810 

JM8X=JS-1 

1320 

00  30  1=2,  IM8X 

1830 

00  30  J=2,  JM8X 

1346 

81  J=-2.  /c:<c  I-1)-XC  I  >  )/(XC  I+D-XC  t  >  >-2.  /C  •'<  J-i)-Y<  J  '■  V 

1350 

* 

CYC  J+D-VC  J)  ) 

i860 

8 1  Ml  Cl,  J>=2.  /CXC  I-D-XC  I  )  )/CXC  I-D-XC  I-l)  )/'8I  J 

1870 

8 1  Pic  I,  J)  =-£.  /cxc  i  i  >  >/cxc :  -i>-xc :  -i  >  .vftu 
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"  '  ■ — - . — 

■  • 

k-  • 

1880 

AJM1C I ,  J>*2.  VC  J-D-VC  J)  >/<Y< J-l)-Y< J*l) )/AI  J 

r 

1890 

AJP1CI,  J)*-2.  /<V< J+l)-Y< J> )/CY< J-D-VC J+l) >/AI J 

X 

1900 

BRHSC I, J)=2.  /AIJ 

1910 

80 

CONTINUE 

1920 

DO  90  I»l,  19 

1930 

DO  90  J=i,  J8 

1940 

PCI,  J)=*0.  0000000000001 

1950 

90 

CONTINUE 

► 

I960 

OMEG-1.  5 

a 

1970 

0M0=1.  -OMEG 

1980 

DO  95  IT*1,  800 

1990 

AMAX=0.  0 

2000 

DO  94  J=2,  JMAX 

2010 

IMAX=I8-1 

2020 

IMIN= I 6+1 

2030 

IF< J.  GT.  Jl)  I M I N= I 4+1 

2040 

IFCJ.  GT.  J2)  I MAX= I 5-1 

■l 

2050 

IF<  J.  GT.  J3)  IMIN-Il't-1 

.1 

2060 

IFCJ.  GT.  J4>  I MI N  =  2 

20?  i 

IFCJ.  GT.  J5)  IMIN=I2*1 

2080 

IFCJ.  GT.  J6)  IMIN=I3-t-l 

2090 

IFCJ.  GT.  J6>  I MftX=  I 9-1 

2100 

IFCJ.  GT.  J?)  IMfi»I7-l 

2110 

DO  94  I*IMIN, IMAX 

2120 

POLD-PC I,  J) 

*  m 

2130 

PNEW=AI M1C  I ,  J)*PCI-1,  J)+AIP1CI,  J>*PCI+1,  J) 

W 

2140 

PNEW=PN£W*RJM1C  I,  J)*P(I,  J-D+-AJP1C  I,  J)*PCI,  J+D+BRHSC  I,  J) 

2150 

PCI,  J)=PNEW*OMEG-t-OMO*POLD 

2160 

VAL=ABSC  CPC  I,  J)-POLD)/POLD> 

21?0 

IFCVAL.  GT.  RMAX)  AMAX=VAL 

2180 

94 

CONTINUE 

2190 

IFCVAL.  LT.  ERR)  GO  TO  110 

2200 

95 

CONTINUE 

rr 

2210 

WRITEC6,  96) 

l 

2220 

96 

FORMAT C2X,  'EXCEEDS  ITERATIONS' > 

2230 

110 

JMAX*J8-1 

2240 

IMAX=I8-1 

2250 

A J=0.  0 

2255 

C 

INTEGRATION  ROUTINE 

2260 

DO  115  1=1,  IMAX 

2270 

DO  115  J=i,  JMAX 

mm 

2280 

Z=P<  I,  J)*PC  I+l,  J)+P(  t+1,  J4-D+PC  I,  J<-1) 

S 

2290 

AJ=AJ+.  25*2*<XCI+1)-X(I  >  )*CVC  J-t-l)-VCJ)  > 

2300 

115 

CONTINUE 

2310 

J=J3*1 

2320 

2330 

AJO=AJ 

2340 

116 

CONTINUE 

2350 

A=E,t‘AIX/G/AJO 

r  m 

2360 

ECC=PL+-VBAR 

2365 

c 

WARPING  FACTOR  CALCULATION  C  RJMEGA) . 

2370 

XT0T*X2-kX4«-X5+X6 

2380 

TlAV=V7*CX7/XT0T)+<  V7+-VS)  *C  C.X6-X?)/XT0T)  +  <  VS+YT-VS) * ( X5/X70T  > 

2398 

* 

+  C  V2+V4)  »C  X4/XTGT)  +-V2«C  X2/XT0T ) 

2400 

YTOT=V9+V10 

2410 

T2AV=X10*CV10/VTOT)-<-C.X9*X10>*C  C  V1-V10)  XVTOT) -t-C  X1  +  X9*X10 ) 

2420 

* 

*C C V2+V2+V4>-CY?*V8> )/VTOT 

2430 

H=V10+V9+TlAV/2. 

w 

2440 

BDIM=XTOT 

2450 

AJMEGA  =  TlAV'**3*8DIM**3/144<-TiAV«'*2*H**3/2S 

2460 

B-AIX*RC**2/AJMEGA 

2470 

FflCl*N**2/A+N*>»,4/6  +  N**2 

2480 

FAC2=N**2-*-l.  /A+N**2/B 

2490 

FAC2  =  1.  *1.  /A-t-N** i/S 

2500 

F AC4=N**2/A+N**4/g- i. 

'  • 

2510 

AMO=£CC*0*RPI  *PG*SIMCN*C*/C2.  *>PP r  >  >/RC 
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2520  PHI=AM0*R'C**27E/AIX*<1.  /<  <N**27A-t-N**47B*l.  >-<l.  +  1.  7A+N**27B> 

2520  *  7<N**2-t-l.  7A-t-N**278)*<N««*27A-i-N**47B+N**2)  >  > 

2540  V=AM0*RC**27E7AIX*(1.  7<FAC1-FAC27FAC2*FAC4> ) 

2550  RSEAL=RC-V7PHI 

2555  C  PRESSURE  MOMENT  CALCULATION 
2560  Fl=2.  *PI*<RL+X7+XS+X9)*V10 

2570  F2=PI*<R0Q**2-R8**2) 

2580  F3=2.  *PI*<RL+X7'i-X8't-X9+X10)*Vl 

2590  F4=2.  *PI*<RL't-X7+XS+X9't-X10'i-Xi)*V2 

2600  F5=2.  *PI*<ROO)*Y3 

2610  F6=PI*<R00**2-<R00-X2)**2) 

2620  F7*2.  *PI*<R00-X3)*V4 

2630  F8=PI*<  <ROO-X2)**2-<ROO-<X44X24T1>  >**2) 

2640  F9=2.  *PI*<.'RH-XX2)*VY3 

2650  F10=F2-F6-FS 

2660  AM1*-F1*<YBAR-V1072.  ) 

2670  R2  =  2.  73.  *<R00**3-RB**3)7< R00**2-RB**2> 

2680  AM2=-F2*<R2-RC) 

2690  AM2=F3'*‘<YBAR-Y172.  ) 

2700  AM4=F4*<  VBAR-c'  V1+V2/2.  >> 

2710  AM5=-F5*<Y1«-V2+V372.  -VBAR> 

2720  R6=2.  72.  *<  R00**3-<.  R00-X2>**2)  7  ( R00**2-<  R00-X2)  **2) 

2720  AM6=F6*<R6-RC> 

2740  AM7*-F7*<Y1+Y24Y3*V472.  -VBAR) 

2750  R8=2.  73.  *< CR00-X3) **2-<R00-X2-X4-Tl> **2> 7< (R00-X2) **2-( ROO-X3- 

2760  *  X4-Ti)**2) 

2770  AMS=F8*<RS-RC) 

2780  AM9=-F9*<Y2+Vl+Y3-t-Y4+VV372.  -VBAR) 

2790  R10=2.  72.  •XX2+-RI 

2800  AM10=F10*<R-"'-RC) 

2810  TM=<AMH-AM2+AM2+AM4*AM5+AM6+AM?+AM8+AM9''-AM10)7C2.  *PI*RC> 

2820  URITE<6,  «)  B 

2820  R0TATE=TM*RC**27E7AIX*PU 

2840  GO  TO  39 

2350  26  URITEC6. 27>A1,  A2, A3- A4, AS, A6 

2960  27  F0RMAT<5X, 'NEGATIVE  AREAS  ARE  PRODUCED  BV  THOSE  DIMENSIONS'*  7,  5X 

2370  *  *  '  Al='  *  F6  4,  2X,  '  A2*'  *  F6  4,  £x,  '  A3*'  ,  F6.  4,  2X,  '  A4='  ,  F6  4,  2X,  '  A5='  , 

2880  *  F6  4,  2X,  '  A6=' ,  F6.  4,  777) 

2890  GO  TO  5 

2900  29  WRITE<6,  40)  XBAR,  VBAR*  AIX,  AIY,  AJR,  AJO,  A*  AMO, V*  ECC, RC,  EIX,  ROTATE 

2910  40  FORMAT ( 7,  5X,  '  XBAR*' ,  F7.  4,  7,  5X,  '  VBAR*'  ,  F7.  4,  7,  5X,  ' IX*' ,  E12.  5,  7,  5X, 

2920  *  '  I  V*' ,  E12.  5,  7,  5X,  '  JR*' ,  E12.  5.  7,  3X<  '  JTHETA*' ,  E12.  5,  7,  5X,  '  A=' ,  E12.  5 

2920  *  ,  7-  5X>  '  MTHETA*'  ,  E12.  5,  7,  5X,  '  V='  ,  E12.  5,  7-  5X,  '  ECC*'  ,  E12.  5,  7,  5X, 

2940  *  '  RC='  ,  F7.  4,  7,  5X,  '  EIX*'  ,  E12.  5,  7,  5X,  '  ROTATION®'  ,  E12.  5> 

2950  URI TE<  6, 45)  RSEAL,  AIXY,  XI, X2, X2, X4,  .  5, X6, X7,  XS, X9, X10, PHI 

2960  URI  TE(6,  50)  VI,  Y2,  Y3,  Y4,  Y5,  V6,  V7,  VS,  V9,  V10,  ROO,  RPI,  PL,  RI,  RL 

2970  45  FORMAT (5X,  '  RSEAL='  ,  E12.  5,  7,  5X,  '  IXV*'  ,  £12.  5, 

2980  *  7,  5X,  '  XI*'  ,  F6.  4,  2X,  '  X2='  ,  F6  4,  2X,  '  X2=' ,  F6.  4,  2X,  '  X4='  ,  F6  4,  2X, 

2990  *  '  X5*' ,  F6.  4,  2X,  '  X6='  ,  F6.  4,  2X,  '  X7*'  ,  F6.  4,  7,  5X.  '  XS='  ,  F6.  4,  2X,  '  X9*'  , 

3000  *  F6  4,  2X.  '  X18*' ,  F6.  4,  7,  5X,  '  PHI*' ,  E12.  5) 

3010  50  FORMAT <  SX,  '  VI*'  -  F6.  4,  2X,  '  V2='  ,  F6  4,  2X,  '  Y2*'  -  F6.  4,  2X,  '  Y4*'  ,  F6.  4 

3020  *  ,  2X,  '  Y5='  ,  F6.  4,  2X,  '  V6*'  ,  F6.  4,  2X,  '  V7='  ,  F6.  4,  7,  5X,  '  V8*'  ,  F6  4 

3020  *  ,  2X,  '  V?*'  ,  F6.  4,  2X,  '  V10*'  ,  F6  4,  7,  5X,  '  ROO*'  ,  F6  4,  2X, 

2040  «  '  RPI='  ,  F6.  4,  2X»  '  PL*'  ,  F 6.  4,  2X,  '  RI='  ,  FS.  4,  2.';,  '  RL  =  '  ,  F5  4,  777) 

2050  GO  TO  5 

3060  END 
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**** 


Program  Listing 


*  *  *  * 


The  fluid  used  is  the  nitrogen  gas. 

Alternating  Direction  Implicit  technique  is  used  to  analyze 
a  moving  wave  gas  seal. 

A  set  of  simultaneous  equations  at  each  time  step  are  solved 
by  a  special  Gaussian  Elimination  Subroutine,  in  which  a 
tridiagonal  coefficient  matrix  is  rearranged  in  a  form  of 
three  columns  matrix.  The  advantage  is  to  save  a  lot  of 
computer  space,  i.e.  n  x  n  storages  decreases  to 
n  x  3  storages. 

Implicit  real*8  is  used  to  convert  all  real  single  precision 
numbers  to  double  precision  numbers. 

The  seal  is  divided  into , three  identical  parts  (3  waves), 
each  part  is  subdivided  into  19x20  nodes,  where  19 
nodes  are  in  circumferential  directions  and  20  nodes 
are  in  radial  directions. 

The  values  used  in  this  program  have  the  following  units: 


C  B 
C  c 
C  DT 
C  ETA 
C  EM 
C  GM 
C  G 
C  HA 
C  PC 
C  PM 

C  PI 
C  PO 
C  PSP 
C  PSBAR 
C  PHI 
C  RB 
C  RC 
C  RI 
C  RO 
C  RE 
C  RHO 
C  RPM 
C  TEMP 
C  UCON 


Balanced  ratio 

One-half  maximum  roughness  height 
Time  step 

Viscosity  of  nitrogen 
Young's  Modulus 
Shear  Modulus 
Gravity 

Ring  centroid  amplitude 
Cavity  pressure 
Pressure,  at  asperity  contact 
(equals  compressive  strength) 

Seal  inside  pressure 
Seal  outside  pressure 
Spring  pressure  on  face 
Shear  strength  of  pressure 
Tilt 

Seal  balance  radius 
Seal  centroid  radius 
Inside  radius 
Outside  radius 
Friction  radius 
Density 

Shaft  angular  speed 
temperature 
Universal  constant 


:  in 
:  sec 

:  lbf- sec/in/in 
:  psi 
:  psi 

: ft/sec/sec 
:  in 
:  psi 
:  psi 

:  psi 
:  psi 
:  psi 
:  psi 
:  in/in 
:  in 
:  in 
:  in 
:  in 
:  in 
:  in 
:  rpm 
:  °R 

:  ft- lbf/lbm-R0 
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The  main  program 


IMPLICIT  REAL *8  (A-H,0-Z) 

DIMENSION  NO (18) 

COMMON  /BLOCKA/C , DA , DEM , DR , ETA , OMEG , P I , PO , RO , WE , IMAX , JMAX 

COMMON  /BLOCKB/A^DT.HA^H^R^WV^NIT 

COMMON  /BLOCKC/H( 18, 20) , P( 18, 20 ) , Q( 18 , 20 ) , RR( 18 , 20 ) 

COMMON  /BLOCKD/AMAXER, RI , NN, IM1 , JM1 

COMMON  /BLOCKE/PFAC , PHLS , PM, TOTLD , WMECH 

COMMON  /BLOCKE/FCONV , WSTAR , ROOC , NWAVE 

IMAX=19 

JMAX=20 

JM1=JMAX- 1 

IM1=IMAX-1 

NWAVE =3 

NN=NWAVE 

RC=1 . 9361DO 

HA=5 . 4D-6 

PIE=3 . 141592654D0 

EM=3.D+6 

GM=EM/2 . 4D0 

AN=2.D0*PIE/NN 

DA=AN/IM1 

AMAXER=6  .'ODO 

B=0 . 75D0 

C=20.D-6 

DT=0 . 5D-6 

PSP=28 . 5DO 

PC=0 .  ODO  '• 

G=32. 174D0 ' 

RI=1 . 900D0 
RO=2 . 0875D0  ' 

WLIMIT=3 . ODO 
DR=(R0-RI)/JM1 
PI=14 . 7D0 
PO=514 . 7D0 

RB=DSQRT(RO**2-B*(RO**2-RI**2) ) 

PHI=500 . D-6 

RHO= . 0023D0 

RPM=1800 . ODO 

OMEG=2 . D0*PIE*RPM/60 . DO 

SSTIFF=23 . DO* . 833D0 

DISP= . 29125D0 

SPNO=12 .DO 

SPLOA=SSTIFF*DISP*S?NO 

ETA=2 . 5236D-9 

TEM?=598. 67D0 

UCON=53 . 3D0*12 .DO 

DEM=-P IE/12 .DO/ETA/UCON/TEMP 
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FPLOA=PIE*(PO*(RO**2-RB**2)+PI*(RB**2-RI**2) ) 

ROOC=RO/C 
PM=38000.D0 
PSBAR=0 . 1D0*PM 

RF=2.DO*(RO**3-RI**3)/3.DO/(RO**2-RI**2) 

WSTAR=SPLOA+FPLOA 
WRITE( 6,1) 

1  FORMAT* 1H1,/,15X, '************************************** 

*********************< ,//21X, 'TYPICAL  INPUT  PARAMETERS 
*AND  OUTPUT  RESULTS",//, 15X, ' _ 

* _ '  ,/) 

WRITE (6,2)  WSTAR, HA, PHI , AMAXER, WLIMIT 

2  FORMAT( IX, //15X, 'REQUIRED  LOAD  SUPPORT  =' , D12 . 5 , 4X, 

*' LBS ' ,//15X, 'WAVINESS  AMPLITUDE  HA  = ' , D12 . 5 , 4X, ' INCH' 
*,//15X, 'TILT  PHI  =' ,D12. 5, 17X, ' IN/IN' ,//15X, 'PRESSURE 
*CONV.  CRIT.  AMAXER= ' , F6 . 3 , 5X, ' PSI ' , //15X, 

* ' LOAD  CONV.  CRITER.  WLIMIT  = ' , F4 . 2 , 3X, ' LBS ' , //) 

WRITE ( 6,3)  B, C , RPM, DT, PI , PO 

3  FORMAT (15X, ' BAL  RATIO  =' , F6 . 3 , //15X, ' HALF  MAX  ROUGH  HEIGHT 
*C  =' ,D12.5,2X, ' INCH'  //15X, 'RPM  = ' , F8 . 1 , //15X, ' HALF 
♦TIME  STEP  DT  =', D12 . 5 , 8X, ' SECONDS ' ,//15X, ' PI  =' 

*,F8.2,7X, 'PSI' ,//, 15X, 'PO  =' , F9 .2, 6X, ' PSI ' ,//) 

RWW=RC -HA/PHI 

RNO= ( RWW-RI ) / ( RO-RI ) 

WRITE ( 6,4)  RWW , RNO , RB 

4  FORMAT(15X, 'RADIUS  OF  CONTACT  POINT  =' ,D12 . 5, 2X, ' INCH' 
*,//I5X, 'DIMENSIONLESS  OF  CONTACT  RADIUS  = ' , F5 . 2 , //15X, 
♦'BALANCED  RADIUS  =  ' , F12 . 5 , 9X, ' INCH ' , //) 

HOB=29 . 5D-6 
CALL  WLOAD ( HOB ) 

WOB=WE 
NITB=NIT 
HOA=33 . OD-6 
CALL  WLOAD (HOA) 

WOA=WE 

NITA=NIT 

WRITE (6,5)  HOA , HOB , WOA , WOB 

5  FORMAT (15X, 'HOA  =' , D12 . 5 , 2X, ' INCH' , 16X, ' HOB  =',D12.5 
*,2X, ' INCH' ,//, 15X, 'WOA  =' , D12 . 5 , 2X, ' L3S ' , 17X, ' WOB  =' 

* , D12 . 5 , 2x, 'LBS' ,//) 

100  HNEW=HOA+ ( HOB-HOA ) * ( WOA-WSTAR ) / ( WOA- WOB ) 

CALL  WLOAD (HNEW) 

WNEW=WE 

NITN=NIT 

IF ( DABS (WNEW- WSTAR) .LT. WLIMIT)  GO  TO  200 

KOB=HNEW 

WOB=WNEW 

GO  TO  100 

200  WRITS (5, 6)  HNEW, WNEW 

6  FORMAT (15X, 'MINIMUM  FILM  THICKNESS  H0= ' , D12 . 5 , IX, ' INCH ' 
*,//15X, ' FLUID  LOAD  SUPPORT  WNEW  = ' , D12 . 5 , 2X, ’ L3S ' , //) 

WRITE(S,7)  NITA, NIT3, NITN 

7  FORMAT (15X, 'NITA  =' , I4,//15X, ' NIT3  =',I4,//15X 
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*, 'NITN  =  ' , 14,//) 

HST=HNEW+0 , 5D-6 
CALL  WLOAD(HST) 

WOST=WE 

STI FF= ( WOST-WNEW ) / ( HST-HNEW ) 

WRITE(6/8)  STIFF 

8  FORMAT(15X, 'STIFFNESS  OF  FLUID  =’ ,D12 . 5, 7X, ’ LBS/IN '  //) 
FMI=O.DO 

DO  600  1=1,11*1 
J=2 

JMM=J-I 

JP1=J+1 

JP2=J+2 

JP3=J+3 

JP4=J+4 

600  FMI=FMI+DEM*RI*H(I,1)**3*( (Q(I, J)-Q(I, JMM))-(Q(I, JP1) 
*-2.D0*Q(I, J)+Q(I, JMM) )/2.D0+(Q(I, JP2)-3.D0*Q(I, JP1) 
*+3.D0*Q( I , J)-Q( I, JMM) )/3.D0-(Q( I, JP3) 

*-4.D0*Q( I, JP2)+6.D0*Q( I, JP1)-4.D0*Q(I, J)+Q(I, JMM) )/4.D0 
*+(Q(I, JP4)-5.D0*Q(I, JP3)+10.D0*Q(I, JP2 ) -10 .D0*Q( I , JP1) 
*+5.D0*Q(I, J)-Q(I, JMM) )/5.D0)/DR 
FMI=FMI *NN/G/RHO 
WRITE( 6, 10)  FMI 
FMO=0 . DO 
DO  700  1=1, IM1 
J=JMAX 
JMM=J-1 
JM2  = J  -  2 
JM3=J-3 
JM4=J-4 
JM5=J-5 

700  FMO=FMO+DEM*RO* (H( I , JMAX) **3)*((Q(I, JMM) -Q( I,J))-(Q(I, JM2 ) 
*-2 -D0*Q( I , JMM) +Q( I , J) )/2.D0+(Q(I, JM3)-3.DO*Q(I, JM2) 

*+3 .d) *Q( I , JMM) -Q( I , J) )/3 . DO- (Q( I , JM4) -4 .DO 

**Q( I , JM3 ) +6 . D0*Q( I , JM2 ) -4 . D0*Q( I , JMM) +Q( I , J ) ) /4 .DO 

*+(Q(I, JM5)-5.D0*Q(I, JM4)+10.D0*Q(I, JM3 ) -10 . D0*Q( I , JM2 ) + 

*5 .DQ*Q( I , JMM) -Q( I , J) )/5 . DO) /DR 
FMO=- FMO*NN/G/RHO 
FMN=FMI-FMO 
WRITE (6,9)  FMO, FMN 

9  FORMAT(15X, 'OUTER  LEAKAGE  QOUT  =’ ,D12 . 4, 7X, ' SCFM'  //  • 

* 15X, 1  MASS  FL  ERR0R=' ,D12 . 4, 13X, ' SCFM' ,//) 

10  FORMAT ( //15X, ' INNER  LEAKAGE  QIN  = ' , D12 . 4, 8X, ' SCFM' , //) 

CALL  FRICT(?SBAR,RF,AMU,FFORC/FFA,TORQ) 

WRITE ( 6 , 11)  AMU , FFORC , PHLS , WMECH, FFA, TORQ 

11  FORMAT ( 15X, ' COEF  OF  FRICT= ' , F12 . 8 , //15X, ' FRICT  FORCE=' , 
*D12.5,26X, 'LBS'//,15X, 'PERCENT  OF  FLUID  LOAD  =',F10  S//' 
*15X, 'TOTAL  MECH  CONTACT  SHEAR  FORCE  = ' , F10 . 5 , 9X, ' LBS '  // 
*1SX, 'TOT  FRICTION  DUE  TO  ASPERITY  CONTACT  =' , F13 . 6 LBS ' // 
*15X, ' TOTAL  TORQUE  = ' , F13 . 6 , 24X, ' IN-L3 ' , ///) 

DO  300  1=1, IM1 
NO ( I ) =1 
800  CONTINUE 
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WRITE(6, 12) 

12  F0RMAT(1H1,48X, ' ******  FLUID  FILM  THICKNESS  *******  /SOX 

**(  INCH  )',////)  ‘ 

WRITE(6, 13 )  ((N0(IR)/H(IR/JC),IR=1/IM1)/ JC=1, JMAX) 

13  FORMAT (9 ( IX, 12 , IX, DIO .4) ,  /  ) 

WRITE(6, 14) 

14  FORMAT( 1H1, 48X, ' ******  PRESSURE  DISTRIBUTION  ******* 

* ,/62X, ' (  PSI  V  ,////) 

WRITE ( 6,13)  ( (NO( IR) , P ( IR, JC) , IR=1 , IM1 ) , JC=1 , JMAX) 

WRITE( 6, 15 )  .  * 

15  FORMAT (//// ,  40X, ************************************* 

********',//) 

STOP 

END 
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c 

c 

C  This  subroutine  is  used  to  find  the  pressure  distribution 
C  all  over  the  seal  surface. 

C  A  transformation  of  tridiagonal  coefficient  into  three 
C  column  matrix  is  also  shown. 

C 

C 

SUBROUTINE  WLOAD(HMA) 

IMPLICIT  REAL *8  (A-H,0-Z) 

DIMENSION  NO( 18) , FOLD( 18,19) 

DIMENSION  RIMH( 18, 19),RIPH (18,19), RJMH( 18,19), RJPH( 18,19) 
DIMENSION  HIMH( 18,19), HIPH( 18,19), HJMH( 18, 19),HJPH(18, 19) 
DIMENSION  A2 ( 18 , 19 ) , A2 JM1 (18,19), A2 JP1 (18,19) 

DIMENSION  A1IP1(18,19),A1(18,19) , A1IM1( 18,19) 

DIMENSION  RJM1 (18,19), RJP1 (18,19)  , QA( 18, 19 ) , PA( 18, 19 ) 
DIMENSION  QIM1(18, 19) ,QIP1(18, 19) ,QJM1(18, 19) ,QJP1(18, 19) 
DIMENSION  2(18,5), R(18) 

DIMENSION  Z2( 18, 3 ) , ZR( 18) 

COMMON  /BLOCKA/C , DA , DEM , DR , ETA , OMEG , P I , PO , RO , W , IMAX, JMAX 

COMMON  /BLOCKB/AN, DT, HA, PHI , RC, WW, NIT 

COMMON  /BLOCKC/H( 18 , 20 ) , P ( 18 , 20 ) , Q( 18 , 20 ) , RR( 18,20) 

COMMON  /BLOCKD/AMAXER, RI , NN, IM1 , JM1 

COMMON  /BLOCKE/PFAC , PHLS , PM , TOTLD , WMECH 

IM2=IMAX-2 

JM2=JMAX-2 

IM=( IMAX+1 ) /2 

IMM=IM-1 

IMP=IM+1 

CON=6 . DO*ETA*OMEG 

VALUE=12 . DO*ETA/DT/DR 

RD=RO-RI 

PD= ( PO-P I ) /RD 

DO  100  J=2,JM1 

DO  100  1=1, IM1 

RJM1 ( I , J ) =RI + ( J-2 ) *DR 

RR( I , J)=RI+ ( J-l ) *DR 

RIMH( I , J)=RR( I , J) 

RIPH( I , J) =RR( I , J) 

RJP1( I , J)=RI+J*DR 
RJMH( I , J) =RI + ( J-l . 5D0 ) *DR 
RJPH( I , J)=RI+( J-0 . 5D0) *DR 
100  CONTINUE 

DO  105  J=1 , JMAX 
DO  105  1=1, IM1 
RR( I , 1 )=RI 
RR ( I , JMAX ) =RO 
105  CONTINUE 

DO  107  J=1 , JMAX 
DO  107  1=1, IM1 
CALL  WEAR ( RR ( I , J ) ) 

H( I , J)=HMA+WW+(HA+  (RR ( I , J ) -RC ) *PHI ) *DCOS ( NN* ( I -1 ) *DA) 

107  CONTINUE 
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DO  109  J=2 , JM1 
DO  109  1=1, IM1 
CALL  WEAR  (RR(I,J)) 

HIMH ( I , J ) =HMA+WW+ ( HA+ { RR ( I , J ) -RC ) * PHI ) *DCOS ( NN* ( I - 1 . 5D0 ) *DA ) 
HIPH( I , J)=HMA+WW+(HA+ (RR( I , J)-RC) *PHI ) *DCOS(NN* ( I- . 5D0) *DA) 
CALL  WEAR(RJMH( I , J) ) 

H JMH ( I , J) =HMA+WW+ ( HA+ ( RJMH ( I , J ) -RC ) *PHI ) *DCOS ( NN* ( I - 1 ) *DA ) 
CALL  WEAR(RJPH( I , J) ) 

HJPH ( I, J) =HMA+WW+ ( HA+ ( RJPH ( I , J ) -RC ) *PHI ) *DCOS ( NN* ( I - 1 ) *DA ) 
109  CONTINUE 

DO  111  KM=1, 1000 
ERR=0 . DO 

IF(KM.NE.l)  GO  TO  22 
DO  125  J=2 , JM1 
DO  125  1=1, IM1 
P(I, J)=PI+PD*(J-1)*DR 
POLD (I,J)=P(I,J) 

Q( I, J)=P(I, J)**2 
QIP1 ( I , J)=Q( I , J ) 

QIM1 (I,J)=Q(I,J) 

QJM1 ( I , J)=(PI+PD* ( J-2 ) *DR) **2 
QJP1 ( I , J)=(PI+PD*J*DR) **2 
125  CONTINUE 

DO  135  J=l, JMAX 
DO  135  1=1, IM1 
P(I,1)=PI 
P(I, JMAX)=PO 
Q( I , 1)=PI**2 
Q  ( I ,  JMAX)=PO**2 
135  CONTINUE 

22  CONTINUE 

DO  150  J=2 , JM1 

NEQ=IM1 

M=2 

LMAX=2*M+1 

DO  200  1=1, IM1 

IF( I .GE. IM)  GO  TO  23 

K=1-(1-I)*M 

GO  TO  21 

23  IF(I.GT.IM)  GO  TO  24 
K=IM1 

GO  TO  21 

24  K=( I MAX- I ) *M 

21  ARE A=RR ( I , J ) *DA  *DR 

A1 ( I , J)=- (HIPH( I , J) **3/RIPH( I , J) +9IMH( I , J) **3/RIMH( I , J) ) /DA 
*-CON*(HIPH(I, J)*RIPH(I, J)-HIMH(I, J)*RIMH(I, J))/DSQRT( 

*Q( I, J) )-VALUE*H( I, J)*AREA/DSQRT(Q(I, J) ) 

A1IP1 ( I , J)=HIPH( I , J)**3/RI?H( I, J)/DA-RIPH( I , J) *CON*KIPH 
* ( I , J)/DSQRT(QIP1( I , J) ) 

A1IM1  ( I ,  J)=HIMH(  I ,  J)  +  *3/RI?*H(  I ,  J)/DA+CON*HIMK(  I ,  J)+RIMH 
* ( I , J ) /DSQRT ( Q INI ( I , J ) ) 

R(K) =DA  /DR/DR* ( ( (RJ?H( I , J) I , J) **3-RJMH( I , J) * 

*HJMH( I,J)**3)*Q(I,J) ) -QJP1 ( I , J) *RJ?H( I , J) *HJPH( I , J) **3 
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*-QJMl ( I # J) *RJMH( I , J) I , J) **3 ) 
*-VALUE*H( I , J) *AREA*DSQRT(Q( I ,  J) ) 

DO  210  L=1 , LMAX 
210  Z(K, L)=0 .DO 

Z ( K , M+ 1 ) =A1 ( I , J ) 

IF(I.NE.l)  GO  TO  27 
Z(K,M+2)=A1IM1( I , J) 

Z(K, LMAX)=A1IP1( I, J) 

GO  TO  200 

27  IF(I.NE.IMl)  GO  TO  28 
Z(K,M)=A1IP1( I, J) 
Z(K,LMAX)=A1IM1(I, J) 

GO  TO  200 

28  IF(I.NE.IMM)  GO  TO  29 
Z(K, 1 )=A1IM1( I ,  J) 

Z(K,M+2)=A1IP1( I ,  J) 

GO  TO  200 

29  IF(I.NE.IM)  GO  TO  31 
Z(K, 1 )=A1IP1( I , J) 

Z(K,M)=A1IM1( I , J) 

GO  TO  200 

31  IF(I.LT.IM)  GO  TO  32 
Z ( K , LMAX ) =A1 IM1 ( I , J ) 

Z(K, 1 )=A1IP1 ( I , J) 

GO  TO  200 

32  Z(K, 1)=A1IM1(I, J) 

Z(K, LMAX)=A1IP1( I, J) 

200  CONTINUE 

CALL  GAUSS (Z,R,M,NEQ) 

DO  250  1=1, IM1 
IF(I.GE.IM)  GO  TO  25 
K=l- ( 1- I ) *M 
GO  TO  19 

25  IF(I.GT.IM)  GO  TO  26 
K=IM1 

GO  TO  19 

26  K=( IMAX-I ) *M 
19  Q(I,J)=R(K) 

P(I,J)=DSQRT(Q(I,J)) 

IF(I.NE.l)  GO  TO  33 
QIP1( I , J)=R(K+M) 

QIM1(I, J)=R(K+1) 

GO  TO  250 

33  IF( I.LE. IM)  GO  TO  34 
QIM1(I, J)=R(K+M) 

IF(I.NE.IMl)  GO  TO  18 
QIP1 ( I , J)=R( K-l ) 

GO  TO  250 

18  QI?1(I, J)=R(K-M) 

GO  TO  250 

34  IF(  I  .LT.  IM)  GO  TO  35 
IF( I.NE. IMM)  GO  TO  36 
QIP1 ( I , J ) =R ( K+l ) 
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QIM1( I , J)=R(K-M) 

GO  TO  250 

36  QIP1(I, J)=R(K-M) 

QIM1(I/ J)=R(K-1) 

GO  TO  250 

35  QIP1(I, J)=R(K+M) 

QIM1{ I , J)=R(K-M) 

250  CONTINUE 
150  CONTINUE 

DO  300  1=1, IM1 
M=1 

LMAX=2*M+1 

NEQ=JM2 

DO  350  J=2 , JM1 
K=J-1 

AREA=RR( I , J) *DA*DR 

A2 ( I , J)=-DA  /DR/DR* (RJPH( I , J) *HJPH( I , J) **3+RJMH( I , J) * 
*HJMH( I , J) **3 ) -VALUE*H( I , J) *AREA/DSQRT(Q( I , J) ) 

A2 JP1 ( I , J) =DA*RJPH( I , J ) *HJPH ( I , J ) * *3  /DR/DR 
A2JM1 ( I , J)=DA*RJMH( I , J) *HJMH( I , J) **3  /DR/DR 
ZR(K)=Q( I , J) * ( (HIPH( I , J) **3/RIPH( I , J) +HIMH( I , J) **3/RIMH( 
*1, J) )/DA- (HIMH( I , J) *RIMH( I , J) -HIPH( I , J) *RIPH( I , J) ) *CON/ 
*DSQRT(Q( I , J) ) )-QIPl(I, J)*(H:PH(I, J) **3/RIPH( I , J)/DA-CON* 
*HIPH( I , J) *RIPH( I, J)/DSQRT(QIP1( I, J) ) ) -QIMl ( I , J) * (HIMH( I , 
*J) **3/RIMH< I , J)/DA+CON*HIMH( I , J) *RIMH( I , J)/DSQRT(QIM1 ( I , J) 
* ) ) -VALUE*H( I , J) *AREA*DSQRT(Q( I , J) ) 

DO  310  L=1 , LMAX 
310  ZZ(K, L)=0.D0 

ZZ ( K , M+ 1 ) =A2 ( I , J ) 

IF( J . NE . 2 )  GO  TO  40 

ZZ (K, 1 )=0 .DO 

ZZ ( K, LMAX) =A2 JP1 ( I , J ) 

ZR(K )=ZR(K' - R2 JM1 ( I , J) *PI *PI 
GO  TO  350 

40  IF(J.NE.JMl)  GO  TO  42 
ZZ(K, 1)=A2JM1( I, J) 

ZZ ( K , LMAX ) =0 . DO 
ZR(K)=ZR(K) -A2 JP1 ( I , J) *PO*PO 
GO  TO  350 

42  ZZ(K, 1 )=A2 JM1 ( I , J) 

ZZ(K,LMAX)=A2JP1(I,J) 

350  CONTINUE 

CALL  ZGAUSS ( ZZ , ZR, M, NEQ ) 

DO  400  J=2 , JM1 
K=J-1 

QA(I,J)=ZR(K) 

PA( I , J)=DSQRT(QA( I , J) ) 

IF( J.NE.2)  GO  TO  43 
Q JP 1 ( I , J ) =ZR ( K+M ) 

Q JM1 (I,J)=PI*PI 
GO  TO  50 

43  IF(J.NE.JMl)  GO  TO  45 
QJ?1( I, J)=PO*PO 
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Q JM1 ( I ,  J ) =ZR ( K-M ) 

GO  TO  50 

45  Q JP 1 ( I , J ) =ZR ( K+M ) 

Q JM1 ( I ,  J ) =ZR ( K-M ) 

50  DEV=DABS(PA( I ,  J) -POLD( I ,  J)  ) 
Q<I,J)=QA(I,J) 

POLD( I , J)=PA( I ,  J) 

P( I , J)=PA( I ,  J) 

ERR=ERR+DEV 
400  CONTINUE 
300  CONTINUE 

IF(ERR.LT.AMAXER)  GO  TO  70 
111  CONTINUE 
70  CALL  LDTOT 
W=TOTLD 
NIT=KM 
RETURN 
END 
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To  find  the  minimum  wear  with  respect  to  the  same  radius 


SUBROUTINE  WEAR(XR  ) 

IMPLICIT  REAL *8  (A-H,0-2) 

DIMENSION  W( 18 ) 

COMMON  /BLOCKA/C , DA , DEM , DR , ETA , OMEG , P I , PO , RO , WE , IMAX , JMAX 

COMMON  /BLOCKD/AMAXER, RI , NN, IM1, JM1 

COMMON  /BLOCKB/AN, DT , HA , PHI , RC , WW , NIT 

DA=AN/IM1 

DO  100  1=1, IM1 

W( I )=(HA+ (XR-RC) *PHI ) *DCOS (NN* ( 1-1 ) *DA) 

100  CONTINUE 
SMALL=W ( 1 ) 

DO  200  1=1, IM1 

IF ( SMALL . LE . W ( I ) )  GO  TO  200 

SMALL=W ( I ) 

200  CONTINUE 

WW=-( SMALL) 

RETURN 


t 
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C 

C 

C  To  solve  simultaneous  equations  by  Gaussian  Elimination 
C  with  tri-column  matrix  coefficients. 

C 

C 

SUBROUTINE  GAUSS ( A, B, MBAND, NEQ) 

IMPLICIT  REAL *8  (A-H,0-Z) 

DIMENSION  A(90)/B(18) 

CALL  ESP IE 
NQ1=NEQ-1 
DO  10  1=1, NQ1 
11=1+1 
I2=I+MBAND 
DO  20  11=11,12 
J=I 

K= ( J- I I + MBAND ) *NEQ+ I I 
IF(A(K) .EQ.O.DO)  GO  TO  20 
KK= ( J- I + MBAND ) *NEQ+ I 
C=A(K)/A(KK) 

J1=J 

J2=J+MBAND 

IF( J2.GT.NEQ)  J2=NEQ 
DO  40  JJ=J1,J2 
KKK= ( J J- I I + MBAND ) *NEQ+ I I 
KKKK= ( J J- I + MBAND ) *NEQ+ I 
A ( KXK ) =A ( KKK ) -C* A ( KKKK > 

40  CONTINUE 

B ( I I ) =B (II) -C*B ( I ) 

20  CONTINUE 
10  CONTINUE 

K=MBAND*NEQ+NEQ 
B ( NEQ ) =B ( NEQ ) /A ( K ) 

DO  50  I 1=1 , NQ1 

I=NEQ-II 

J1=I+1 

J2=I+MBAND 

S=0.0 

IF( J2.GT.NEQ)  J2=NEQ 
DO  60  J=J1 , J2 
K=( J-I+MBAND) *NEQ+I 
60  S=S+A(K) *B( J) 

KK=MBAND  *N£Q+ 1 
50  B(I)=(B(I)-S)/A(KX) 

RETURN 

END 
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C  To  solve  simultaneous  equations  by  Gaussian  Elimination 
C  with  tri-diagonal  matrix  coefficients. 

C 

C 

SUBROUTINE  ZGAUSS ( A ,  B , MBAND , NEQ ) 

IMPLICIT  REAL *8  (A-H,0-Z) 

DIMENSION  A(54)/B(18) 

CALL  ESP IE 
NQ1=NEQ-1 
DO  10  I=1,NQ1 
11=1+1 
I2=I+MBAND 
DO  20  11=11,12 
J=I 

K= ( J- I I + MB AND ) *NEQ + I I 
IF(A(K) .EQ.O.DO)  GO  TO  20 
KK= ( J- I + MB AND ) *NEQ+ I 
C=A(K)/A(KK) 

J1=J 

J2=J+MBAND 
IF(J2.GT.NEQ)  J2=NEQ 
DO  40  JJ=J1,J2 
KKK= ( J J- I I + MB AND ) *NEQ+ I I 
KKKK= ( J J- I + MBAND ) *NEQ+ I 
40  A(KKK)=A(KKK)-C*A(KKKK) 

B ( I I )=B (II) -C*B ( I ) 

20  CONTINUE 
10  CONTINUE 

K=MB AND  *  NEQ +NEQ 
B ( NEQ ) =B ( NEQ ) /A ( K ) 

DO  50  11=1 , NQ1 
I=NEQ- I I 
J1=I+1 
J2=I +MBAND 
S=0 . DO 

IF( J2.GT.NEQ)  J2=NEQ 
DO  60  J=J1,J2 
K= ( J- I + MB AND ) *NEQ+ I 
60  S=S+A(K)*B( J) 

KK=MBAND*NEQ+I 
50  B(I)=(B(I)-S)/A(KK) 

RETURN 

END 


To  find  the  total  load  support  including  hydrodynamic, 
mechanical  and  shear  supports,  where, 

WMECH  is  the  total  mechanical  load  support 
TOTLD  is  the  total  load  support 
PHLS  is  the  percentage  of  fluid  load 
W  is  the  total  fluid  load  support 


SUBROUTINE  LDTOT 
IMPLICIT  REAL*8  (A-H,0-Z) 

COMMON  /BLOCKA/C , DA , DEM , DR , ETA , OMEG , P I , PO , RO , WE , IMAX , JMAX 
COMMON  /BLOCKC/H( 18,20),P(18,2Q),Q(18,20), RR( 18,20) 

COMMON  /BLOCKD/AMAXER , RI , NWAVE , IM1 , JM1 

COMMON  /B LOCKE/P F AC , PHLS, PM, TOTLD, WMECH 

W=0 . DO 

WMEC=0 . DO 

DO  100  J=2, JMAX 

DO  100  1=1, IM1 

JM=J-1 

HT=H(I, J)/C 

CALL  ASPRAT(HT,BI) 

IF(I.EQ.l)  GO  TO  10 
IF(I.EQ.IMl)  GO  TO  20 
I IM=I-1 
IIP=I+1 
GO  TO  30 
10  I IM=IM1 
IIP=I+1 
GO  TO  30 
20  IIM=I-1 
IIP=1 

30  W=W+ ( P ( I IP, J) +P ( I IP, JM) +P ( I IM, J) +P ( I IM, JM) )/4 . DO* ( RR( I , J) 

* -DR/2 . DO ) * D A * DR * B I 

WMEC=WMEC+PM*  ( 1 .  DO-B I )  *DA*DP.*  ( RR  ( I ,  J )  -DR/2  .  DO ) 

100  CONTINUE 

TOTLD= ( W+WMEC ) * NWAVE 
WMECH=WMEC* NWAVE 
PHLS=W*NWAVE/TOTLD* 100 . DO 
RETURN 
END 
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c 

c 

C  To  find  the  frictional  force  and  the  coefficient  of 
C  friction,  where, 

C  FF  is  the  total  viscous  force 

C  FA  is  the  mechanical  friction  force 

C  AMU  is  the  total  coefficient  of  friction 
C  TORQ  is  the  torque  of  the  driving  shaft 
C  FFORC  is  the  total  combined  friction  force 
C 

SUBROUT I NE  FR I CT ( P SBAR , RF , AMU , FFORC , FFA , TORQ ) 

IMPLICIT  REAL *8  (A-H,0-Z) 

COMMON  /BLOCKA/C , DA , DEM , DR , ETA , OMEG , P I , PO , RO , WE , IMAX , JMAX 
COMMON  /BLOCKC/H( 18,20),P(18,20),Q(18,20) , RR( 18 ,20) 

COMMON  /BLOCKF/FCONV , WSTAR , ROOC , NWAVE 
COMMON  /BLOCKG/DEL, EH, E10H, EH3 
FF=0 . DO 
FA=0 . DO 

DEL=RO *  OMEG*  ETA/C/P  SB AR 

IM1=IMAX-1 

DO  100  J=2 , JMAX 

DO  100  1=1, IM1 

R=RR ( I , J ) -DR/2 . DO 

HT=H( I , J)/C 

CALL  EXPT(HT) 

CALL  ASPRAT ( HT , B I ) 

JM=J-1 

IF(I.EQ.l)  GO  TO  10 
IF(I.EQ.IMl)  GO  TO  20 
I IM=I-1 
IIP=I+1 
GO  TO  30 
10  I IM=IM1 
IIP=I+1 
GO  TO  30 
20  I IM=I - 1 
IIP=1 

30  DP=(P( IIP, J)-P( IIM, J)+P(IIP, JM)-P( IIM, JM) )/4. DO/DA 
FF=FF+R*R*  <R*ETA*0MEG*E10H+ . 5D0/R*DP*EH) *BI/RF*DA*DR 
FA=FA+ ( 1 . D0-3I ) *PSBAR*R*R*DA*DR/RF 
100  CONTINUE 
FSUM=FF+FA 

AMU=FSUM/WSTAR*NWAVE 

FFA= FA* NWAVE 

FF0RC=FSUM*NWAVE 

TORQ=FSUM*RF*NWAVE 

RETURN 

END 
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onono 


To  find  the  expected  values  of  film  thickness 
where  E10H  is  the  expected  value  of  E(l/H) 

SUBROUTINE  EXPT(HT) 

IMPLICIT  REAL *8  (A-H,0-Z) 

COMMON  /BLOCKA/C , DA , DEM , DR , ETA , OMEG , PI , PO, RO, WE, IMAX, JMAX 
COMMON  /BLOCKG/DEL, EH, El OH, EH3 
EH=HT*C 

IE(EH.LT.O.DO)  EH=0.D0 
EH3=HT**3*C**3 
I F ( EH3 . LT . 0 . DO )  EH3=0.D0 
IF(HT.LE.l.DO)  GO  TO  25 
ElOH=l .DO/HT/C 
I F ( HT . GT . 4 . DO )  GO  TO  50 

E10H=( 1 . D0-HT*HT) **3*DLOG( ( HT+ 1 . DO ) / ( HT- 1 . DO ) ) 
E10H=(E10H+2 . D0*HT* ( 33 . D0+HT*HT* ( 15 . D0*HT*HT-40 . DO ) ) 

*/15 . DO ) *35 . DO/32 . DO/C 
GO  TO  50 
25  CONTINUE 

E10H=-405 . DO+HT* ( 60 . DO+147 . D0*HT) 

E10H=-55 . DO+HT* ( 132 . DO+HT* ( 345 . DO+HT* ( - 160 . D0+HT*E10H) ) ) 
E10H= ( E10H/60 . D0+ ( 1 . DO-HT*HT ) *  *3  *DLOG ( ( 1 . DO+HT )/DEL ) ) * 
*35. DO/32. DO/C 
50  RETURN 
END 


c 

c 

C  To  find  the  asperity  ratio  of  the  surface  where 
C  BI  is  the  fraction  of  seal  subject  to  fluid  pressure 
C 
C 

SUBROUTINE  ASPRAT(HT, BI ) 

IMPLICIT  REAL *8  (A-H,0-Z) 

H=HT 
BI=1 .DO 

IF(H.GE.l.DO)  GO  TO  10 

BI=( 16 .DO+H* (35 . D0+H*H* ( -35 . D0+H*H* (21 .DO- 5 .D0*H*H) ) ) ) 
*/32 . DO 
10  RETURN 
END 


********  The  end  of  the  whole  programs  ******** 


c 


TYPICAL  1WUT  PARAMETERS  AMU  OUTPUT  RESULTS 


RE0U1RED  LOAD  SUPPORT  *  0.962300  03  L6S 

MAV1NES&  AMPLITUDE  HA  *  0.  3641000-03  INCH 

TILT  PHI  •  0.300000-03  IN/IN 

PRESSURE  CONV.  CR1T.  AMaXER*  6.000  PA I 

LOAD  CONV.  OUTER.  MLiMiT  >3.00  La  A 


6AL  RATIO  *  0.730 

HALP  MAX  ROUGH  HEIGHT  C  »  0»2tt*00D-4>6 
RPN  •  1600.0 

HALF  TIME  STEP  OT  a  0.300000-4* 

PI  *  16.70  PAX 

PO  a  316.70  PA1 

RAOIUA  OF  CONTACT  POINT  a  0.1S233U  01  INCH 

oinensionles*  of  contact  radios  •  o.is 
balanced  kaoiua  a  i.96«s~  inch 

HOA  •  0.330000-06  INCH  HOb  a  0.2*300(1-06  INCH 

M0A  a  0.979210  03  l3A  mU»  a  0.963060  03  LbA 

HINDI  UN  FILH  THICKNESS  HOa  0.30 17  oO -06  INCH 
FLUID  LOAD  SUPPORT  M«E6  a  0.962330  03  LBS 

N1TA  a  619 
N1T8  a  662 
NXTN  a  636 

STIFFNESS  OF  FLUID  »— C. 116660  07  L3A/IN 

INNER  LEAXAGE  01N  ■  -0.12090  00  ACFN 

OUTER  LEAKAGE  UOUT  a  -0.36600-01  SCFM 

HAAS  FL  ERR  UR*  -0.66330-01  SCFM 

C06P  OF  FRXCT*  0.00006673 
FRICT  FORCE*  0.639*10-411 
PERCENT  OF  FLUID  LOAD  a  1GO.OUOOO 
TOTAL  HECh  CONTACT  Eh  EAR  FORCE  a  OM 
TOT  FRICTION  CUE  TO  ASPERITY  CONTACT  * 

TOTAL  TOR CUE  a  0.067673 


LdS 

LOS 

O.u  LSS 

IN -La 


INCH 

SECONDS 
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